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1.  A  Statement  of  the  Problem  Studied 


The  main  objective  of  this  research  is  to  apply  recent  advances  in 
kinematics,  utilizing  fundamental  models  of  dynamics,  coupled  with  experi¬ 
mental  observations  in  order  to  develop  methods  capable  of  predicting  the 
performance  as  well  as  motion  characteristics  of  certain  basic,  finite-dwell 
type  of  intermittent-motion  mechanisms  including  Geneva  mechanisms,  ratchets 
and  escapements,  which  are  useful  in  industry  and  military  applications. 

The  contract  was  originally  awarded  for  a  two-year  period.  Later,  it 
was  awarded  an  additional  year  to  extend  the  theory  developed  earlier  under 
this  contract  to  two  related  areas  of  mechanism  design  problems,  namely,  the 
inomentary-dwel 1  type  of  intermittent-motion  mechanisms  and  an  investigation 
on  the  workspace  of  a  robotic  manipulator  system. 

2.  A  Summary  of  Significant  Accomplishment 

As  it  was  mentioned  in  the  preceeding  paragraph,  the  project  involves 
an  investigation  of  the  dynamic  as  well  as  the  kinematic  characteristics  of 
mechanisms  of  the  following  three  types: 

•  Intermittent-motion  mechanisms  of  the  finite-dwell  type 

•  Intermittent-motion  mechanisms  of  the  moment a ry -dwe 1 1  type 

•  Robotic  Manipulators 

Results  obtained  from  this  project,  which  are  believed  to  be  new  and 
significant  pertaining  to  the  above  three  catagories,  are  described  next. 

(a)  Kinematic  and  Dynamic  Analysis  of  Intermittent-Motion 

Mechanisms  of  tne  finlte-'Dwel  1  type 

A  mathematical  model  of  impact  has  been  formulated  for  the 
intermittent-motion  mechanism.  The  model,  which  includes  clearance, 
damping,  mechanism  elasticity  and  material  compliance  is  basic  to  the 
determination  of  the  dynamical  response  such  as  force  amplification,  dynamic 
load  and  the  stresses  of  mechanisms  with  intermittent-motion . 

The  theoretical  model  developed  has  been  applied  to  several  practical 
mechanisms,  the  Geneva  mechanism,  the  ratchet  and  the  escapement  mechanism. 
These  are  popular  Indexing  devices  and  have  been  widely  used  in  many  pro¬ 
duction  machinery  and  automatic  weapon  systems.  A  computational  procedure 
incorporating  the  dynamic  model  with  the  finite-element  approach  Is  deve¬ 
loped  for  the  Geneva  mechanism.  Dynamic  load  can  be  calculated  and  Is  sub¬ 
sequently  used  for  the  estimate  of  stresses.  An  automatic  mesh  generator 
for  finite-element  approach  Is  presented  for  the  Geneva  wheel,  although  not 
simple,  due  to  the  complicated  geometry  of  the  wheel.  The  computer-aided 
procedure  makes  possible,  for  the  first  time,  to  accurately  predict  the 
performance  as  well  as  durability  of  an  intermittent-motion  mechanism  with 
respect  to  its  oprational  speed. 
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In  the  theoretical  development,  a  new  approach  in  the  modeling  of  sys¬ 
tem  damping  is  presented.  Instead  of  using  damping  ratio,  which  is  diffi¬ 
cult  to  estimate  accurately,  a  new  damping  function  is  introduced,  which 
characterizes  the  speed  and  load  dependent  nature  of  damping.  Two  types  of 
damping  functions  are  proposed  and  both  of  their  corresponding  damping 
forces  satisfy  the  expected  hysteresis  boundary  conditions.  A  comparative 
study  of  the  present  model  with  conventional  dynamic  models  is  performed. 

It  demonstrates  the  characteristics  of  the  proposed  model  and  its  useful¬ 
ness  for  the  study  of  the  dynamics  of  intermi ttent-motion  mechanism. 

Fxperimental  data  on  stresses  of  the  Geneva  mechanism  under  both  dyna¬ 
mic  and  static  loadings  have  been  obtained  using  MTS  machine  on  a  4.14“  O.n. 
Geneva  wheel  (AISl  2340  Steel).  Ream  stresses  at  the  tip  and  root  are 
determined  using  strain  gauges.  The  results  check  well  with  the  theoretical 
prediction  using  finite-element  methods. 

The  result  of  this  part  of  the  project  has  technological  appl ications . 
Since  most  automatic  weapon  systems  involve  intermittent-motion  mechanism  as 
their  main  indexing  device,  the  work  represented  here  may  be  used  to  predict 
the  dynamic  response  of  the  weapon  system  and,  ultimately,  used  for  a  dura¬ 
bility  and  life-cycle  analysis  of  the  system. 

( b )  Analysis  and  Synthesis  of  Intermittent-Motion  Mechani sms  of  the 

Momentary-Dwell  Type 

A  general  approach  for  the  kinematic  synthesis  of  the  momentary-dwel 1 
mechanism  has  been  presented.  The  approach  involves  using  the  solution  of 
first-order  dwell  criterion  as  an  initial  estimate  and  the  development  of  a 
computer-aided  procedure  to  subsequently  readjust  the  mechanism  proportions 
for  a  closer  match  to  other  dwell  criteria  by  optimization.  The  proportions 
thus  obtained  provide  an  approximation  of  the  high  orders  of  dwell  and,  for 
most  practical  purposes,  the  solution  is  useful  and  acceptable.  The  ap¬ 
proach  has  been  illustrated  by  two  examples,  a  six-bar  linkage  and  a 
chaln-and-coupler-1 ink  mechanism.  These  mechanisms  are  selected  because 
they  represent  the  general  characteristics  of  most  momentary-dwel 1  mecha¬ 
nisms  which  are  difficult  to  analyze  and  design.  It  has  been  shown  in  this 
investigation  that  on  the  design  of  these  momentary-dwel 1  mechanisms,  ana¬ 
lytical  approach  alone  is  not  possible,  and  a  joint  consideration  using  a 
numerical  optimization  technique,  therefore,  becomes  necessary. 

( c )  On  the  Kinematic  Character! sti cs  of  Robotic  Manipulators 


This  part  of  investigation  concerns  an  analytical  investigation  of  the 
kinematic  characteristics  of  robotic  manipulators.  It  includes  two  parts. 
The  first  part  deals  with  the  study  of  manipulators  constructed  with 
revolute  joints  in  series  and  having  open  kinematic  chain  configuration. 

The  second  part  treats  the  problem  of  a  maneuverable  platform  mechanism  with 
closed  kinematic  chain  configuration. 

In  the  first  part,  the  kinematic  relationship  between  consecutive  re¬ 
volute  joints  with  both  unlimited  and  limited  rotatabilities  are  formulated, 
and  subsequently,  the  manipulator  workspace  is  analytically  derived  and 
represented  by  a  set  of  recursive  equations.  Two  of  the  basic  character¬ 
istics  on  manipulator  workspace,  l.e.,  the  existence  of  hole  and  void,  are 


Investigated  in  great  detail.  Based  on  these  theories,  algorithms  and 
computer  program  for  quantitative  evaluation  of  the  volume  of  workspace  and 
for  outlining  the  boundary  of  the  workspace  are  developed. 
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Furthermore,  an  important  theorem  on  manipulator  workspace  is  pre¬ 
sented.  It  is  found  that  for  a  given  manipulator  structure,  the  ratio  of 
the  volume  of  the  workspace  to  the  cube  of  its  total  link  length  is  a 
constant.  This  theorem  leads,  therefore,  to  the  introduction  of  an  effec¬ 
tive  manipulator  performance  index.  Rased  on  this  index,  an  algorithm  on 
the  optimization  of  manipulator  workspace  has  been  developed.  It  is  shown 
that  for  a  popular  commercial  robot  configuration,  an  optimum  structure  can 
be  found  through  both  analytical  and  numerical  means.  An  interesting  com¬ 
parative  study  with  human  arm  proportions  is  also  provided. 

In  the  second  part,  the  number  synthesis  of  the  n-SPS  platform  mecha¬ 
nism  is  first  performed  (S  and  P  denote  the  spherical  and  prismatic  joints, 
respectively).  It  is  shown  that  the  6-SPS  mechanism  is  the  only  platform 
structure  which  has  six  controllable  degrees-of-freedom  and  can  be  used  as  a 
general  manipulating  device.  Analytical  formulation  for  the  kinematic 
analysis,  synthesis  and  control  of  the  mechanism  is  presented.  Cross- 
sections  of  the  workspace  on  three  primary  planes  are  constructed.  Conse¬ 
quently,  several  basic  kinematic  problems  are  Investigated.  These  are: 
the  analysis  of  a  planar  4-bar  with  adjustable  link  lengths;  the  rotata- 
bility  of  a  ball-and-socket  joint;  and  the  possible  interference  or  locking 
phenomenon  during  the  motion  of  the  platform.  An  algorithm  and  a  computer 
program  to  outline  these  cross-sections  have  been  developed. 

The  results  of  this  research,  it  is  believed,  have  contributed  toward  a 
basic  understanding  of  the  kinematic  characteristics  and  design  of  robotic 
manipulators. 
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On  The  Dynamics  of 
Intermittent-Motion  Mechanisms 

Part  I.  Dynamic  Model  and  Response 

This  paper  deals  with  a  baste  problem  regarding  intermittent-motion  mechanisms, 
namely,  how  to  formulate  a  predicative  model  for  the  study  of  the  dynamics  of 
these  mechanisms.  A  mathematical  model  is  developed  in  this  investigation.  The 
model,  which  includes  clearance,  damping,  material  compliance,  anti  mechanism 
elasticity,  is  basic  to  the  determination  of  the  dynamical  response  such  as  force 
amplification  and  motion  characteristics  r  mechanisms  with  intermittent  motion. 
.-1  new  approach  in  the  modeling  of  svso  am  ping  is  presented.  Instead  of  using 
damping  ratio,  which  is  difficult  to  estir  accurately,  a  new  damping  function  is 
introduced,  which  characterizes  the  spe <  'd  load  dependent  nature  of  damping. 
Two  types  of  damping  functions  are  p  >osed  and  both  of  their  corresponding 
damping  forces  satisfy  the  expected  resis  boundary  conditions,  i.e..  zero 
damping  force  at  zero  and  maximum  displacement  of  contact.  A  com¬ 

parative  study  of  the  present  model  w  *r  vent.  na!  dynamic  models  is  per¬ 
formed.  It  demonstrates  the  characters, and  the  usefulness  of  the  proposed 
mode!  for  the  study  of  the  dynamics  of  intermittent -motion  mechanisms. 


I  Introduction 

Intermittent-motion  mechanisms  play  an  important  role  in 
modern  technology  and  industry.  Ihcv  are  the  essential 
elements  in  many  machinery  and  automatic  weapon  systems. 
A  variety  of  intermittent-motion  mechanisms  have  been 
described  in  the  literature  [l|.  Basically,  there  arc  two  types: 
One  gives  a  finite  dwell,  such  as  the  standard  external  Geneva 
mechanism  |2),  the  other  provides  instantaneous  or 
momentary  dwell  (3|.  From  the  point  of  view  of  mechanism 
structure,  intermittent-motion  can  be  generated  either  from 
single  or  compound  mechanisms.  Geneva  mechanisms,  cam- 
follower  systems,  and  linkage  |4)  are  examples  of  the  single 
type;  compound  Geneva  mechanisms  (5,  6)  chained  linkage 
(7],  and  geared  linkages  [3),  belong  to  the  compound  in¬ 
termittent-motion  mechanisms. 

Of  particular  interest  and  concern  in  this  investigation  is  the 
study  of  the  dynamics  of  high-speed  intermittent-motion 
mechanisms  of  the  finite-dwcll  and  of  the  single  type  which 
are  more  difficult  to  analyze.  There  are  many  reasons  at¬ 
tributing  to  this  fact.  First,  the  mechanism  is  characterized  by 
inherent  shock  loading,  that  is,  the  presence  of  essentially 
discontinuous  forces,  masses,  velocities  and  potential 
energies;  and  second,  because  of  this  special  loading  con¬ 
dition,  as  well  as  the  inertia  loading  due  to  high  speed 
operation,  the  effects  of  material  compliance,  mechanism 
elasticity,  clearance  and  damping  have  to  be  taken  into 
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consideration.  Damping  which  has  an  important  role  on  the 
dynamics  of  the  mechanism  is  especially  difficult  to  model.  It 
is  generally  known  that  damping  depends  upon  load  as  well  as 
speed.  The  conventional  way  of  modeling  requires  an  accurate 
estimate  of  (he  damping  ratio  which  is  difficult  to  attain.  The 
problem  is  therefore  complicated  and  needs  special  attention. 

The  purpose  of  this  investigation  is  to  develop  methods 
capable  of  predicting  the  performance  as  well  as  the  motion 
characteristics  of  certain  basic  mechanisms  which  generate 
intermittent-motion  and  are  of  general  interest  in  practice.  In 
particular,  the  objectives  are:  I)  to  formulate  an  analytical 
model  for  a  certain  class  of  intermittent-motion  mechanisms, 
and  2)  to  present  a  new  damping  function  in  the  model,  which 
characterizes  system  damping  and  which  can  be  quantitatively 
estimated.  The  result  of  this  investigation,  it  is  hoped,  will  not 
only  provide  a  predictive  dynamic  model  and  a  computational 
package  for  the  dynamics  of  intermittent-motion 
mechanisms,  but  also  an  added  advantage  -  a  deeper  un¬ 
derstanding  of  the  dynamic  load  on  intermittent -motion 
mechanisms  and  a  physical  insight  to  the  dynamic  charac¬ 
teristics  of  these  mechanisms.  Such  an  investigation  is 
believed  to  be  trustworthy,  and  to  the  authors’  knowledge, 
has  not  been  reported  in  the  literature. 

We  begin  by  discussing  the  background  of  the  dynamics  of 
the  intermittent-motion  mechanism  and  its  relationship  with 
the  well-investigated  field  of  clearance  in  the  mechanical 
joints.  An  analytical  model  suitable  for  intermittent-motion 
mechanisms  is  formulated  and  its  solutions  are  derived.  A 
comparative  study  of  the  proposed  model,  with  the  con¬ 
ventional  dynamic  models,  is  performed,  flic  results, 
whenever  possible,  aie  ilh.siiatcd  by  graphs  which  may  offer 
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X, 


Fig  i  The  dynamic  model  lor  intermittent-motion  mechanisms 


some  physical  msighi  into  the  behavior  of  the  intermittent  - 
motion  mechanism.  Applications  of  the  present  model  are 
extended  to  practical  mechanisms  such  as  Geneva  mechanisms 
and  ratchets,  and  the  results  are  reported  in  Part  II  of  this 
investigation.  The  result  of  this  investigation,  it  is  believed, 
will  represent  a  contribution  towards  the  more  efficient  and 
economical  design  and  analysis  of  mechanical  systems  in¬ 
volving  high-speed  intermittent  motion. 


2  Mathematical  Model 

All  intermittent-motion  mechanisms  share  some  common 
characteristics,  first,  there  exists  generally  shock  or  impact 
loading  which  may  be  inherent  in  the  mechanism,  or  as  a 
result  of  external  loading.  In  the  dynamic  analysis  of  the 
mechanisms,  therefore,  not  only  the  kinematics  but  also  the 
load  and  driving  characteristics  have  to  be  taken  into  con¬ 
sideration.  Secondly,  clearance  or  backlash  is  unavoidable  in 
any  mechanical  joint  or  connection  which  makes  possible  a 
loss  ol  contact  between  the  members,  llus  may  lead  to 
subsequent  impacts  which  would  give  rise  to  vibration.  The 
study  of  the  dynamic  response  of  intermittent-motion 
mechanisms  at  high  speed  is  therefore  analogous  to  the  in¬ 
vestigation  of  ihe  dynamics  of  mechanical  systems  with 
clearances  -  a  field  which  has  become  important  and  active  in 
the  past  decade. 

Most  of  the  investigations  relating  to  mechanism  clearance 
involve  the  formulation  of  mathematical  models  of  impact  as 
a  basis  of  their  study.  Essentially  simple  systems  of  springs 
and  dashpots  wi.h  clearance  are  used  to  simulate  the  complex 
and  nonlinear  phenomena  of  mechanical  connections.  Some 
of  the  notable  models  were  presented  by  1. anger  (8j,  Johnson 
|9|,  kobrinskiy  and  Babu/ky  (I0|,  Dubowsky  and 
Freudenstein  (II),  Hunt  and  Crossley  (I2|,  Herbert  and 
McWhannell  |I3|.  etc.  A  useful  outline  ol  the  background  lo 
this  field  of  investigation  is  given  by  Dubowsky  and  Gardner 
[14).  Recently  Haines  [151  gave  a  comprehensive  review  of  the 
subjeei . 

In  this  study  a  mathematical  model  has  been  formulated  to 
investigate  the  dynamic  response  of  the  intermittent-motion 
mechanism.  The  physical  model  itself  is  basically  taken  from 
the  one  proposed  by  Dubowsky  and  Freudenstein  (1 1 1  on  their 
though) ful  study  of  clearance  in  the  mechanical  joints  The 
differences,  however,  lie  on  the  mathematical  modelin"  of  the 
stillness  and  damping  of  (he  system.  These  are  mechanism 
elasticity  and  most  importantly,  a  new  approach  which  in¬ 
volves  a  damping  function  introduced  in  the  modeling  of 
system  damping.  Two  types  of  damping  functions  arc 
proposed  and  their  corresponding  damping  forces  both  satisfy 
the  expected  hysteresis  boundary  conditions,  i.e.,  zero 
damping  force  at  zero  and  maximum  relative  displacement  of 
contact. 

The  physical  model  used  in  this  investigation  is  shown  in 
Fig.  I  The  constants  At,  and  A/j  represent  the  equivalent 
masses  of  the  driven  and  the  driving  elements,  respectively. 
For  instance,  they  may  correspond  to  the  Geneva  wheel  and 


pm  mi  ,i  iti  nn.i  me, Ii.iiiimii  A  .mil  /*  i cpi t'sci 1 1  ili<  ii  iu 
si  1 1 1  lies*,  .mil  (1. imping,  icspci  livdv .  .mil  t  ili-iinics  I  lie 
clearance.  I  lie  external  loading  applied  on  the  driving  clement 
is  icpicscilli’d 

Iii  this  investigation  the  system  stillness,  A  ,  is  modeled  by 
two  single  springs  in  series,  in  addmoii  to  the  suilacc  com 
pliuncc  (A  )  which  is  included  in  Dubowskv's  model  |ll|, 
another  spring  represented  by  the  mechanism  elasticity  (A  )  is 
considered,  following  Dubowsky  and  Freudenstein  |ll|.  the 
surface  compliance  is  modeled  using  the  Imeai  izcd  expression 
ol  Hertz  contact  theory.  For  instance,  iii  the  case  ol  a 
cylindrical  pin  contacting  with  a  Hat  surface.  " c  have  the 
deflection. 

|  Si i'l>  1 

o  -  (A  ,  *  A;)(/*/u)  In  -In/*  2  tl) 

1  (  A  ;  '  A  -  )  A  J 

where  the  subscripts  I  and  2  refer  to  the  parameters  associated 
with  the  elements  I  and  2,  or  the  driven  and  driving  members, 
respectively.  Parameters  u  and  h  arc  the  hall  length  ot  the  pin 
and  the  natural  base  The  linearized  stillness  ol  Hertz  con¬ 
tact,  A'.  .  represents  the  slope  on  the  load  and  deflection  curve 
characterizing  the  contact  phenomena  I  lie  system  stillness  is 
tllCliloiC. 

I 

=  C) 

ax) 

During  sontacl  stages  the  relationship,  represented  hv 
equation  (2).  although  it  is  nonlinear,  is  noi  loo  Far  away  irotu 
linear.  Therefore,  the  system  stiffness  has  csscntiallv  a 
piecewise  linear  characteristic  which  may  be  used  to  model 
complicated  nonlinear  bchavioi  ol  intermittent  motion 
mechanism.  Oil  the  other  hand,  equation  (2)  is  simple  enough 
hv  estimate  and  can  be  used  lo  verily  the  response  ot  the 
dynamic  sv  stem  with  simple  input,  such  as  harmonic  force. 

Tlte  damping  of  this  system  is  modeled  through  a  new 
damping  function  introduced  in  this  investigation,  as  /A  and 
will  represent  the  major  content  of  our  discussions  in  the  next 
few  sections.  The  equation  of  motion  for  the  model,  show  n  in 
l  ig.  1,  may  he  easily  w  ritten  as, 

MX,  +  UX,  +  K..X,  =/>(/)  (?) 

where  Af  represents  the  equivalent  mass  of  Af ,  and  Af ; ,  I)  and 
A't.  denote  the  damping  function  and  equivalent  system 
stiffness,  respectively,  and  A",  represents  the  relative 
displacement  between  A,  and  X,:-  A  simple  numerical  in¬ 
tegration  scheme,  for  instance,  the  Rungc-Kutta  scheme,  can 
be  used  to  yield  the  solution  of  equation  (3). 

In  the  modeling  of  damping,  it  is  convenient  to  model  it  by 
viscous  damping  using  a  damping  ratio  which  is  generally 
known  to  be  speed  and  load  dcpendcnl.  However,  an  accurate 
estimate  of  this  ratio  is  difficult.  This  is  apparent  in  many  ol 
the  previous  investigations  on  the  study  of  the  dynamics  ol  a 
mechanical  system,  such  as  mechanical  joints  [II). 

In  this  investigation,  one  of  our  purposes  is  lo  model 
damping  in  a  way  that  some  of  the  modeling  parametci  s  could 
be  estimated  from  either  empirical  or  published  cxpci  micntul 
data.  This  is  done  by  introducing  a  damping  function,  rather 
than  the  damping  ratio,  and  lo  express  it  in  terms  ol  the 
product  of  a  damping  factor  and  another  function.  Ihc 
former  can  be  derived  from  the  coefficients  of  restitution 
which  can  be  lound  experimentally  and  such  ■|ila  arc 
available  in  the  literature  (16).  The  latter  depends,  more  or 
less,  on  a  heuristic  choice  such  that  the  resulting  damping 
force  would  represent  the  required  dynamic  characteristics  ol 
the  system,  such  as  the  boundary  conditions  on  the  hysteresis 
loop.  Two  types  ol  damping  functions  are  presented  in  this 
study.  It  will  be  shown  later  lhal  ihc  modified  viscous 
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damping  approach  is  preferred  Wc  iK'gin  with  the  Itrsl 
dumping  t unction. 

3  Kirsl  Damping  Function,  />, 


where  /,  corresponds  to  the  time  at  the  moment  that  two 
In  idles  are  separating.  Kewriimg  equation  (X).  we  obtain. 


1  A.  i,  ’  i 

A  ! 

i  ( 1 

5A.  /.  '  v 

h  > 

1  .  (  /  ) 

'  12  37  > 

123/  ' 

1237  7 

1  et  l),  represent  the  first  damping  function,  we  then 
define, 

/),  =  <  ,  -  / ,  (-V.  )  (4) 

where  (',  represents  a  damping  factor  and  r,(.V, )  is  a 
displacement  function. 

3.1  Displacement  Function,  7',(.V.  I.  The  displacement 
function  is  essentially  a  surface  deflection  function.  Follow  ing 
the  approach  of  Hunt  and  Crossley  (121.  instead  of  using  a 
displacement  function  which  is  a  power  series,  we  simply 
choose  a  linear  function,  i.e., 

7  , (.V.  I  =  V,  ,  (5) 

The  function,  7, |.V, ),  represents  the  simplest  form  of  a 
surface  deflection  function  whose  corresponding  damping 
force  can  satisfy  the  expected  boundary  conditions  of  the  zero 
damping  force  at  .V,  -Oand.V,  =  .V,)lkW . 

3.2  I  he  Damping  Factor,  (',.  Hie  damping  factor,  C,,  is 
denied  ftom  llic  coefficient  of '  restitution,  h  lly  representing 
/>  as  a  polynomial  in  terms  of  the  relative  approaching  velocity 
o!  two  impact  bodies  obtained  from  the  experimental  data 
and  equating  it  to  the  calculated  value  from  the  equation  ot 
motion,  an  expression  of  C,  can  be  derived.  The  procedure  is 
outlined  in  the  following. 

1)  By  correlating  the  data  from  a  good  account  of  many 
experiments  (I6|  conducted  with  different  initial  velocities  of 
impacting  bodies,  we  cun  express  h  by  means  of  a  least 
squares  fit  as: 


h  =  i 


(6) 


faking  only  the  first  term  on  the  series,  we  then  have, 

h  =  <t,  -  u,  t 


(7) 


where  1  represents  the  relative  approaching  velocity  of  two 
impacting  bodies  and  u*  denotes  the  coefficients  of  the  ap¬ 
proximating  polynomial. 

2)  From  the  definition  of  classical  mechanics  the  coefficient 
of  restitution  can  be  calculated, 

7>=-.V,  ‘  (8) 

where  A'.,  represents  the  relative  departing  velocity,  and  it  can 
be  voh  ed  from  the  equation  of  motion. 

The  equation  of  motion,  equation  (3),  can  now  be  rewritten 
as, 

3/.V.  +  (C,.V  ).V,  +  K.  X  =  /'(()  (9) 

t  he  analytical  solution  of  equation  (9)  is  difficult  to  obtain 
due  to  the  nonlinear  term  associated  with  the  damping. 
However,  a  homogeneous  solution  of  equation  (9)  can  be 
obtained  using  the  polynomial  approximation,  i.e., 

..  ,  (<  I  +A,  ) 

A.ltlil.l- 


6  M 


i  y 


I4C,  I  .  +  A,  )((  ,  I  +K,  )  I 


/  +  .. 


1 203/: 

The  relative  velocity  alter  impact  may  then  be  obtained  as. 
It,  1+ A.  ) 

23/ 


(10) 


A',  (f,  )±l 


I  t, ; 


<4C,  r,  +  A,  )(<’,  I  ,  +  A,  M 
_  24  M: 


(III 


(12) 


l.qualing  the  coellieienls  ol  like  powers  ol  I  ,  liom  equations 
<7)and  1 1 2),  we  obtain. 


2...  3/ 

C,  =  _  A',  11.1) 


Equation  (13)  shows  that  the  damping  factor,  t'|,  is 
proportional  to  surface  stiffness,  and  the  coefficient  o,  which 
agrees  with  the  findings  reported  by  Hunt  and  C'rossley  (I2| 
uses  an  energy  balance  based  on  an  estimated  dissipated 
energy  of  impact  bodies. 

The  first  damping  function  can  therefore  be  represented  as 
/),  ±«).75<,,A'.  I  V.  (14) 

Because  ol  the  displacement  Itmction,  the  dettved  damping 
lunction,  equation  (14).  will  give  a  damping  (otce  satisfying 
the  requited  hysteresis  conditions. 


4  Second  Damping  Function,  f)2 

Viscous  damping  has  been  eommonlv  used  in  the  modeling 
of  system  damping,  such  as  Dubousky  and  Frcudenstein’s 
model  |ll).  I  his  approach,  although  convenient 
mathematically,  has  certain  limitations,  f  irst,  it  would  lead  to 
a  damping  force  which  is  maximum  at  the  beginning  of  the 
impact  process  because  of  the  use  of  a  constant  damping 
coefficient,  and  second,  this  would  give  a  negative  sum  of 
damping  and  spring  forces  during  the  return  stoke  of  the 
hysteresis  loop.  In  addition,  the  value  would  reach  its  peak 
when  the  relative  displacement  between  two  bodies  vanishes 
(Fig.  3).  These  situations  are  not  realistic  and  often  make  the 
viscous  model  ineffective.  It  is  therefore  desirable  to  develop  a 
new  approach  which  would  satisfy  the  hysteresis  conditions  of 
impact  as  well  as  preserve  the  basic  characteristics  of  viscous 
dumping. 

The  proposed  second  damping  function,  D:,  is  a  modified 
viscous  damping  factor.  It  is  the  product  of  a  damping  factor 
and  a  smoothing  function,  or  a  transition  function.  Instead  of 
treating  the  damping  factor  as  a  constant  throughout  the 
impact  process,  like  the  conventional  approach,  a  function 
called  the  transition  function,  denoted  as  7\,  is  introduced.  It 
is  essentially  a  smoothed  step  function.  When  it  is  multiplied 
by  the  damping  factor,  C;,  the  resulting  damping  function 
would  represent  the  behavior  of  system  damping,  especially  in 
the  initial  impact  period.  The  derivation  of  the  damping 
factor.  02.  follows  a  similar  approach  as  the  Inst  damping 
factor,  Ol . 

The  transition  function,  shown  in  Fig.  2,  is  defined  in  a 
transition  /one  specified  by  the  region  Os  A.  <>  where  the 
lower  and  upper  bounds  correspond  to  the  initial  and  final 
stages  of  the  impact  process.  There  may  be  various  ways  to 
define  the  transition  function,  one  way  to  represent  is: 

T:( A.)  =  'V'  1  exp[|(A,  -<)-  I.V,  .ll«//<)j  (15) 

where  the  parameter  <,  which  has  a  unit  of  displacement, 
defines  the  width  of  transition  /one,  and  q  is  a  parameter 
specifying  the  curve  path  within  the  transition  /one. 

Figure  2  represents  a  graphical  form  of  equation  (15).  The 
transition  function  is  zero  when  A",  sO.O  and  is  one  when 
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,V.  *i  I  hr  value  nl  >  may  he  uilmraiily  clmscii,  bin  it  must 
vadvty  (lie  following  condition',. 

O-.-.V,,,,,  (lb) 

and 

C:r;(A\  )X,  +  A..A.  20  (17) 

Equations  (15)  and  (17)  provide  (he  conditions  that  guarantee 
the  sum  of  the  damping  and  spring  for*  *o  be  positive,  i.c., 
the  hysteresis  loop  is  within  the  first  quadrant  (Tig.  3). 

4.1  Damping  Factor,  The  second  damping  factor  is 
supposed  to  have  the  same  characteristics  as  a  viscous 
damping  coefficient,  which  can  be  calculated  by  knowing  the 
homogenous  solution  of  equation  (3)  with  zero  transition 
zone,  and  then  following  the  same  procedures  as  the 
calculation  of  the  first  damping  factor.  It  gives 


C'. 


(In  />) - 
(In  b):  +  n: 


(18) 


Equation  (18)  can  be  used  to  estimate  the  damping  factor  of 
impacting  bodies  due  to  translation.  1  arger  coefficient  of 
restitution  means  lower  damping  factor  and  stiller  surface 
compliance  leads  to  a  larger  damping  factor  or  a  lower 
coefficient  of  test  it  in  ion.  The  second  damping  function  can 
be  written  as 


(.v  ) 


(19) 


The  new  damping  function,  /),,  gives  a  damping  force 
similar  to  the  one  predicted  by  Dubowsky  (1 1|,  except  in  the 
transition  zone.  The  damping  varies  exponentially  in  the 
transition  zone  and  it  is  invariant  elsewhere.  The 
corresponding  damping  force,  therefore,  is  zero  at  the 
beginning  and  the  end  of  impact;  and  it  reaches  maximum  or 
minimum  value  at  A \-t  depending  upon  whether  the 
damping  process  is  in  its  working  or  returning  stroke.  With 
this  approach,  the  resulting  hysteresis  boundary  conditions 
are  satisfied. 

the  major  dilleience  between  the  (wo  purposed  damping 
functions  is  that  the  first  damping  function  is  a  nonviscous 
model  while  the  second  is  a  viscous  model.  Consequently,  the 
damping  force  derived  from  the  first  function  gives  a  much 
thinner  hysteresis  loop  as  compared  to  the  second  and  it  is 
particularly  thin  during  the  high  speed  zone.  This  is  a 
drawback  and  w  ill  become  obv  ious  in  our  later  discussions. 


5  Dynamic  Load 

In  the  design  of  high-speed  intermittent-motion 
mechanisms,  stress  distribution  on  machine  parts  as  well  as 
the  contact  stress  between  two  impacting  parts  are  of  major 
concern.  In  order  to  evlauate  these  stresses,  the  load  must  be 
known,  but  the  determination  of  this  load  is  not  simple.  Due 
to  the  inherent  shock  loading  of  the  intermittent-motion 
mechanism  ami  other  factors  such  as  backlash,  clastic 
deformations  under  load,  imbalance,  load  and  the  inertias  of 
mechanism  components,  the  resulting  action  is  a  variable 
dynamic  load.  Not  many  investigations  on  the  dynamic  load 
of  intermittent-motion  mechanisms  are  on  the  record  117). 
Lee  1 1 7]  presented  a  simple  dynamic  model  using  the  classical 
mechanics  theory  to  predict  the  dynamic  load  of  a  Geneva 
mechanism,  and  showed  that  the  effect  of  the  dynamic  load  is 
significant,  especially  at  higher  speeds.  It  is  therefore  the 
prim?'  •  concern  in  this  investigation  to  apply  the  proposed 
model  in  order  to  predict  the  dynamic  load  of  a  general  class 
of  intermittent-motion  mechanisms. 

It  is  expected  that  the  result  would  represent  an  im¬ 
provement  over  Lee’s  [1 7)  since  most  of  the  dynamic 
characteristics  of  intermittent-motion  mechanisms  have  been 


considered  in  the  present  model.  The  result  of  that  com¬ 
parison  will  be  presented  in  118). 

The  dynamic  load,  P,h  is  defined  as, 

P,i  =  A',.  ‘X,  (20) 

The  spring  deflection,  X,,  is  the  solution  of  the  equation  of 
motion  (equation  (3))  and  can  be  obtained  numerically, 
depending  on  the  geometrical  parameters  of  the  mechanism, 
the  load  condition  P(t),  and  the  initial  conditions,  A",  (0)  and 
A-,  (0).  The  spring  deflection  is  also  an  implicit  function  ol 
many  factors  such  as  clearance,  equivalent  mass,  spring 
stiffness  and  damping,  etc. 

The  dynamic  load  ratio  is  a  useful  parameter  and  is  in- 
•roduced  and  defined  as  follows: 

Pr«  ~  (21) 

This  parameter  is  dimensionless  and  is  used  to  measure  the 
dynamic  effect  of  the  mechanism.  Almost  similar  in  function 
to  the  classical  "design  acceleration  factor”  (8),  which  is  a 
function  of  the  weight  of  the  part  and  its  natural  frequency  of 
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“starling  velocity"  concept  which  is  considered  to  be  more 
accurate  in  the  modeling  of  high-impact  shock  problems.  The 
model  is  assumed  to  he  accelerated  up  lo  some  finite  velocity 
in  a  negligible  length  of  time  and  then  maintains  that  velocity 
during  the  initial  and  most  important  phase  of  the  shock 
motion.  The  parameter,  Pr,  is  useful  in  the  design  of  in¬ 
termittent-motion  mechanisms. 

6  Comparative  Study 

The  model  proposed  in  this  investigation  has  certain 
characteristics  which  make  it  especially  suitable  for  the 
dynamic  study  of  intermittent-motion  mechanisms.  In  this 
section  a  comparative  study  of  the  present  model  with  some  of 
the  known  conventional  dynamic  models  is  presented.  The 
purpose  is  to  demonstrate  these  characteristics. 

t  hree  types  of  models  which  are  representative  and  well- 
known  are  chosen  for  comparative  study.  These  are: 
Dubowsky  and  Freudenstein’s  impact  pair  model  [11);  Hunt 
and  Crossley’s  model  [12]  and  Herbert  and  McWhannell’s 
model  1 1-1);  and  the  classical  model  of  Johnson  [9|.  Each  of 
these  selected  models  has  its  own  characteristics  and  merits. 
Our  investigation  owes  much  to  these  previous  investigators 
whose  results  have  been  an  inspiration  and  have  prompted 
this  study. 

Realizing  the  loading  condition  varies  with  different 
mechanisms,  for  comparison  purposes  a  harmonic  forcing 
function  is  assumed  for  all  models.  In  addition,  initial  con¬ 
ditions  arc  the  same  (zero  initial  displacement  and  an  initial 
velocity  of  25.4  cm  s).  We  are  concerned  essentially  with  the 
behavior  of  the  following  quantities  which  characterize  a 
dynamic  model.  These  are:  the  output  motion  characteristics; 
the  hysteresis  loop;  and  the  dynamic  load.  We  begin  with  the 
motion  characteristics. 

6.1  Motion  Characteristics.  Displacement  curve  of 
various  models  obtained  numerically  from  the  equation  of 
motion  are  plotted  in  Fig.  4.  They  represent  the  solution  of  a 
complete  cycle  of  a  sinusoidal  forcing  function  with  same 
initial  conditions.  We  have  the  following  conclusions: 

1.  Because  of  the  impact  forcing  function  and  the  initial 
conditions,  the  displacement  curves  of  various  models  have 
the  following  three  different  modes:  (</)  contact  mode  -  this  is 
found  on  the  first  half  period  of  the  driving  force  due  to  low 
initial  velocity,  and  consequently,  low  kinetic  energy,  the 
mass,  AT; ,  in  this  mode,  generally  making  contacts  with  the 
mass.  A/, .  (A)  free-flight  mode  -  it  occurs  when  (he  driving 
force  changens  its  direction.  In  this  mode,  the  mass  AT;  loses 
contact  with  the  mass  A/,  and  is  ready  lo  impact  on  the  op¬ 
posite  side  of  the  mass  AT, .  (i)  impact  mode  -  it  occurs  in  the 
second  half  of  the  period.  A  characteristic  of  this  mode  is  that 
the  impact  phenomena  dominate  and  there  are  frequent 
impacts  and  bounces.  This  is  due  to  the  high  kinetic  energy 
resulting  from  the  previous  free-flight  mode. 

2.  Besides  having  three  modes,  the  displacement  curve 
consists  also  of  two  components  of  frequency  response.  One 
component  corresponds  to  the  forcing  frequency  and  the 
other,  which  is  superimposed  on  the  first,  has  a  frequency 
relating  on  the  surface  compliance,  the  equivalent  mass  and 
the  damping  of  the  system. 

3.  A  comparative  study  of  various  displacement  curves,  as 
shown  in  Fig.  4,  demonstrates  an  interesting  remark.  The 
solution  obtained  from  using  the  second  damping  function 
(i.c..  the  viscous  model)  is  more  stable  than  the  ones  obtained 
by  nonviscous  models,  namely,  the  first  damping  function, 
and  the  models  of  Hunt  |I2|  and  Herbert  |I3|.  This  instability 
is  primarily  due  to  the  displacement -dependent  nature  of  the 
nonviscous  models  which  predict  a  low  damping  force  at  or 
near  the  high  speed  region  of  the  impact  process,  i.e.,  near 
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Fig.  4  Motion  characteristics  ol  various  mortals 
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Fig.  5  Hystarasis  diagram  using  the  first  damping  function  with 
different  damping  factor* 

A',  =  0.  This  shows  that  the  second  damping  function  is  more 
realistic  than  the  first  damping  function. 

6.2  Hysteresis  Diagram.  A  comparative  study  on 
hysteresis  diagram  is  of  importance  since  it  demonstrates  the 
damping  characteristics  of  the  dynamic  model.  The  following 
arc  the  conclusions  drawn  from  this  study: 

I.  The  first  damping  function  has  certain  characteristics: 
First,  the  hysteresis  loop  is  generally  thinner  than  the  one 
predicted  using  the  viscous  model.  This  implies  that  lower 
dissipation  energy  is  expected  from  the  model.  Second,  the 
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Fig.  6  Surface  hyitvrevl*  loops  predicted  by  different  damping 
functions 

hysteresis  diagram  is  thinner  in  the  high  speed  region  of  the 
impact  process,  i.e.,  near  the  origin,  compared  with 
elsewhere.  It  remains  to  be  thin  even  with  the  increasing 
damping  factor,  as  shown  in  Fig.  $.  This  is  not  realistic  and 
may  be  considered  to  be  a  drawback  of  the  first  damping 
function,  or  the  nonviscous  model  in  general. 

2.  The  second  damping  function  is  a  viscous  model  which 
relates  closely  to  Dubowsky's  model  except  in  the  transition 
zone.  The  required  boundary  conditions  are  now  satisfied, 

i.e.,  the  calculated  damping  force  vanishes  at  the  origin. 

3.  The  difference  between  the  first  and  the  second  damping 
functions  is  apparent  (Fig.  6).  The  former  is  not  a  viscous 
model  and  the  dissipation  energy  predicted  is  always  much 
lower  than  the  one  predicted  using  the  second  damping 
function. 

Another  way  to  see  the  difference  between  the  two  damping 
functions  is  to  consider  a  special  case  when  the  clearance,  r, 
vanishes  in  the  physical  model  and  ask  what  would  happen  to 
the  hysteresis  diagram.  As  r  approaches  zero,  the  nonviscous 
model  (the  first  damping  function)  would  present  a  number  of 
eight-shaped  hysteresis  diagrams  while  the  second  damping 
function  would  predict  the  expected  elliptical  shape  (Fig.  7) 
which  conforms  to  the  classical  Kelvin-Voigt  model  where 
there  is  no  impact  at  all. 

6.3  Dynamic  Load.  From  the  designer’s  point  of  view, 
information  on  the  dynamic  load  is  of  ultimate  concern  which 
leads  to  the  evaluation  of  stresses  on  machine  parts.  In  the 
following  a  comparative  study  on  dynamic  load  using  various 
models  is  given  and  the  comparison  is  on  three  areas:  the 
dynamic  load  profile,  the  effect  of  initial  conditions  on  the 
dynamic  loads,  and  extreme  values  of  the  dynamic  load.  We 
first  begin  with  the  comparison  of  the  dynamic  load  profile. 

6.3.1  Dynamic  Load  Profile.  The  dynamic  load  profile 
will  reveal  much  information  regarding  the  dynamics  of  the 
impact,  such  as  the  maximum  impact  force  or  the  dynamic 
load  h.  the  compressive  displacement  (X  =  Kr/F),  duration  of 
contact,  and  the  distribution  of  dynamic  force  during  impact. 
A  comparison  of  various  dynamic  load  profiles  is  given  in 
Fig.  8.  The  following  conclusive  remarks  may  be  drawn: 


FOKX 


Fig.  7(a)  with  clearance 


Fig.  7  Surtaca  hyatareaia  diagrams  o!  various  modals 


1.  Two  groups  of  curves  can  be  seen,  one  corresponds  to 
the  viscous  model  and  the  other  belongs  to  the  nonviscous 
model.  It  is  shown  that  the  first  damping  function  gives  a 
dynamic  load  profile  between  Hunt's  112)  and  Herbert’s  (I3J 
models,  and  all  give  a  larger  dynamic  force  compared  to  the 
viscous  model  during  most  of  the  contact.  This  is  due  to  the 
fact  that  less  energy  is  dissipated  and  much  is  being  used  to 
contribute  to  impact. 

2.  The  rate  of  change  of  the  dynamo  load  is  higher  at  the 
initial  period  of  impact  for  the  present  model  using  the  second 
damping  function  than  for  the  nonviscous  models. 

3.  The  contact  time  for  the  second  damping  function  is  less 
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IH'iinu-iilally.  As  nunc  cx|Kiimcnlal  data  become  .iv.til.iHc,  a 
picdicativc  dynamic  model  can  be  established. 

3.  Out  of  the  two  proposed  damping  functions,  the  second 
damping  I  unci  ion  (viscous  model)  which  involves  a  transition 
function  is  lavored.  It  oilers  solutions  which  are  more  viable 
and  realistic  as  compared  to  the  first  damping  function  or  the 
nonviscous  damping  model  in  general.  Moreover  the  second 
damping  function  would  predict  a  hysteresis  diagram  which 
conforms  with  the  classical  Kelvin-Voigt  model  of  nonimpact 
as  the  clearance  vanishes.  Basic  theory  and  model  presented 
here  will  be  applied  to  three  practical  uncimilicm  motion 
mechanisms,  i.e..  Geneva  mechanisms,  ratchets,  and 
escapements,  the  results  of  which  are  presented  in  Part  II. 


than  the  one  predicted  by  the  nonviscous  models.  This  is 
because  much  of  the  energy  is  dissipated,  and,  therefore,  the 
system  has  less  energy  to  do  work  in  order  to  maintain 
contact.  Separation  therefore  occurs  earlier. 

4.  The  dynamic  load  profile  of  the  second  damping  func¬ 
tion  approximates  most  closely  a  half  sine  except  at  the  initial 
phase  and  the  finishing  stage  of  contact. 
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6.3.2  Extreme  I  alues  of  the  Dynamic  Load.  It  is  in¬ 
teresting  to  examine  the  dynamic  load  predicted  by  various 
models  under  an  extreme  condition,  namely  the  condition  that 
one  of  the  impacting  body,  say,  body  I,  is  considered  to  be 
fixed,  i.e.,  A/,=oo.  This  is  essentially  the  infinite-mass 
assumption  used  by  Johnson  [9|.  The  purpose  is  to  give  an 
estimate  of  the  maximum  possible  dynamic  load  range.  The 
follow  ing  remarks  can  be  made  referring  to  Table  1 : 

1 .  Comparing  all  models,  Johnson's  model  gives  the 
highest  maximum  dynamic  load  range.  This  is  due  to  the  fact 
that  his  model  does  not  consider  damping. 

2.  The  maximum  dynamic  loads  predicted  by  the  reset  of 
the  models  fall  into  two  distinct  categories.  On  the  category  of 
viscous  models,  the  value  predicted  by  the  second  damping 
function  is  slightly  higher  than  Dubowsky's  model.  However, 
both  give  much  smaller  predictions  than  the  nonviscous 
models. 

7  Summary  of  Conclusions 

I  A  dynamic  model  which  is  especially  formulated  for 
intermittent-motion  mechanisms  is  presented.  Dynamic 
characteristics  of  this  model  are  demonstrated. 

2.  A  new  approach  to  model  the  system  damping  is 
presented  by  introducing  a  damping  function  rather  than  the 
conventional  damping  ratio  used  in  the  viscous  model.  Two 
types  of  damping  functions  arc  proposed:  One  represents  a 
nonviscous  model;  the  other  represents  a  viscous  model. 
There  are  many  advantages  of  the  present  approach:  f  irst, 
both  of  these  functions  predict  a  damping  force  which 
satisfies  the  expected  hysteresis  boundary  conditions.  Second 
and  more  important,  unlike  the  damping  ratio  which  is 
difficult  to  estimate,  these  functions  can  be  estimated 
systematically  through  a  procedure  outlined  in  this  in¬ 
vestigation.  They  are  derived  in  terms  of  physical  parameters, 
such  as  coefficient  of  restitution  which  can  be  obtained  ex- 
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The  theory  developed  in  Part  /  has  been  applied  to  the  determination  of  the 
dynamic  response  of  three  typical  types  of  intermittent-motion  mechanisms, 
namely,  Geneva  mechanisms,  ratchets,  and  escapements.  Insight  into  the  behavior 
of  such  mechanisms  can  be  obtained  by  studying  the  characteristics  of  the 
mathematical  model  and  by  numerical  experiments.  The  results,  whei  "ver  possible, 
are  illustrated  by  computer-plotted  graphs.  A  comparative  study  with  some  limited 
current  investigations  on  the  same  subject  is  provided. 


Possibly  the  most  significant  result  of  this  investigation  is 
that  a  dynamic  model  and  its  associated  new  approach  on 
system  damping,  proposed  in  Part  1,  will  exhibit  the  com¬ 
plicated  behavior  common  to  intermittent-motion 
mechanisms.  It  is  the  purpose  of  this  part  of  the  investigation 
to  demonstrate  this  finding.  In  the  following,  three  in¬ 
termittent-motion  mechanisms  -  Geneva  mechanisms,  rat¬ 
chets.  and  escapements  which  are  believed  to  be 
representative  and  of  practical  interest,  arc  chosen  for  study. 
In  each  of  them,  we  are  concerned  with  the  motion  and 
dynamic  characteristics  of  its  components  as  well  as  its 
dynamic  load.  Information  such  as  this  is  essential  for  design. 
We  begin  first  with  the  Geneva  mechanism. 

Geneva  Mechanism 

I  Introduction  to  Geneva  Mechanism.  The  Geneva 
mechanism  is  a  popular  indexing  device.  It  converts  a  uniform 
rotary  motion  to  an  intermittent  rotary  motion.  The  main 
advantage  of  the  mechanism  lies  in  its  simplicity.  However, 
the  Geneva  mechanism  has  two  undesirable  characteristics: 
First,  the  output  motion  starts  and  ends  with  nonzero  ac¬ 
celerations  which  usually  induce  shock  loading  and  can 
produce  vibration  in  the  driven  system;  and  secondly,  there 
exists  high  contact  stress  between  the  drive  pin  and  the  wheel 
slot,  which  can  cause  pin  wear  and  wheel  breakage.  Both  of 
these  factors  adversely  affect  performance  and  life,  and 
impose  severe  constraints  on  high-speed  operation.  Despite  its 
inherent  limitation,  the  Geneva  mechanism  continues  to  serve 
as  a  common  and  useful  indexing  device,  even  at  high  speed. 

There  are  several  reasons  to  choose  this  mechanism  for 
study.  First,  it  is  a  popular  indexing  device  which  has  been 
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widely  used  in  industry.  Second,  there  are  not  too  many  in¬ 
vestigations  of  the  dynamics  of  this  mechanism  on  the  record, 
except  in  a  very  few  cases  (1,  2].  The  mechanism  has  been 
operated  at  high  speeds  regardless  of  its  inherent  dynamic 
limitations,  and  as  a  result,  pin  wear  and  wheel  breakage  has 
often  occurred.  It  is,  therefore,  necessary  to  gain  deeper 
understanding  of  the  mechanism  and  to  obtain  some  design 
guidelines.  Consequently,  the  study  ol  the  dynamics  ot 
Geneva  mechanisms  for  high-speed  application  is  both  of 
important  and  of  practical  interest. 

The  Geneva  mechanism  design  has  been  the  subject  of 
investigation  for  many  years  [3J.  Much  of  the  work  concerns, 
however,  the  kinematics  of  the  mechanism.  This  is  mainly  due 
to  the  lack  of  a  mathematical  model  which  can  adequately 
represent  the  difficult  impact  phenomena  involved  in  in¬ 
termittent-motion  mechanisms.  Recently  Lee  (A),  using  classic 
mechanics  theory,  presented  a  set  of  performance  parameters 
and  applied  optimum  design  on  Geneva  mechanisms.  A 
simple  dynamic  model,  similar  to  the  one  used  on  gear 
dynamics,  was  first  introduced  to  estimate  the  Geneva 
dynamic  load.  In  this  paper,  the  model  developed  in  Pari  I  is 
applied  to  provide  a  more  predictive  dynamic  model  on  the 
dynamics  of  Geneva  mechanisms. 

2  Dynamic  Response  of  Geneva  Mechanisms.  This  section 
begins  with  a  description  of  the  kinematics  and  the  dynamic 
modeling  of  the  Geneva  mechanism.  Then,  numerical  ex¬ 
periments  on  the  application  of  the  dynamic  model  are 
presented  and  finally  a  discussion  is  given. 

Figure  I  shows  the  motion  graphs  of  a  typical  four-station 
Geneva  mechanism.  In  the  study  of  the  dynamics  ol  Geneva 
mechanisms,  the  model  of  Part  I  is  used.  Referring  lo  the  l  ig. 

I  of  Pari  I,  the  parameters  M,  and  AT  represent  in  this  ease 
the  eqnis  alenl  mass  of  the  ( ieneva  wheel  and  pin.  tespcclively. 
The  spring  coefficient.  A',.,  denotes  the  equivalent  stillness  ol 
the  material  compliance  and  the  mechanism  elasticity.  A 
simple  cantilever  beam  approach  which  was  employed  by  Lee 
(4)  was  adopted  again  to  represent  the  (ieneva  mechanism 
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Fig.  1  The  kinematics  ol  the  Geneva  mechanism 

elasticity.  The  modified  viscous  damping  approach  (i.e,  /):) 
presented  in  I’art  I  is  used  for  ihc  system  damping  function, 
D.  The  forcing  function  P  refers  to  the  applied  load  at  the 
center  of  the  Geneva  pin.  The  clearance  r  denotes  the 
clearance  between  the  Geneva  pin  and  wheel. 

In  the  numerical  experiments  of  the  dynamic  model,  the 
equation  of  motion  is  solved  starting  at  the  contact  mode  with 
a  given  pin  initial  velocity,  say,  2.5-1  em/s.  The  dynamics  of 
the  pin  and  Geneva  w  heel  as  well  as  the  dynamic  load  ratio  are 
obtained  and  given  in  the  following. 

2. 1  Pin  Dynamics.  The  pin  dynamics  can  be  illustrated 
from  a  time-displacement  plot.  Figure  2  shows  the  pin 
movement  which  consists  of  vibratory  effect  due  to  Geneva 
load.  During  the  working  stroke,  the  pin  is  usually  in  contact 
with  the  wheel  surface.  During  the  return  stroke,  however, 
there  is  an  observable  period  of  free-flight  mode  (i.e.,  non¬ 
contact).  This  is  due  to  the  high  impact  velocity  introduced  by 
the  clearance  and  the  Geneva  static  load. 

Certain  observations  can  be  made  from  Figure  2. 

1  The  pin  in  general  follows  the  behavior  of  the  applied 
load,  P{t).  The  motion  of  the  pin,  as  shown  from  the 
displacement  plot,  characterizes  the  superposition  of  two 
effects,  namely,  the  applied  load  and  the  vibration  which  is 
due  to  stiffness  and  system  damping. 

2  The  pin-displacement  plots  are  characterized  by  a 
successive  sequence  of  free-flight  and  contact  modes  at  the 
beginning  of  both  the  working  and  the  return  strokes. 

(a)  The  free-flight  mode  is  caused  by  the  bounce  of  the  pin 
which  depends  on  the  impaeting  velocity  and  the  static  load 
applied  upon  it.  During  the  free-flight  mode,  the  pin  is  driven 
by  the  load  P  only,  therefore,  the  pin  velocity  increases  with 
an  increasing  slope  until  it  reaches  a  position  where  the  load 
reaches  maximum. 

(b)  At  contact  mode,  the  maximum  relative  displacement, 


iclmn  viioke 

3  The  relative  velocity  curve  shows  that  the  velocity 
changes  abruptly  at  the  initial  stage  of  both  the  working 
stroke  and  the  return  stroke,  and  is  then  followed  by  a  stable 
region  of  contact  mode. 

2.2  M  heel  Dynamics.  The  dynamics  of  the  Geneva  wheel 
are  shown  in  Fig.  3.  Dynamic  effects  can  be  seen  from  the 
figure,  which  cause  certain  distortions  of  the  kinematic 
profiles  of  wheel  displacement  and  velocity  (Fig.  I).  During 
the  noncontact  period,  the  wheel  separates  from  the  pin,  and 
moves  essentially  at  a  constant  velocity  due  to  inertia.  In  the 
contact  region,  acceleration  remains  almost  constant  while  the 
velocity  of  the  wheel  changes  abruptly.  Clearance  has  an 
influence  on  the  separation  of  the  pin  and  the  wheel.  Con¬ 
sequently.  it  has  a  dominating  effect  on  the  dynamics  of  the 
Geneva  wheel. 

2.3  Geneva  Dynamic  Load  Ratio.  In  this  section,  the 
dynamic  load  ratio  of  the  Geneva  mechanism  is  investigated. 
The  discussion  includes:  the  effects  of  pin  velocity,  surface 
compliance,  and  clearance  on  the  dynamic  load  ratio.  The 
following  gives  a  detailed  discussion. 

1  Using  different  initial  pin  velocities,  the  calculated 
dynamic  loads  are  shown  in  Fig.  4(u).  It  shows  that  the 
dynamic  load  is  affected  by  the  initial  pin  velocity  on  the 
working  stroke  of  the  Geneva  motion,  while  on  the  return 
stroke  the  effect  is  not  apparent.  This  is  because  of  the  viscous 
damping  function  which  can  effectively  absorb  the  impact. 
Therefore,  the  starting  velocity  of  the  working  stroke  has 
essentially  no  effect  on  the  return  stroke. 

2  The  effect  of  surface  compliance  is  shown  in  Fig.  4(h). 
Two  different  values  of  stiffness  are  used,  one  using  a  surface 
compliance  of  3502.5  x  10"*  N/m.  the  other  using  a  surface 
compliance  of  7005.0 x  |04  N/m.  The  graph  shows  that  the 
surface  compliance  has  a  dominant  effect  on  the  return  stroke 
of  Geneva  motion.  This  implies  that  the  dynamic  load  ol  a 
softer  surface  is  generally  lower  than  the  dynamic  load  of  the 
stiffer  surface,  and  this  is  of  course  obvious  from  intuition. 

3  Figure  4(e)  shows  the  effect  of  clearance  on  the  dynamic 
load  ratio.  It  has  an  important  influence  on  the  return  stroke 
of  the  motion,  but  has  no  effect  on  the  working  stroke,  and 
this  is  reasonable.  Moreover,  it  is  of  interest  to  note  that  in 
increasing  clearance,  the  peak  of  the  dynamic  load  not  only 
increases,  but  also  delays  on  the  return  stroke.  This  is  because 
the  clearance  allows  the  pin  to  increase  the  kinetic  energy. 
Therefore,  the  larger  the  clearance,  Ihc  larger  the  dynamic 
load  on  the  return  stroke. 

From  the  foregoing  discussions,  several  conclusions  can  be 
made:  1)  The  dynamic  load  between  the  pin  and  the  wheel  is 
considered  to  be  an  important  design  factor,  which  depends 
on  many  variables  such  as  initial  conditions,  surface  com¬ 
pliances  and  clearance.  2)  Clearance  has  a  dominating  role  on 
the  dynamic  load.  Therefore,  both  controlling  clearance  and 
using  a  material  with  lower  surface  compliance  are  ad¬ 
vantageous. 

2.4  Phase-Plane  Plot.  Figure  5  contains  the  phase-plane 
plot  of  the  motion  of  the  Geneva  mechanism  and  is  a  result  of 
numerical  calculation.  The  phase-plane  analysis  presents  a 
convenient  manner  of  displaying  and  interpreting  the  tran¬ 
sient  behavior  of  the  nonlinear  system.  The  following  ob¬ 
servations  can  be  made: 

1  In  the  initial  stage  of  the  contact  mode,  the  pm  generally 
exhibits  bouncing.  The  bouncing  depends  upon  the  initial 
velocity  of  impact  and  the  driving  load. 

2  The  driving  load  plays  an  important  role  on  the  dynamic 
responses  of  the  Geneva  mechanism.  Referring  to  the  phase- 
plane  plot,  the  graph  is  unsymmctric  with  respect  to  the  line  of 
zero  velocity  and  is  shifted  toward  the  direction  ot  the  driv  mg 
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FI*.  2  Th«  dynamic*  of  tho  Geneva  mechanism:  pin  displacement  and 
velocity 
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FI*.  3  The  dynamic*  of  the  Geneva  mechanism  wheel  dtaplacamant 
and  velocity 


force.  Furthermore,  the  slopes  of  the  curve  at  the  two 
boundary  points  H  and  C  differ  from  each  other.  The  phase- 
plane  curve  above  the  zero  velocity  line  refers  to  the  outward 
stage  of  contact,  while  the  lower  one  refers  to  the  inward  stage 
of  the  contact.  At  the  moment  the  pin  impacts  on  the  surface 
of  the  wheel  (point  B),  its  velocity  reduces  from  1 .4  to  0  m 's. 
during  the  inward  stage  of  contact,  which  converts  the  impact 
kinetic  energy  to  potential  energy.  During  the  outward  stage, 
the  relative  velocity  increases  from  0.0  to  its  maximum 
velocity  of  0.75  m/s  at  point  C,  but  decreases  to  0.56  m/s  at 
the  boundary  point  D.  This  is  due  to  the  fact  that  the  diiving 
force  acts  in  the  opposite  direction  of  the  pin  motion  which  is 
at  the  verge  of  bouncing  off  the  surface. 

3  There  are  two  limit  cycles  of  the  phase  plane  plot  due  to 
the  clearance.  The  size  of  the  limit  cycle  depends  on  many 
factors,  such  as  the  impact  velocity,  driving  force,  clearance, 
and  damping  system,  etc.  It  is  of  interest  to  see  the  limit  cycle 
shown  in  Fig.  5.  The  graph  shows  a  bounded,  quickly  damped 
limit  cycle.  This  implies  the  Geneva  mechanism,  under  proper 
design  and  operation,  has  small  uncontrolled,  but  bounded 
vibrational  characteristics  between  the  pin  and  the  wheel.  This 
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unstable  vibration  can  be  quickly  stabilized  after  a  few 
bounces.  Therefore  the  system  is  a  stable  one  globally. 


Ratchets 

3  Introduction  to  Ratchets.  The  ratchet  is  a  simple 
mechanism  and  is  also  a  versatile  device.  A  schematic  diagram 
of  a  ratchet  mechanism  is  shown  in  Fig.  6.  The  mechanism  has 
many  advantages,  such  as  simplicity,  low  cost,  and  reliability, 
and  it  is  also  noted  for  its  ability  to  carry  a  large  load  in 
relation  to  its  si/c.  But  there  arc  disadvantages  also.  It  is 
mainly  an  impacting  mechanism  and  therefore  has  very  poor 
dynamic  characteristics.  Although  there  were  ways  to  reduce 
the  impacts  in  certain  versions,  impact  will  almost  always  be 
present  to  some  extent,  and  this  can  lead  to  noise,  wear, 
control,  and  stability  problems  unless  the  rest  of  the  system  is 
properly  designed.  A  major  problem  lies  in  the  fact  that  the 
forces  resulting  from  impacts  are  amplified  throughout 
various  members  of  the  mechanism,  which  may  be  well  in 
excess  of  the  subsequent  drive  forces.  The  dynamics  of  rat¬ 
chets  are  therefore  important  and  warrant  our  investigation. 
In  this  section,  the  model  developed  in  Part  1  is  applied  to  a 
ratchet  in  order  to  predict  its  performance  and  motion 
characteristics.  To  the  authors’  knowledge,  a  comprehensive 
study  on  the  dynamics  of  ratchets  has  not  been  performed  and 
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reported  in  literature,  although  some  kinematics  on  ratchets 
are  available  (6].  This  investigation,  therefore,  represents  a 
first  attempt  to  treat  the  dynamics  of  ratchets.  We  begin  by 
discussing  the  basic  kinematics  of  the  ratchet  mechanism 
using  an  analytical  approach. 

4  The  Kinematics  of  Ratchets.  The  kinematic  motion  of  a 
ratchet  mechanism  generally  can  be  classified  into  two  modes: 
(I)  the  free-flight  mode,  and  (2)  the  contact  mode.  Each  of 
these  modes  is  discussed  in  the  following. 

4.  t  f'ree-FNght  Mode.  The  free-flight  mode  occurs  when 
the  pawl  does  not  make  contact  with  the  ratchet  wheel  (Fig. 
6 (ft)),  and  consequently  does  not  drive  the  wheel.  Using 
complex  numbers,  we  can  write  the  closure  equation  as 

ce"*  +  be"*  +  he"1  -  ae*  +  d  =  0  (11 

with  h  ss,  and  0  =  0„  +  ^  as  the  constraint  conditions. 
Rewriting  equation  ( I ),  we  have, 

ccosttt  +  bcosa  +  hcos0-acos^  +  d=O  'l 

r  (2) 

csin<l+ hsinar+ hsin0  -  asintj/ =  0  ) 

Expressing  in  terms  of  h, 

h2  =  a2  +  b2  +  c2  +  d2  +  2bcos(<t>  -  a)  -  2accos(0  -  <//) 

-  2abcos(^  -  o)  +  2rf|ccosC>  +  bcosa  -  acos^]  (3) 

The  parameter  h  is  a  distance  parameter  which  can  be  used 
to  define  the  region  of  various  modes. 

Referring  to  Fig.  6(b)  and  equation  (3),  if  the  parameter  h 
equals  teto,  the  contact  mode  starts.  At  this  mode  the  wheel 
contacts  with  the  pawl,  and  therefore  drives  the  wheel.  When 
h  is  larger  or  equal  to  5,  the  pawl  shifts  to  the  next  ratchet 
tooth. 


4.2  Contact  Motto.  When  the  pawl  comes  into  contact 
with  the  ratchet  wheel,  the  contact  mode  starts.  Referring  to 
Fig.  6(b) ,  the  closure  equation  of  this  mode  can  be  written  as: 

ce1*  +  be"' ~  oe  *  +  d  =  0  (4) 


Rewriting  equation  (4)  into  scalar  equations, 

ccos<t>  +  bcosa  -  ocos^  +  d  =  0  3 

csin<t>  +  hsinn  -osin^  =  0  J 

Rearranging  equation  (5),  we  obtain: 
a2  -  b!  +  c2  +  d2  -  2aecos«i  -  ^)  +  2r/[cos<t>  -  ocosi^l  =  0 


(5) 

(b> 


«/(t„)=vi„  (9) 

integrating  equation  (7)  gives  the  angular  displacement  of  the 
wheel. 

,  T  -  ■ 

'!'=  if  +  ^,.(f-r„)  +C„  (101 


The  equation  (6)  is  equivalent  to  the  displacement  equation 
of  a  planar  four-bar  linkage.  However,  instead  of  being  a 
constant  rotation,  the  input  motion  for  ratchets  undergoes 
mostly  a  nonuniform  oscillating  motion,  liquations  (I)  and 
(4)  or  equations  (3)  and  (6)  represent  the  motion  charac¬ 
teristics  of  the  ratchet  wheel.  They  are  useful  in  our  study. 

5  Various  Models  of  Ratchet  Mechanism.  5.1  The 

Classic  Model.  In  this  section,  the  dynamics  of  ratchets  are 
studied  using  classic  mechanics  approach.  It  is  based  on  the 
assumption  that  the  high-order  effect  caused  by  the  elasticity 
of  individual  members  is  neglected.  Two  types  of  modes  are 
considered:  the  contact  mode,  and  the  free-flight  mode. 

5.1.1  The  Free-flight  Mode.  During  the  free-flight 
motion,  the  wheel  moves  independently  of  the  pawl. 
Therefore  as  a  constant  torque  T  acts  on  the  wheel,  the 
equation  of  motion  of  the  wheel  is, 

i  =  T/lw  (7) 

where  /„  denotes  the  wheel  inertia. 

The  applicable  initial  conditions  are: 

i(t„)=i„  (8) 

ami 


5.1.2  The  Contact  Mode.  The  equation  of  motion  of  the 
ratchet  wheel  can  be  expressed  as 

/„v£k  -  «x  P„n  -  (  Ti.)  (II) 

where  n  and  k  are  unit  vectors  and  a  represents  the  moment 
arm  vector. 

The  pawl  load,  P„,  is  the  normal  force  which  acts  on  the 
wheel  and  the  tangential  force  is  neglected.  Equation  (ID 
gives  the  equation  of  motion  of  the  ratchet  wheel,  where  the 
pawl  force,  P„,  is  assumed  to  be  known.  The  input  force  P„ 
may  be  generally  expressed  in  the  following  form: 

P„=A  +  2  (/usinmr  +  gmcoswr)  (121 

m  -  I 

where  K,/m,  and  gm  are  Fourier  scries  coefficients. 

The  solution  of  equation  (II),  may,  therefore,  be  expressed 

as: 

2  »  * 

<1/  -  — •  [o/f  -  71  +  C|  /  +  (':  -  57  (  cosmA 

2/„  *  „  m _  |  \trr  m-  / 

where  C,  and  C:  are  constants  depending  on  the  initial 
conditions. 

Figure  7  shows  the  angular  displacement  (equation  1 1),  and 
the  plot  of  the  aiigulat  velocity,  rf,  versus  tunc. 
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5.2  The  Dynamic  Model.  The  dynamic  model  presented  in 
Part  1  is  applied  to  ratchets.  However,  some  physical  com¬ 
ponents  of  the  model  need  to  be  modified  in  order  to  fit  the 
special  requirement  of  ratchets.  Again,  Mg.  I  of  Part  I  can  be 
used  as  the  ratchet  model  for  this  investigation.  In  this  model, 
M |  and  M:  represent  the  equivalent  mass  of  the  wheel  and  the 
pawl,  respectively;  K,.  represents  the  equivalent  stiffness  due 
to  surface  compliance  and  mechanism  elasticity  of  the  wheel 
and  the  pawl;  D  represents  the  damping  function  which  is 
used  to  model  the  system  damping;  P„  denotes  the  driving 
force  that  the  pawl  exerts  to  the  wheel,  and  the  constant  load, 
Pw,  refers  to  the  load  applied  on  the  wheel  in  order  to 
maintain  contact  between  the  pawl  and  the  wheel.  The  driving 
force,  P„  can  be  expressed  as  a  Fourier  series  depending  on 
the  speciality  of  the  ratchet.  The  equation  of  motion  of  the 
contact  mode  can  therefore  be  expressed  as 


MX.  +  DX,  +  K,X,  =  T  1 P.  -  [  ~  ~~  P„  (U> 

L  M,  +  Mi  J  l  Mi  +  .Vf, 


where  X ,  represents  the  relative  displacement  between  the 
driving  pawl  and  the  wheel. 

During  the  period  of  noncontact,  the  equation  of  motion 
can  be  expressed  as; 

•V/,  X,  =  /'.  1 


(15) 


M:X:=P„ 


A  numerical  integration  scheme,  such  as  the  Rungc-Kutta 
method  can  be  used  to  yield  solution.  The  following  sections 
give  a  discussion  of  the  dynamic  responses  of  the  ratchet. 


A  Dynamic  Responses  of  Ratchets.  The  motion  curves  of 
the  ratchet  wheel  vary  widely  depending  upon  the  means  used 
to  power  the  drive  pawl,  the  input  force  to  drive  the  wheel 


Fig.  8  Motion  curves  ot  the  ratchet  wheel  under  various  loads 


can  be  a  spring,  a  solenoid,  a  cam  or  other  pulse  drive,  and 
therefore  each  corresponding  output  motion  of  the  wheel  will 
be  different.  Some  input  forces  may  give  a  high  impact  force 
to  the  wheel  and  thus  cause  losing  its  control,  therefore. 
certain  additional  design  refinements  are  required.  In  this 
section,  output  motions  of  the  ratchet  wheel  predicted  by  the 
dynamic  model  are  presented.  The  wheel  dynamic  load  is  also 
estimated. 

6  1  The  Ratchet  Wheel  Output  Motion.  For  numerical 
experiments  the  following  conditions  are  used:  1)  A  periodic 
forcing  function,  P„  =  444.8  (sin(80.0t  -  1 .0)  1 .0),  is  used  to 

drive  the  system;  2)  a  viscous  damping  function  based  on  a 
given  coefficient  of  restitution  of  0.5  is  used  to  represent  the 
system  damping;  and  J)  different  wheel  loads,  ,  are  used  to 
study  their  effects  on  the  wheel  output  motion.  Based  on  these 
conditions,  the  following  observations  can  be  made  from  Fig. 
8  which  provides  the  time-displacement  and  the  time-velocity 
plots  of  the  ratchet  w  heel. 

1  The  curves  of  the  wheel  displacement  and  velocity  are 
piecewise  nonlinear.  This  is  due  to  the  tree-flight  mode 
associated  with  the  contact  mode. 

2  The  dynamic  effects  of  the  model  can  be  seen  from  the 
wheel  displacement  and  velocity  curves,  which  cause  certain 
distortions  of  the  kinematic  profiles  of  the  w  heel  displacement 
and  velocity  (Fig.  7). 

3  The  wheel  velocity  increases  rapidly  at  the  contact  mode 
and  decreases  linearly  at  the  free-night  mode  referring  to  Fig. 
9.  Moreover,  a  small  amount  of  velocity  decrease  is  observed 
before  the  free-night  mode  starts.  This  is  due  to  the  effect  of 
the  wheel  load,  P „ . 

4  The  wheel  load,  ,  has  certain  effects  on  the  duration 
of  Irec-night,  the  maximum  wheel  velocity  during  contact, 
and  the  slope  of  the  wheel  velocity  during  the  f re*'- flight. 

5  The  wheel  tends  to  remain  in  contact  wm>  the  pawl  at 
the  final  stage  of  the  working  stiokc.  I  hen-fore,  this  con¬ 
tributes  lo  a  stable  inn  casing  ot  wheel  velocity . 
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Fig.  S  Th*  r»tch*l  dynamic  loads:  (a)  the  altect  ol  pawl  Initial  velocity: 
and  (h|  the  effect  ol  wheel  load 


From  ihe  foregoing  observaiions,  it  is  found  that  the  wheel 
motion  is  generally  instable  at  ihe  beginning  stage  of  the 
working  stage  and  is  stable  at  the  middle  and  final  stages  of 
the  working  stroke. 

6.2  Ratchet  Wheel  Dynamic  Load.  The  dynamic  load  of 
ratchet  between  wheel  and  pawl  is  important  and  is  of 
practical  concern.  Our  discussion  here  includes:  1)  ihe  effect 
of  initial  velocity  on  the  dynamic  load  of  ratchet  mechanism, 
and  2)  the  effect  of  the  wheel  load,  ,  on  the  dynamic  load 
of  Ihe  ratchet  mechanism.  Each  of  these  effects  is  discussed  in 
the  following. 

1  Using  different  pawl  initial  velocities,  the  dynamic  loads 
calculated  numerically  from  equation  (7)  are  shown  in  Fig. 
9(d).  It  is  found  that  the  dynamic  load  is  affected  significantly 
by  the  pawl’s  initial  velocity  at  the  initial  stage  of  the  working 
stroke.  The  higher  the  pawl’s  initial  velocity,  the  higher  the 
dynamic  load  as  well  as  the  longer  the  duration  of  the  I'rcc- 
flight  mode.  Therefore  it  may  be  advantageous  if  the  initial 
velocity  of  the  pawl  is  set  low  in  order  to  minimize  the 
dynamic  load  on  the  ratchet  wheel. 

2  The  ratchet  wheel  dynamic  loads  predicted  from  dif¬ 
ferent  wheel  loads,  /*„  arc  plotted  in  l  ig.  ')(/>).  It  is  in¬ 
teresting  to  find  that  the  higher  the  wheel  load,  the  shorter 
the  duration  of  the  frce-flight  mode  and  there  appears  to  be 
no  increase  on  the  dynamic  load.  Especially,  at  the  end  of  the 
working  stroke,  the  higher  the  wheel  load.  P, .  the  shorter  the 
duration  of  the  free-fiight  mode  and  a  lower  dynamic  load. 

6.3  The  Phase-Plane  Plot.  The  phase-plane  plot  of  the 
motion  of  the  ratchet  is  given  in  Fig.  10.  The  plot  starts  from 
the  point  .-I  of  the  contact  mode  with  a  pawl  initial  velocity 
0.25  m/s.  During  the  first  contact  stage,  the  pawl  speed 


REL  DISPLACEMENT. (CM) 
Fig.  10  Phita  plane  plot  ot  the  ratchet 


changes  from  0.25  to  -  0.12  m/s.  At  point  W.  a  free-fiight 
mode  is  followed  and  impact  then  occurs  at  point  C.  Finally, 
the  pawl  reaches  point  D.  From  the  phase-plane  plot,  certain 
conclusions  can  be  made. 

1  The  classic  definition  of  the  coefficient  of  restitution 
may  be  quantitatively  evaluated  using  the  phase-plane  plot. 
For  instance,  the  relative  velocity  at  points  .1  and  H  are  0.25 
m/s  and  -0.12  m/s,  respectively.  This,  therefore,  gives  a 
coefficient  of  restitution  of  0.5  for  the  impact  of  the  ratchet. 

2  The  curvature  of  the  phase-plane  plot  changes  rapidly 
within  the  contact  stage  near  point  B.  This  is  due  to  the  wheel 
load,  Pw,  which  emends  the  duration  of  contact.  Con¬ 
sequently,  the  pawl  speed  is  increased  before  the  separation  of 
pawl  and  wheel. 

3  In  the  free-fiight  mode,  the  pawl  speed  starts  from 
-0.12  m/s  and  increases  up  to  0.2  m  s  at  point  C'.  The 
change  of  the  pawl  speed  depends  upon:  I )  the  speed  at  point 
B,  and  2)  the  pawl  force,  P„ . 

4  The  relative  velocity  at  the  terminal  position  of  the  pin 
stroke,  point  D,  reaches  a  value  of  approximately  zero.  This 
demonstrates  that  the  damping  function  is  effective  and  it 
therefore  essentially  stabilizes  the  vibration  between  the  pawl 
and  the  wheel  at  the  final  stage  of  the  working  stroke. 

Escapement  Mechanisms 

7  Introduction  to  Escapemmls.  Escapement  is  a 
mechanism  or  device  which  is  designed  to  alternately  grab  and 
release  a  rotating  shaft,  wheel,  slide,  etc.  |6).  In  general,  it  can 
he  classified  into  three  types  according  to  Ms  functions  and  the 
manners  of  its  operation.  Ihese  types  are:  the  watch 
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escapement,  the  machine  escapement,  and  the  inverse 
escapement.  They  usually  belong  to  a  class  of  mechanisms 
designed  for  particular  purposes,  e.g..  a  lime  keeper,  or  as  a 
mechanical  controller  of  a  rotating  shaft,  or  as  a  ratchet 
performing  similar  or  better  functions.  Most  of  the 
escapements  are  very  simple  kinematically:  they  are  low-cost, 
rugged,  and  reliable.  The  mechanism  offers  many  uses  |6,  7). 
They  are  used  in  clocks  and  watches,  bomb  fuses,  typewriter 
composer  systems,  counters,  precision  stepping  drives,  etc. 

In  the  design  of  escapements,  certain  general  functional 
requirements  of  these  mechanisms  must  be  considered.  Some 
of  the  requirements  are:  1)  stability:  the  escapement  has  to  he 
fast  and  capable  of  operating  at  high  speeds.  For  instance  the 


Tl  ME  .(SEC) 

Fig.  1 2  Angular  displacement  ot  the  escape  wheel 


HIM  sclcclric  composer  requires  lour  Ihousandms  of  a  second 
to  reach  its  maximum  displacement  |x|.  Besides  being  fast,  u 
musf  be  positive  and  reliable,  and  have  minimum  rebound  at 
the  end  of  an  escapement  cycle.  2)  Minimum  impact  force;  in 
order  to  design  the  parts  for  adequate  strength,  the  extreme 
impact  force  must  be  determined;  and  3)  a  predictable  con¬ 
trolled  escapement  time. 

Since  the  function  of  the  escapement  depends  highly  on  the 
preciseness  of  its  motion,  usually  involving  high  speed 
operation,  therefore,  an  investigation  of  the  dynamics  of  the 
escapement  is  necessary.  There  are.  however,  not  many  in¬ 
vestigations  on  the  record  [8.  9|.  I.owen  and  lepper  |9|  did  a 
thoughtful,  and  believed  to  be  the  first  of  its  kind  in¬ 
vestigation  on  the  dynamics  of  the  pin  pallet  runaway 
escapement.  In  their  formulations,  rigid  body  and  momentum 
balance  approaches  are  used.  In  addition,  the  classic  coef¬ 
ficient  of  restitution  was  employed  to  formulate  the  impact 
phenomena.  It  is  thought  worthwhile,  in  this  investigation  to 
apply  the  general  dynamic  model  of  the  intermittent-motion 
mechanism  developed  in  Fart  I  to  the  escapement  mechanism 
In  our  model,  system  damping  and  material  compliance  are 
included.  A  comparative  study  with  the  results  reported  in  [91 
on  the  dynamic  responses  and  motion  characteristics  is 
performed. 

8  Dynamic  Responses  of  the  Escapement.  In  this  paper,  a 
watch  escapement  shown  in  Fig.  II  is  investigated.  The 
kinematics  of  the  mechanism  and  a  useful  computer  program 
have  been  developed  in  |4| 

The  dynamics  of  the  mechanism  involves  three  stages. 
These  are:  I)  The  coupled  motion:  In  this  phase,  the  escape 
wheel  is  driven  by  a  constant  torque  in  the  positive  direction 
of  rotation.  2)  The  free  motion:  This  stage  follows  the  end  of 
the  coupled  motion  or  when  separation  of  contact  occurs  alter 
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the  stage  of  impact.  In  (his  phase,  the  wheel  and  the  pallet 
move  independently  of  each  other  in  free  motion  where  the 
motion  of  the  pallet  depends  only  on  its  initial  conditions.  3) 
The  impact  stage:  This  stage  follows  free  motion  and  usually 
involves  a  reversal  ol  the  motion  of  the  pallet.  Under  certain 
circumstances,  temporary  motion  reversal  of  (he  wheel  is  also 
possible. 

In  this  investigation,  the  first  two  stages  of  (he  motion 
follow  essentially  the  formulations  of  |9|.  In  the  important 
stage  of  the  impact,  however,  our  dynamic  model  is  applied. 

Referring  to  Fig.  I  of  Part  I,  M,  and  M:  denote  (he 
equivalent  masses  of  the  escape  wheel  (the  driver)  and  the 
pallet,  respectively.  The  driving  Torce,  P,  is  calculated  from 
the  given  torque  with  variable  moment  arms.  Following  |9),  a 
coefficient  of  restitution  equal  to  0.3  is  used. 

The  initial  conditions  for  the  numerical  experiments  are 
specified  as  follows: 


relative  initial  displacement  =  0.0 
relative  initial  velocity  -  0.0 
initial  driving  lorcc  1.53  N 


>Thc  coupled  motion 
S  Stage 


relative  initial  displacement  0.0 
relative  initial  velocity  =  1.55  m  s 


The  impact  stage 


Dynamic  responses  of  both  the  pallet  and  the  wheel  are 
given  in  the  following  sections. 

8.1  Motion  Characteristics  of  the  Escapement.  F  igures  12 
and  13  give  the  output  angular  displacement  and  velocity 
curves,  respectively,  of  the  escape  wheel.  A  comparative  study 
with  the  result  of  Lowen  and  Teppcr  |8)  is  also  shown.  We 
discuss  separately  the  following  two  motion  Mages. 

I  The  Coupled  Motion 

The  displacement  predicted  by  this  model  and  the  one  in  |S| 
agree  well.  However,  if  the  relative  velocity  between  the  pallet 
and  the  wheel  increases,  dynamic  effects  due  to  system 
damping  as  well  as  material  compliance  will  he  significant.  In 
that  case,  the  deviation  between  the  two  apptoaches  would  he 
observable. 

The  velocity  of  the  wheel  is  found  to  be  oscillating  during 
the  initial  stage.  This  is  due  to  the  sudden  application  of  the 
driving  load  and  this  oscillation  dies  out  as  the  nine 
progresses.  Figure  13(h)  shows  an  enlargement  ol  the  com 


Journal  of  Machanlsms,  Transmissions,  and  Automation  In  Dsslgn 


SEPTEMBER  I'JHd,  Vol  U)'>/549 


10  Cb 


REL  DISPLACEMENT (  A*1' 


Fig.  t  S  Phase-plane  plots  ol  the  escapement  mechanism 


parison  of  the  result  of  this  investigation  with  that  of  (9). 
Resides  oscillations,  there  is  general  agreement  on  the  velocity 
curve  m  this  stage  ol  motion. 

2  I  he  Impact  Stage 

The  output  displacement  curve  shows  a  significant  decrease 
of  displacement  due  to  impact.  This  agrees  with  our  physical 
intuition.  The  time  duration  of  impact  is  calculated  to  be  80  ns 
which  accounts  for  about  3  percent  of  the  total  coupled 
motion  duration.  This  time  cannot  be  estimated  with  the 
approach  ol  [V|  and  consequently  it  was  neglected.  Figure 
13(f))  shows  there  is  a  sudden  oscillation  in  the  velocity  curve, 
and  a  decrease  of  velocity  as  a  result  of  the  impact 
phenomena. 

H.2  The  Escape  H  heel  Dynamic  Toad.  The  dynamic  load 
between  the  escape  wheel  and  the  pallet  has  a  significant  effect 
on  the  motion,  the  wear,  and  fatigue  of  the  parts.  Figure  14 
gives  the  dynamic  load  during  the  coupled  motion  as  well  as 
the  impact  motion.  The  following  conclusions  can  be  drawn. 

1)  The  dynamic  load  during  the  coupled  motion  is  found 
to  be  of  no  major  concern  since  it  appears  to  be  low.  The 
maximum  dynamic  load  is  only  1.78  N;  the  load  reaches  a 
final  value  of  1 .22  N. 

2)  I  lie  dynamic  load  dining  llic  impact  stage,  on  the  ollici 


hand,  is  significant.  The  maximum  value  is  '56.2  N  which 
would  be  the  main  source  causing  severe  vibrations  to  the 
wheel  and  the  pallet. 

S.3  The  Phase-Plane  Phi.  Figure  15  shows  the  phase- 
plane  plot  of  the  escapement.  We  have  the  following  ob¬ 
servations. 

1  During  the  coupled  motion,  the  graph  shows  that  the 
limit  cycle  is  small  and  vanishes  immediately  after  a  lew 
cycles.  No  bounce  has  been  observed. 

2  The  phase-plane  plot  of  the  impact  stage  is  unstable. 
This  is  because  the  relative  velocity  during  impact  is  large.  A 
significant  bounce  follows  immediately  after  the  impact. 

Conclusion 

The  dynamic  responses  of  Geneva  mechanisms,  ratchets, 
and  escapements  have  been  analyzed.  It  is  shown  that  the 
proposed  dynamic  model  and  its  associated  damping  function 
approach  exhibits  a  variety  of  dynamic  characteristics  which 
arc  believed  to  be  representative  of  intermittent-motion 
mechanism.  The  dynamic  load,  obtained  using  the  present 
approach,  is  an  important  factor  for  the  stress  analysis  and 
the  design  ol  these  mechanisms. 
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SYNOPSIS  This  paper  presents  a  general  approach  for  the  kinematic  synthesis  of  norentary-dwell  mech- 
anisms.  The  approach  involves  using  the  solution  of  first-order  dwell  criterion  as  an  initial  estimate 
and  the  developnc-nt  of  a  computer-aided  procedure  to  subsequently  readjust  the  mechanism  proportions 
for  a  closer  match  to  other  dwell  criteria  by  optimization.  The  proportions  thus  obtained  provide  an 
approximation  of  the  high  orders  of  dwell  and,  for  most  practical  purposes,  the  solution  is  useful  and 
acceptable.  An  example  of  a  six-bar  linkage  is  used  to  illustrate  this  approach. 


INTRODUCTION 

1.  Dwell  mechanisms  are  important  for  their 
many  uses  such  as  in  machine  tools,  packaging  and 
textile  machinery.  Basically,  there  are  two  types. 
One  gives  a  finite  owell,  such  as  the  standard 
external  Geneva  mechanism;  the  other  provides 
instantaneous  or  momentary  dwell.  The  latter  can 
be  used  for  performing  functions  while  machining 

is  in  motion,  for  example  in  flying  shears, 
labeling,  closing  of  cans,  etc.  The  fundamental 
difference  between  these  two  types  lies  in  the 
output  motion.  Mathematically,  the  displacement 
function  for  the  finite-dwell  mechanism  is  a 
constant  within  a  defined  dwell  region.  There¬ 
fore,  the  higher  order  derivatives  of  the  motion 
such  as  velocity,  acceleration,  and  shock  are 
zero  within  this  region.  In  addition  to  these, 
at  the  boundary  of  starting  or  ending  of  a 
finite-dwell,  the  motion  are  generally  discon¬ 
tinuous.  For  the  motion  of  the  momentary -dwell 
mechanism,  the  displacement  is  a  nonlinear 
♦unction  of  time  or  input  angle.  The  higher 
order  derivatives  of  motion  generally  do  not  equal 
to  zero,  and  equal  or  close  to  zero  only  at  an 
instant,  i.e.,  the  point  of  dwell.  Furthermore, 
the  motion  curves  are  continuous,  therefore  there 
is  no  exact  boundary  of  starting  or  ending  of  a 
m<v  ientary-dwel  1 .  Physically,  the  finite-dwell 
mechanism  is  stopped  to  have  dwell  by  a  locking 
device,  such  as  at  the  Geneva  mechanism  the 
locking  position  of  the  crank  can  stop  the  wheel 
from  moving.  On  the  contrary,  momenta ry-dwel  1 
mechanism  does  not  depend  on  any  locking  device 
to  have  dwell. 

2.  Because  of  this  basic  difference  on  motion 
between  the  two  types  of  intermittent -mot ion 
rxchanisns,  the  synthesis  nethods  on  these 
mechanisms,  therefore  differ.  In  the  synthesis 
of  finite-dwell  mechanisms,  much  emphasis  is  on 
the  minimization  of  shock,  I.e.,  on  the  irpreve- 

'  ent  of  the  kinematic  and  dynamic  characteristics 

in 


of  the  mechanism  for  longer  life  and  higher  speeds 
of  operation.  There  have  been  many  thoughtful 
investigations  in  this  respect,  and  a  summary  of 
them  can  be  found  in  [1,2]. 

3.  Of  particular  interest  in  this  investigation 
is  on  the  kinematic  synthesis  of  the  moment a ry- 
dwell  type  of  mechanisms.  The  purpose  is  to 
develop  analytical  and  computer-aided  procedure 

on  the  synthesis  of  mechanisms  with  prescribed 
dwell  characteristics.  The  importance  is  on  the 
design  of  dwell  mechanisms  which  can  possess 
higher  orders  of  dwell  at  one  or  more  instants 
during  the  motion  cycle.  The  development  of 
design  methods  in  this  regard,  however,  has  not 
been  easy.  One  primary  reason  is  that  the  motion 
of  momentary-dwel 1  mechanism  is  usually  an 
algebraic  equation  of  high  order,  having  coeffi- 
cents  which  are  nonlinear  functions  of  mechanism 
parameters.  Dwell  occurs  due  to  the  cancellation 
of  relative  motions  of  one  part  of  the  mechanism 
with  respect  to  the  other.  Since  the  momentary- 
dwell  mechanism  is  usually  made  of  compound 
mechanisms  (i.e.,  the  combination  of  two  or  more 
simple  mechanisms)  whose  notion  characteristics 
are  generally  not  readily  known,  the  derivation 
of  an  explicit  motion  expression  in  terms  of 
design  parameters  and,  subsequently,  finding  its 
solution  can  be  difficult. 

4.  There  have  teen  several  investigations  on 
the  kinematic  synthesis  of  rorentary  dv/ell  mech¬ 
anisms  and  most  of  them  involving  analytical 
treatment  [3,4].  As  the  complexity  of  the  mech¬ 
anism  grows  higher,  algebraic  methods  of  synthe¬ 
sis  become  more  difficult  to  develop  in  order  to 
satisfy  precisely  the  dwell  characteristics.  On 
the  other  hand,  for  most  practical  purposes, 
high  order  dwell  criteria  need  not  to  be  satis¬ 
fied  exactly  and  they  can  be  met  approximately. 

The  following  represents  an  attempt  to  develop 
an  approach  for  the  approximate  s>nthesis  of 
monentary-dwell  mechanisms. 
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A  general  aporouch  for  the  kinc  atic 
synth-sis  of  o  er.i  »ry-dUell  mechanisms  is 
Ores*?’’ ted.  The  method  involves  using  the  solu- 
ion  of  first-order  dwell  criterion  as  an  initial 
.  sti-tte  «*»*  i  the  Jevelop-ent  of  a  computer-aided 
procedure  to  subse ;  ientTy  readjust  the  mechanism 
nrc£-:.r tiers  for  a  closer  patch  to  other  dwell 
ritfia  by  optimization.  The  proportions  thus 
btainel  preside  an  approximation  to  high  orde^ 
dwe’ls  and,  for  cost  practical  purposes,  the 
solution  is  useful  and  acceptable. 

The  oo-;ute>--aidel  synthesis  orocedure  is 
based  on  an  elective  optimization  algorithm, 
“a-ely,  the  HCT  code  (Heuristic  Optimization 
zchnig.e  of  Lee  and  fre^denstein  [5])  to  reach 
.  solution. 

7  To  illustrate  this  basic  approach,  the 
neatic  synthesis  of  a  six-bar  dwell  linkage 
s  presented,  in  the  following,  we  begin  with  an 
analytical  investigation  on  the  motion  charac¬ 
teristics  of  the  rechanism.  The  kinematic 
githesis  of  the  mechanism  with  prescribed  dwell 
•oaracteristics  follows.  It  is  found  that  the 
first-order  dwell  criterion  is  a  second-order 
"quation  in  terms  of  the  half-tangent  of  the  out- 
it  angle  and  its  solution  can  be  readily  obtained 
...aly tical ly.  As  for  the  second-  and  high-order 
dwells,  nurerical  solution  involving  optimization 
•"cores  essential.  This  example,  therefore, 
•-onstrates  the  potential  usefulness  of  the 
.Kproach. 


9.  The  motion  of  the  pivot  point  bet.-.e- ■>  t*o 
couolers  BCE  and  EF  is  of  major  con-.e'-n.  The 
couplers  point,  E,  traces  an  essentially  o/al- 
shaped  path.  This  path  generally  has  high-order 
curvature,  and  therefore,  its  cori  espor  iir.g  radius 
is  not  a  constant.  However,  if  diversion-,  of  the 
mechanisms  are  so  designed  that  the  radius  of  tre 
curvature  of  the  point  E  is  equal  or  clo,-  to  a 
constant  during  certain  ra<’-va  of  the  -s.*ion  cxcle, 
and  if  the  link  length  rg  is  -’ade  eq_.fl  to  this 
instantaneous  radius  of  curvature,  then  t*e  ojtnut 
crar-x,  rg,  has  mcr enta-'y  dwell.  In  the  following, 
we  first  begin  with  the  '-otion  of  a  fc-r-bar 
couoler  point. 

10.  Referring  to  rig,  2,  the  position  of  the 
coupler  point  E  can  be  expressed  as 

nE  *  ( 1 -' ) r B  *  \rc  (1) 

where 

rg  -  1SI,  rg  s  JT5,  r^  -  Ac  and  '  -  e’J 


c 


Fig.  2 


"•’splace-ent  Analysis 

...  A  six-bar  dwell  linkage  [6],  shown  in  its 
schematic  view  in  Fig.  1,  consists  of  six  links 
-•'noted  as  rj ,  r2,  rg,  r$,  r$,  and  rg.  The  in- 
t  is  usually  a  continuous  rotary  motion  which 
..  applied  to  link  r^  and  -Ji  refers  to  as  the 
input  displacement  angle,  the  output  shaft, 
r  ntered  at  the  pivot  G,  would  then  provide  a 
:  c-entary  dwel  1  •t'otion  if  the  linkage  is 
p. uteri/  proportioned. 


Fig.  1 

Ten  design  para  eters  are  required.  These  are 
six  link  lengths;  the  coupler  bent  angle,  o’; 
t  •  extended  coupler  length,  a;  and  the  position 
p  a  »>N>rs  n(  thw  output  crank ,  namely,  the 
Inclination  tnjle  ird  the  fixed  distance  between 
the  rocker  pivot  a-d  the  outwit  pivot,  8  and  h. 


Rearranging  eg.  (1),  we  have  the  position  of  the 
coupler  point  E  in  the  Cartesian  coordinate  system, 


X 

X 


E 

E 


XB  + 


YB  + 


COS  n(Xc-XB) 

cos  »(Yc-Yb) 


sin,(Yc-YB) 

sin,(Xc-XB) 

(2) 


where 


XB  -  r2  cos  e, 
YB  =  r2  Sin  '1 


Xc  =  r4  cos  a2  ♦  r1 

YC  1  r4  sin  -’2 


(3) 


The  velocity,  acceleration  and  shock  equations  of 
the  coupler  point  E  can  be  easily  derived  by 
differentiating  equations  (?)  and  (3). 

11.  Two  forms  of  the  di spl ac*--  en t  equation  of 
the  six-bar  1 inkage  can  be  derj.od.  One  which  is 
derived  from  an  inverted  mechanism  (Fig.  3)  is 
explicit  in  terms  of  the  coxirdinsf.es  of  the  floating 
pivot,  F,  of  the  outlet  crank;  the  Other  is  an 
angular  displacement  equation  relating  the  input 
and  output  angular  relationship.  We  give  the 
results  of  derivation  as  follows; 


Vf  *  Vf  5  co 


(4) 


where 
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Fig.  3 


the  first-order  1~t.ll  's  o'  *  < 

velocity  ratio  which  r  , .■••U  tc  i-  *'0  at  an  iostnrt 
and  the  sat.isfact  ion  V  e  hi •.  •  *  atf  *.-nt  eq.  (5), 

i  .e. , 

A'  sin  ;  ♦  E‘  cos-.  ♦  C*  -*  0  (E) 

A  sin  ;  r  B  c  os ,  *  C  -  0  (5) 


-  V 

A  -  V>G 

-  V 

E  1  VXG  , 

( AA  *  *  EE 1 ) 

-(A?‘E2 

r  - 

r6 

2r6 

For  the  first-order  dwell,  eq.  (8)  and  (9)  should 
t>e  satisfied  simultaneously.  The  solution  is: 


Bq  *  BtRjSin^-uJ-RjSin^-a)] 
C0  = 


A  sin  4*  +  B  cos  *  ♦  C  =  0 


B  =  XE-XG 

C  -  -(A2+B?+r2-r2)/(2r6) 

Equation  (5)  can  also  be  rewritten  in  terms  of  the 
half- tangent  form  of  * 


1  A+ZA^+B  -Cr 


4=2  tan  [- 


■  2  tan'' 


These  two  solutions  of  eq.  (6)  correspond  to  two 
different  configurations  of  a  six-bar  linkage.  It 
is  worth  noting  that  the  displacement  equations, 
eq.  (A)  and  eq.  (5)  have  similar  form.  For  the 
study  of  high-order  motion  analysis,  the  latter  is 
chosen,  since  algebraically,  it  is  more  convenient 
to  handle. 

Dwell  Criteria 

12.  In  a  momentary  dwell  mechanism,  the  dwell 
is  defined  when  the  first  or  second  derivative  of 
the  output  motion  with  respect  to  the  input  notion 
vanishes  at  an  instant.  The  order  of  the  conseu- 
tive  derivatives  which  are  zero  at  that  instant 
refers  to  the  degree  or  order  of  dwell.  In  the 
following  the  dwell  criteria  for  the  first-order 
is  derived.  Criteria  for  higher  orders  can  be 
derived  in  a  similar  manner.  The  definition  of 


tin  ,  .  C'B-B'C 

*D  '  «fV-«r  (10) 
A'C-C'A 
cos  'n  =  AE'-BA' 

For  the  first-order  dwell,  eq.  (10)  should  always 
be  satisfied.  It  can  tie  seen  that  the  motion 
characteristics  of  the  coupler  point,  E  and  link 
lenoth  r,-  and  r,  play  i-portant  roles  in  the  first- 
order  dwPll, 

KINEMATIC  SYNTHESIS 

13.  In  the  previous  sections,  notion  equations 
and  criteria  of  dwell  have  been  developed.  These 
equations  are  nonlinear  and  need  to  be  solved 
simultaneously.  Analytical  solution  is  possible 
only  for  the  velocity  dwell  and  generally  is  not 
possible  for  hiqher  order  dwell.  Therefore  the 
kinematic  synthesis  with  prescribed  high-order 
dwell  characteristics  becomes  a  rather  tedious  and 
involved  task  in  which  numerical  approach  is 
essential.  The  synthesis  procedure  of  the  six-bar 
linlage  with  prescribed  dwell  involves  the 
determination  of  ten  design  parameters,  and  the 
satisfaction  of  specified  dwell  conditions.  These 
desiqn  parameters  include  six  link  lengths, 

(rl,  »"2»  h3»  r4.  r5 ,  and  rg),  two  extended  coupler 
lengths  (a  and  b),  one  coupler  anole  (a1),  and  one 
position  anqle  (p). 

14.  Generally,  in  the  rase  of  first-,  second-, 
and  third-order  dwells,  there  are  fewer  equations 
than  the  number  of  parameters  to  be  determined. 

For  instance,  for  the  third-order  dwell,  there 
are  four  conditions  including  one  displacement 
equation  and  three  dwell  equations.  Additional 
specifications  from  the  desinner,  therefore,  are 
necessary.  In  the  following,  we  first  introduce 
a  rotatability  criterion  of  the  six-bar  linkaoe, 
which  will  define  the  feasible  range  of  design 
parameters  of  the  mechanism.  Then,  a  method  of 
analytical  synthesis  for  the  first-order  dwell  is 
introduced.  Two  design  parameters  (rg  and  rg)  can 
be  specified  throughout  thr  synthesis  procedure. 

For  hiqher  order  dwells,  solution  of  dwell  equa- 


tions  explicit  in  design  paraneters  is  difficult 
to  obtain.  Therefore,  a  numerical  synthesis  proce¬ 
dure  up  to  the  third-order  dwell  is  presented. 

The  procedure  corbines  a  method  of  analytical 
synthesis  with  an  optimum  design  technique,  i.e., 
the  HOT  algorithm  [5],  In  the  following  section 
we  begin  with  a  discussion  on  the  rotatability 
criteria  of  the  six-bar  dwe’l  linkage. 

Rotatability  Criterion 

15.  The  rotatability  study  of  a  six-bar  linkage 
uses  the  well-known  Grashof  criterion.  Referring 
to  Fig.  4,  the  six-bar  linkage  is  a  compound 
linkage  consisting  of  a  four-bar  linkage  and  a 
double-crank.  It  is  desirable  that  the  four-bar 
linkage  is  of  the  crank-rocker  type.  The  rota¬ 
tability  Study  involves  the  following  two  stages: 

1)  Rotatability  of  the  input  crank 

Grashof  criterion  can  be  used  directly,  we  have 

i  ♦  s  <  P+q  (11) 

where  i  and  s  represent  the  longest  and  the  short¬ 
est  link,  respectively;  the  intermediate  links 
are  denoted  as  P  and  Q. 

For  a  crank-rocker  linkage,  S  must  be  the  crank 
and  the  frame  being  either  one  of  the  adjacent 
1  inks. 


Fig.  4 

)  Rotatabil  i  ty  of  the  output  crank 

This  concerns  the  links  rs,  rg  and  rg  and  a 
iriable  length  rr  which  is  the  distance  between 
lints  0  and  E.  Referring  to  eq.(9),  we  must  have 

A2  ♦  B2  >  C2  (12) 


C  *  -(A2*B2*r2-r2)/(2r6) 
•tting  d  =  /A  ♦{}  we  have 


(d+r6)  >  r5 


r5  *  r6  >  d 


r5  *  d  *  r6 


Equations  (13),  (14),  (15),  end  (16),  therefore, 
are  the  rotatability  criteria  of  the  output 
crank.  Since  the  parameter  d  depend*  on  the 
motion  of  the  coupler  point  E,  the  coupler  curve 
of  point  E  plays  an  important  role  on  the 
synthesis  of  the  output  cranks  rg  and  rg. 

Transmiss ion- Angle  Character-'s t>'cs 

16.  The  transmission  angle  of  the  S'x-bar 
linkage  is  denoted  as  the  sraller  angle  between 
links  rg  and  rg  (Fia.  3).  The  transmission-angle 
behavior  of  the  six-bar  linkage  differs  from  that 
of  the  slider-crank  and  four-bar  mechanism,  there¬ 
fore,  it  warrants  an  investigation.  The  trans¬ 
mission  angle  x  is 

cos  x  =  ( r2+rg-d2)/(2r^rg)  (17) 


d  >  |  EG  |  *  /(xe-xg)^(ye-y(.)1 

X^  =  R2  cosf^-x) 

Yr  =  Ry  sin(j»-i) 

E  2  2  (18) 
Xg  =  ^1  +  r2 

Yg  *  R1  sinffy.,) 

(R2+R2+r2-r2-r2)-2R,R2cos(e2-3,-x+u) 
C0ST= - -Tr?-  -  - .  -  - 

r2[R,cos(yu)-a2cos(o  2"1)! 
r5  r6  ~ . . 

Equation  (19)  gives  the  transmission  angle.  It  is 
desirable  that  the  maximum  variation  of  this  angle 
from  the  right  angle  on  the  entire  rction  be 
minimized. 

Analytical  Synthesis 

17.  The  analytical  synthesis  of  a  six-bar 
linkage  with  velocity  dwell  is  introduced  in  this 
section.  The  result  provides  explicit  expression 
of  rg  and  rg  in  terms  of  other  desiyn  parameters. 
The  velocitv  dwell  equations  are  from  eqs.  (8) 
and  (9) 

A'sin  *  ♦  B'cos  *  ♦  C'  -  0 
A  sin  $  ♦  8  cos  *+0=0 


nation  (14)  obviously  holds  from  the  argument 
tflat  rg,  rc  and  d  form  a  triangle.  In  addition, 
we  Can  write 


A  -  Ye-Yg 

n  -  xt  *g 


4S 


C  -  -{A7*F?*rJ-r^)/2r6 
A'  =  Ye‘ 

B’  »  XE* 

C‘  -  -(AA,*B8,)/r6 

Solving  for  and  r^,  we  have 


r  =  V*R  2*l  AWB"Vr,"AWB 

r5  B  'A'sin.  <R'c(>Sir  A^sin^B’cnsc  ) 

-  (A  sin,  *B  rosij )] 


(?0) 


(A'A'B’Bj/fA'sinv+B'cos, ) 


(21) 


The  Optimization  Problem 

18.  Using  the  analysis  developed  before,  it  is 
possible  to  optimize  the  kinematic  design  of  the 
six-bar  linkage  for  higher-order  dwell  by  a  proper 
choice  of  the  geometric  parameters.  The  design 
problem  can  now  be  converted  to  an  optimi zation 
problem  as  follows  for  dwell  up  to  the  third  order: 


Minimize:  f  -  WjE^  +  W2E£  +  W3E3 
Subject  to: 


(22) 


(a)  dimensional  constraints  on 
desiqn  parameters 

(b)  Constrains  on  entire-motion 
characteristics,  such  as 
rotatability,  transmission- 
angle  characteristics, 
limits  on  velocity  fluctua¬ 
tions,  etc. 


In  eq.  (22),  f  represents  the  objective  function; 
W, ,  W-,  and  W,  are  weighting  factors  chosen  by 
the  designer  depending  on  design  applications;  E-j . 
E?,  and  E,  denote  as  the  error  functions 
corresponding  to  the  deviation  of  dwell  positions 
of  the  velocity  acceleration  and  shock, 
respectively,  and  are  given  as  follows 


1  3.0  cm,  r^  1  1.0  c~ ,  r^ 
a  =  2.2  cm,  b  =  2.4  c ,  »' 


3.6  cm,  r4  -  2.3  f, 
i  .  y  a'li  f  --  40'  . 


Solution:  The  dt-siun  stages  ore  qiven  in  the 
folldwinq. 

1.  Analytical  synthesis  fro erd^rp.  Equations 
(20)  and  (21)  oivc  r<,  =  1.75  <r  and  r,  =  1.75  cm. 
With  this  set  of  ri'Sign  poravters,  the  first- 
order  dwell  is  ensured. 

2.  The  optimum  synthesis  picucdure  usina  HOT 
gives  the  solution  set  as  follows:  r.  =  2.98  cm, 
r2  =  0.98  cm,  -  3.6?  cm,  r4  =  ?.3z  cm, 

rc  1  1.7  cm,  rg  =  1.8  cm,  a  =  2.2  cm,  b  =  2.4  cm, 
Q?  *  124“  and  6  -•  41,4°. 

This  set  of  design  pa r«  -Mrs  will  provide  an 
approximate  third-order  .Vel  1 .  Potion  charac¬ 
teristics  of  the  linkage  arc-  shown  in  Fig.  5. 


Fig.  5  •****• »..« 
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Automated  Dynamic  Analysis  of 
Chain-Driven  Mechanical  Systems 

A  genera t  approach  to  analyze  the  dynamics  of  chain-driven  systems  subjected  to 
transient  loads  is  developed  and  applied.  The  method  suitable  for  many  automated 
dynamic  analysis  techniques  involves  the  simulation  of  the  dynamic  effect  of  chain 
by  a  displacement  function  and  the  introduction  of  this  function  as  a  kinematic 
constraint  to  couple  with  the  system  equations  of  motion.  A  general  purpose 
dynamic  analysis  algorithm ,  the  DADS  code  ( Dynamic  Analysis  and  Design 
Systems ) .  is  then  employed  to  generate  the  set  of  system  equations  and  to  provide  a 
computer-aided  dynamic  analysis  of  the  overall  chain-driven  system.  Two  wavs  of 
formulating  the  chain  displacement  functions  are  described.  One  provides  the 
displacement  of  the  chain  based  on  the  pitch  circles  of  chain  sprockets;  the  other 
includes  a  consideration  of  the  polygonal  effect  of  the  chain  which  contributes 
essentially  to  the  dynamics  of  the  chain.  The  latter  involves  the  use  of  the  principle 
of  kinematic  equivalency,  t.e.,  modeling  the  chain  dynamic  effect  by  a  four-bar 
linkage.  Using  the  proposed  displacement  function,  the  kinematic  motion  of  the 
chain  can  be  taken  into  account.  This  procedure,  therefore,  makes  the  system 
adaptable  to  conventional  dynamic  analysis  code  in  which  the  chain  is  usually  not 
included  as  one  of  the  standard  elements.  Moreover,  pulsation  and  dynamic  load  of 
the  chain  as  well  as  the  system  dynamic  response  due  to  chum  effect  may  be 
estimated.  A  typical  large-scale  chain-driven  system  which  is  an  externally  powered 
machine  gun  is  investigated  to  illustrate  the  potential  usefulness  of  the  approach. 


I  Introduction 

Dynamic  analysis  of  any  mechanical  system  of  realistic 
complexity  is  not  a  simple  task.  These  systems  usually  involve 
a  large  number  of  mechanical  components,  such  as  gears, 
shafts,  springs,  cams,  linkages,  etc.,  and  the  interactions 
between  components  are  nonlinear.  In  addition,  loading  and 
operational  speed  may  impose  additional  complexities,  i.c., 
the  presence  of  discontinuous  forces,  mass,  velocities,  and 
potential  energies,  which  characteri/e  the  mlci  mil  lent  mol  ion 
and  inertia  loading  Irom  the  high  speed.  I  lie  problem  ol 
analyzing  the  dynamics  of  a  large  mechanical  system  is, 
therefore,  complicated. 

There  have  been  significant  contributions  recently  on  the 
study  of  dynamic  response  of  mechanisms  and  mechanical 
systems  involving  computer-aided  techniques. '  Most  of  them 
take  an  approach  which  is  general,  involving  the  treatment  of 
Lagrangian  formulation  of  the  system  equations  of  motion 
and  the  use  of  numerical  methods  to  yield  solution.  This  is 
useful  in  many  applications  concerned  particularly  with  large 
and  complicated  mechanical  systems.  System  components 
may  have  various  degrees  of  freedom  and  the  generation  of 
system  equations  in  this  case  is  not  simple,  l.arge-scale 

1  I  hr  work  ol  I  J  Haug  and  his  to- wot  kefs;  the  work  ol  M.  A  t'hacc.  V 
Orlandes.  J.  J.  locker,  etc.;  R.  E.  Beckeli.  K  C  Pan  and  S.  (  Shu.  ASM  I 
Journal  of  kngineerint  Jvr  Industry,  Aug.  1977,  pp  Mi6-67). 

Contributed  hy  the  Design  Automation  Commmee  and  presented  ai  Ihr 
Design  and  Production  Engineering  Technical  l  onference,  Washington,  l>  <\. 
Sepfemhei  12  IV  PM2.  of  Tiif  AMteicAu  Snnitv  or  Mmimshm 
I oo.inims  Miiiioscripl  reensed  al  ASMI-  lleadtiuailers.  luni-  V  IMH2  Paper 
No  *2  Ol  I  M 


kinematic  and  dynamic  analysis  codes  such  as  ADAMS, 
DADS,  DRAM.  DYMAC,  IMP,  KINSYN.  UNCAGES,  and 
others  have  significantly  facilitated  the  implementation  of 
modern  analytical  developments  in  mechanisms  to  the 
solution  of  realistic  problems  which  may  be  otherwise  dif¬ 
ficult  to  obtain. 

Of  particular  interest  and  concern  in  this  investigation  is  the 
study  of  llic  dynamics  ol  clt.iin  dmon  nicch.mic.il  systems 
I  hoc  arc  several  reasons  motivating  t his  study.  I  irslly,  there 
is  a  real  need  for  the  dynamic  analysis  of  systems  of  this  kind. 
In  recent  times,  primarily  because  of  low  cost,  high  efficiency, 
and  capability  to  operate  at  high  loads,  chain  drives  have 
become  one  of  the  most  popular  and  economic  means  of 
power  transmission.  Second  and  most  importantly,  methods 
and  approaches  are  not  available  to  treat  chain-driven  systems 
and  there  is  not  much  work  reported  in  the  literature  regard¬ 
ing  chain  drives.  Moreover,  modern  kinematic  and  dynamic 
codes,  to  the  author's  knowledge,  have  not  yet  included  chain 
as  a  basic  clement  in  the  codes,  and  consequently,  these  codes 
are  not  directly  applicable  to  chain-driven  systems.  The 
problem,  therefore,  warrants  our  investigation. 

We  begin  with  a  dynamic  analysis  of  the  chain  drive.  A 
dynamic  model  of  the  chain  is  presented  and  it  is  basic  to  the 
determination  of  the  dynamic  response,  such  as  dynamic  load 
and  motion  characteristics  of  the  chain.  A  typical  large-scale 
mechanical  system,  an  externally  powered  modern  machine 
gun  involving  chain  drive,  is  then  investigated.  It  is  used  as  an 
illustration  of  the  proposed  approach.  Ihc  kinematic 
modeling  and  dynamic  modeling  using  an  effective  code. 
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DADS,  ate  given.  I  mulls,  iHimuic.il  results  ate  pi  esciilcd  and 
conclusions  are  drawn. 

2  The  Dynamics  of  Chain  Drives 

In  the  following,  we  give,  first,  a  discussion  on  the  chain 
polygonal  action  which  affects  the  velocity  fluctuation  of 
chain  based  on  a  dynamic  model.  Secondly,  we  give  the 
derivation  of  the  position,  velocity  and  acceleration  equations 
of  chain  drive.  A  chain  displacement  function  which 
represents  essentially  the  motion  characteristics  ol  the  chain  is 
presented  and  is  used  in  our  study  of  the  dynamics  of  the 
chain-driven  system. 

2.1  Polygonal  Action.  The  action  of  a  chain  as  it  runs  over 
a  sprocket  can  be  compared  to  a  nonslipping  thread  running 
over  a  regular  polygon  prism.  When  the  chain  engages  a 
sprocket,  each  link  falls  on  a  chord  of  the  sprocket  pitch 
polygon  rather  than  the  true  pitch  circle.  Therefore,  the  active 
pitch  diameter  of  the  sprocket  varies  between  the  limits  of  a 
circle  inscribed  in  the  pitch  polygon  and  a  circle  cir¬ 
cumscribing  the  pitch  polygon.  The  chain  lying  on  the 
sprocket  forms  a  polygon  rather  than  a  circle  resulting  in  the 
“polygonal  action"  of  the  chain  drive. 

The  speed  of  the  chain,  therefore,  fluctuates  even  with  the 
constant  driving  speed  of  the  sprocket  due  to  this  "polygonal 
action."  Impact  occurs  because  of  the  difference  in  the 
velocities  of  the  point  on  the  roller  and  the  point  which  comes 
into  contact  with  the  sprocket.  This  results  in  cyclic  ac¬ 
celeration  and  deceleration  which  creates  dynamic  loads  in  the 
chain  and  often  causes  excessive  noise  and  vibration.  The 
effects  of  impact  become  complex  and  significant  at  higher 
speeds.  Especially  if  the  frequency  of  the  velocity  fluctuation 
coincides  with  a  resonant  frequency  of  the  chain  system,  large 
stresses  may  occur  which  can  cause  fatigue  failure. 

Some  aspects  of  polygonal  action  have  been  discussed  in  the 
literature  1 1  -5] .  Several  interesting  approximations  have  been 
made  of  the  driven  shaft  angular  velocity  charactciistics 
because  ol  the  polygonal  died  ol  sptockcls.  A  closet 
examination  of  the  chain  drive  components  shows  that  the 
motion  of  the  system  is  actually  transmitted  through  a  series 
of  four-bar  linkages  (Tig.  I).  Using  this  concept,  Morrison  |3| 
analyzed  the  variation  of  the  angular  speed  of  (he  driven 
sprocket  and  found  that  by  increasing  the  number  of  teeth  in 
the  sprockets  the  difference  between  pilch  polygon  and  pitch 
circle  will  reduce  and  so  does  the  polygonal  effect.  Bremer  |?| 
recommended  that  the  sprocket  should  at  least  have  30  teeth. 
In  the  particular  case  that  the  drive  sprockets  have  a  one-to- 
one  ratio  of  their  pitch  circles  and  their  center  distance  equals 
a  whole  number  of  link  pitches,  the  chain  will  transmit 
uniform  angular  motion  through  the  driven  sprocket.  Then, 
there  is  no  speed  fluctuation  in  the  driven  sprocket.  The  chain 
links,  however,  are  still  subject  to  pulsations.  There  is  a 
polygonal  effect  on  the  connecting  chain  links,  which  is  due  to 
the  speed  fluctuation  of  each  link  along  the  linear  direction  of 
its  chain  path.  In  this  case,  the  motion  of  the  chain  links  and 
their  masses  become  significant  in  the  design  consideration  of 
the  chain  drives. 

In  chain  drives,  there  are  two  ways  of  transmitting  the 
power.  In  the  common  case,  load  is  transmitted  through  the 
driven  sprocket.  In  other  cases,  the  load  may  he  carried  by  a 
stud  on  the  chain  which  makes  (he  driven  sprockets  serve 
essentially  as  idlers.  In  the  latter  case,  usually  the  impact 
effect  between  chain  rollers  and  sprocket  teeth  can  be 
neglected.  1  he  masses  of  chain  links  and  the  masses  attached 
to  the  chain  links,  however,  would  still  be  subject  to  velocity 
fluctuation  and,  therefore,  generate  dynamic  loads.  The 
motions  characteristics  of  the  chain  links  and  dynamic  loads 
are  analyzed  and  presented  in  the  next  section. 

2.2  I'uKatinn  and  Dynamic  Iziadi.  As  mentioned  in  the 
previous  section,  the  chain  drive  can  be  ircaicd  by  an 


fOUR-BAH  LM KA',t 


Fig.  1  Equivalent  tour-bar  linkage  ol  the  chain  drive 


equivalent  mechanism  of  a  four-bar  linkage.  The  position, 
velocity,  and  acceleration  of  any  point  in  the  chain  can. 
therefore,  he  expressed  mathematically  in  terms  of  the  four 
bar  parameters  and  the  driving  crank  angle.  A  convenient 
method  to  derive  these  motion  expressions  is  to  use  complex 
numbers  to  express  the  vector  loop-closure  equation  of  the 
equivalent  linkage  and  to  obtain  its  time  derivatives.  I  lie 
following  rept events  such  an  apptoach. 

Take  a  point.  T,  for  instance,  in  the  chain  link  as  shown  in 
Tig.  2.  We  are  interested  in  deriving  its  motion  expressions,  in 
particular,  its  linear  velocity  and  acceleration  in  the  direction 
of  the  chain  path.  The  point,  I’,  may  be  considered  as  a  point 
on  the  coupler  of  the  equivalent  four-bar  linkage.  Referring  to 
Tig .  2  and  denoting  cz,  h.  <\  and  d  as  the  link  length  of  the  four 
bars,  wc  have  the  position  vector  of  P  in  complex  number 
rotation. 

h  —  tif  "  i  (1) 

where  the  parameters  0, ,  0,,  and  0,  denote  the  angular 
positions  of  the  links  as  shown  in  Tig.  2;  the  vectors  a  =  0.1 
and  i(  =  .-f/>. 

Differentiating  equation  (I)  with  respect  to  time,  we  obtain 
the  velocity  of  the  point.  P.  as 

f  =  ftfu.v"'1  +  (2 ) 

The  real  and  imaginary  components  are  then, 

R ,  =  —  trw,  sin  0,  -gu;;sin0: 

/,  =  flu.', cos#,  4-gu!;COS0;  (3) 

Turther  differentiating  equation  (3)  with  respect  to  time,  we 
obtain  the  acceleration  of  the  point,  P,  Denoting  the  angular 
accelerations  as  </u.v  f/z  =  «i  andrfu =  <«- ,  we  have 

a  —  (  —  u.v  +  «>i  hie'"1  ♦  (  —  «-■  *  ««; (4) 
The  real  and  imaginary  parts  are. 

Rd  =  -flfwv’cos  0,  +  a |  sin  0,) 

-£( u>  - ' cos  +  or- sin  0; ) 

/.,  -  -o(u.ysin  0,  r»|COsfl,) 

-  n  (<*);-’  sin  -  rt-COS  0;)  (5) 
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ketcrrmg  to  I  ig.  3,  it  is  appaient  that  the  real  parts  of 
equations  (3)  and  (3)  represent  the  linear  velocity  and  ac¬ 
celeration  of  the  point,  P,  parallel  to  the  A1 -direction  of  the 
chain  path,  and  the  imaginary  parts  are  the  velocity  and 
acceleration  perpendicular  to  the  x -direct  ion. 

In  the  case  that  the  driving  and  the  driven  sprockets  having 
a  one-to-one  ratio,  the  linear  velocity  (V , )  and  acceleration 
(a, )  equations  along  the  direction  of  the  chain  path  of  the 
contact  point  between  the  driving  sprocket  and  the  chain  links 
become  very  simple.  They  are 

y,  =  -tiu),  sin  9, 

a i  =  —  a(wi:cos  8,  +  a, sin  9, )  (6) 

Furthermore,  if  the  driving  sprocket  is  running  at  constant 
speed,  we  have  <»,  =  Oand 

V ,  =  -Uo>|Sin  0| 

a,  =  -Uuvcos  9,  (7) 

The  chain  pulsates  with  a  period  which  is  equal  to 

I  =  6/0/ Nn  s  (8) 

where  N  denotes  the  speed  of  the  driving  sprocket,  and  n 
denotes  the  number  of  teeth  on  the  sprocket.  Thus,  the 
frequency  of  this  pulsation  is 

/  =  A/n/60  pulses  per  s  (9) 

If  this  frequency  is  close  to  the  natural  frequency  of  the 
chain  system,  the  amplitudes  in  the  T-direction  can  reach  a 
dangerous  magnitude. 

Now,  consider  that  one  chain  link  engages  with  the  driving 
sprocket  in  the  range  of  angle  90  deg  -d  to  90  deg  +  d(Fig.  3). 
By  using  equation  (7),  the  accelerations  at  the  instants  that  a 
particular  link  engages  the  sprocket  while  an  adjacent  one 
leaves  the  sprocket  can  be  written  as 


^rn|il|K  —  tlU)  Sin  d 


0|c,v,  =uu),:sin  d  (II) 

where,  w,  denotes  the  angular  velocity  of  the  sprocket  and  it  is 

=  30  radA 


d  =  sin  '  (  ~~  )  ,12) 

The  instantaneous  increase  in  chain  acceleration  as  a  new 
sprocket  tooth  engages  the  chain  is,  therefore,  equal  to 

d<r  =  ,<'cn»  -tficvt*  =2tftti,’sind  (13) 

After  substituting  w,  and  d  into  equation  (13),  we  have 

/  wN  \  -  w 
Au  =  2ol  -  I  sin  — 
v  30  /  n 

t  \  ; 

=  (,4) 

It  is  this  change  of  acceleration  which  results  in  an  in¬ 

stantaneous  dynamic  load.  We  have 

-)**"£]  (,5) 


where  L,  IT,  denotes  the  sum  of  the  weight  of  all  chain  links 
plus  the  equivalent  weight  of  all  the  following  cascaded 
mechanical  components. 

Equations  (l)and  (15)  which  represent  the  general  form  of 
the  chain  displacement  equation  and  the  instantaneous 
dynamic  load,  respectively,  are  useful  in  the  dynamic  study  of 


chain-driven  systems.  Depending  upon  the  particular  system 
in  question,  expressions  similar  to  these  can  be  derived  and 
this  will  be  demonstrated  in  a  later  section. 

3  Dynamic  Modeling  of  an  Automatic  Weapon 
System— The  Model  X  Gun 

In  this  section,  we  would  like  to  illustrate  the  use  of  chain 
displacement  function  in  connection  with  an  effective 
dynamic  analysis  code  to  study  the  dynamics  of  a  complicated 
large-scale  mechanical  system  with  chain  drive.  It  is  the  study 
of  a  new  externally  powered  machine  gun  hereinafter  referred 
to  as  the  Model  X  gun. 

Before  a  description  of  the  kinematics  of  the  gun,  it  is 
essential  to  give  a  brief  introduction  of  the  dynamic  computer 
code,  DADS,  which  is  chosen  in  this  investigation. 

The  DADS  code  (Dynamic  Analysis  and  Design  System) 
developed  by  Haug,  Wchagc,  and  their  co-workers  |b,  7]  is  a 
powerful  computational  package  for  the  analysis  and  design 
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ol  mechanical  systems.  The  code  di tiers  from  other  computer 
packages  not  only  in  its  computational  method,  but  also  in  its 
approach  to  treat  mechanisms  or  mechanical  systems  with 
intermittent  motion.  The  code  utilizes  logical  spring  and 
dampers  |8,  9|  to  simulate  (he  jump  conditions  which  often 
occur  for  such  systems.  This  makes  the  code  particularly 
attractive  to  the  dynamic  study  of  mechanical  systems  with 
intermittent  motion,  such  as  the  automatic  weapon  system.  In 
the  following  sections,  the  DADS  code  is  applied  to  studs  the 
dynamics  of  the  Model  X  gun. 

The  kinematics  of  the  Model  X  gun  is  determined  essen¬ 
tially  by  the  motion  characteristics  of  the  bolt -carrier 
assembly  (i.e.,  bolt-carrier  group)  shown  in  Fig.  4,  which  is  a 
mechanical  system  with  intermittent  motion. 

The  motion  of  the  bolt-carrier  contains  a  sequence  of  events 
involving  finite  masses  capture  and  release  as  well  as 
mechanical  stops.  The  bolt-carrier,  driven  by  the  stud  of  a 
chain  which  moves  along  a  definite  near-square  path  (l  ig.  5), 
moves  in  a  reciprocating  motion.  Various  gun  operational 
stages  are  performed  as  the  chain  undergoes  a  cycle. 

Let  us  consider  now  the  modeling  of  the  driving 
mechanism,  in  particular,  the  bolt-carrier  and  its  associated 
components,  in  order  to  be  able  to  adapt  to  the  DADS  code. 
We  give  first  a  terminology  of  components.  The  basic  ap¬ 
proach  of  modeling  is  presented  next.  Then,  the  kinematic 
motion  cycle  of  the  bolt-carrier  assembly  is  outlined  and 
graphically  displayed. 

Terminology.  The  bolt-carrier  assembly  consists  of  the 
following  components  (Fig.  4): 

Mass  No.  Part  Name  and  Function 

1  gun  body;  taken  as  reference  frame 

2  bolt-carrier;  the  part  that  carries  the  bolt,  striker, 
buffer  spring  and  retainer  and  plug  in  the  motion 

.1  striker;  striking  the  fire  pin  to  start  propellant 
ignition 

4  striker  plug  and  retainer;  retaining  the  buffer 
spring  in  compression 

5  buffer  spring;  the  lumped  mass  of  butter  spring 
that  functions  as  buffer  and  supplies  the  energy 
for  striking 

f>  sear;  the  portion  of  the  mechanism  which  controls 
the  firing  of  the  gun 

7  bolt;  the  body  that  locks  the  cartridge  into  and 
extracts  it  from  the  chamber  through  camming 
action  relative  to  the  bolt-carrier 

In  the  dynamic  model,  the  friction  loice  between  all 
lianslalioiial  |omls  ate  assumed  to  he  zero.  Mass  5 
corresponds  to  the  lumped  mass  of  the  buffer  spring.  It  is 
assumed  to  be  located,  initially,  at  the  center  of  bodies  )  and  4 
and  is  precompressed  to  the  length  of  6.045  cm  or  2. .18  in. 

Basie  Approach.  I  lie  driving  mechanism  of  the  Model  X 
gun  is  a  typical  example  of  a  large-scale,  complex  in¬ 
termittent-motion  mechanism.  It  is  difficult  to  model.  The 
application  of  DADS  to  analyze  its  dynamic  response  is  not  a 
simple  task.  In  the  following,  the  basic  approach  of  modeling 
is  given. 

(1)  The  model  includes  only  the  bolt-carrier  assembly. 
Other  components  of  the  power  train  are  not  taken  into 
consideration. 

(2)  The  kinematic  motion  of  the  boll -carrier  sci  ves  as  input 
of  the  code,  which  is  supplied  as  an  externally  user-supplied 
constraint  function. 

The  displacement  function  of  the  bolt-carrier  obtained 
from  the  kinematic  analysis  computer  code  of  the  Model  X 
gun  |I0|  is  applied  to  DADS  as  a  constraint.  The  code  in¬ 
corporates  this  constraint  with  the  system  equations  of 
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Fig  S  Operational  cycle  ol  the  Model  X  gun 


motion  which  are  then  solved  simultaneously  and  iteratively 
at  each  time  increment  to  match  the  value  of  stare  variables 
corresponding  wiih  the  constrained  displacements. 

(3)  Logical  springs  and  dampers  ate  used  to  simulate  the 
mass  discontinuity  and  the  impact  phenomena.  The  selection 
of  the  magnitudes  of  these  parameters  and  the  proper 
locations  to  place  them  are  olten  a  heuristic  approach  defers 
to  Fig.  4).  Details  will  be  discussed  later. 

The  Motion  Cycle.  I  he  motion  cycle  ol  the  bolt-carrier 
assembly  involves  two  rectilinear  motion  periods  (forward 
and  rearward  motions')  and  two  dwell  periods.  The  dwell 
period,  referring  to  Fig.  5.  corresponds  esentially  to  the  in¬ 
dexing  period  of  the  Geneva  mechanism  (stage  I),  and  the 
firing  and  the  gas  escaping  period  (stage  >).  In  i his  study ,  only 
the  forward  and  rearward  motions  are  of  concern.  A  detailed 
description  of  these  motions  is  given  in  the  following: 

(I)  Forward  Motion.  The  forward  motion  involves  tiie  gun 
operational  Stages  2.  3.  and  4  of  l  ie  5  In  terms  of  the 
progressing  of  time,  we  have  the  following  events: 

(<z)  At  i  t),  the  boh  carrier  is  ready  to  move  rightward, 
and  bodies  3.  4.  and  '  are  in  the  extieme  letr  positions  I  he 
velocity  and  displacement  ol  each  body  are  at  its  initial  value. 
The  sear  is  stationary  against  the  gun  body  (body  I )  (Fig.  6). 

(/>)  At  l  =  i,  >0.  the  bolt-carrier  starts  to  move  rightward 
with  the  bolt  (body  ')  and  the  striker  (body  T );  however,  the 
striker  plug  and  retainer  (body  4)  remains  stationary  (l  ie.  '). 

(<•)  At  l: ,  the  bolt-carrier  captures  the  sinker  plug  and 
retainer  (body  4)  and  both  move  toward  the  right.  At  this 
moment,  all  the  components,  except  seat,  .lie  moving  at  the 
same  speed  and  diicclion  (I  ig.  S) 

(r/)  At  r,,  the  striker  (body  3)  makes  contact  with  the  sear 
(body  6).  and  the  buffer  spring  is  under  compression  due  to 
this  effect.  The  rest  of  the  masses  (bodies  2.  7,  and  4)  con¬ 
tinue  to  move  rightward  (Fig.  9). 

(e)  At  r4,  the  bolt  (body  ')  is  locked  on  the  gun  barrel,  and 
the  boll-carrier  (body  2)  is  in  contact  with  sear  (body  6)  I  lie 
spring  between  the  sear  and  the  gun  bodv  (body  It  is  in 
compression  (Fig.  10). 

(f)  At  the  striker  is  reads  for  firing.  All  other  com¬ 
ponents  (bodies  2,  4.  6.  and  7)  are  stopped  and  are  going  into 
the  dwell  stage  (Fig.  1 1 ). 

The  lime  events  described  above  are  graphically  displayed 
in  Figs.  6  through  II.  The  logical  spring-damper  pairs 
(referring  to  Tig.  4)  arc  used  lo  simulate  the  finite  masses 
capture  and  release,  which  arc  defined  between  body  2  and 
bodies  3.  4,  6,  and  7.  In  addition,  they  are  used  to  simulate  the 
mechanical  stop  (i.e..  impact),  which  are  defined  between 
body  I  and  bodies  2,  3,  4.  6.  and  7. 

1  Oiw.mi  ntonoii  )x  .»v  tlir  molmn  i  \»  l«*  K  loir  Tim  up 
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Fig.  8  Initial  position  of  tha  bolt-carrier  in  forward  motion 


Fig.  10  Tha  position  that  tha  bolt  carrier  contacts  the  sear 


Fig.  7  Bolt-carrier  assembly  starts  to  move  In  the  rightward  position  Fig  11  The  position  ot  gun  tiring 


3M/ Vol  105,  SEPTEMBER  1983 


T ransactlons  ot  the  ASME 


r ' 

i 


(2)  Kearward  Minion.  I  lie  real  waul  motion  corresponds 
to  the  gun  operational  stages  6,  7,  and  8  (Fig.  5).  This  part  of 
the  motion  cycle  is  simulated  with  the  consideration  that  (he 
bolt  (bods  7)  and  the  striker  (body  3)  are  keyed  together  as  a 
single  body,  denoted  as  body  3.  This  is  because  there  is  no 
relative  motion  between  them  during  the  period  of  motion. 
The  time  events  arc  given  as  follows: 

(u)  At  /  -II,  the  bolt-carrier  (body  2)  and  its  associated 
components  are  ready  to  make  leftward  motion  while  the 
buffer  spring  is  compressed.  The  velocities  and  displacements 
of  all  members  are  at  their  initial  conditions  (Tig.  12). 

(b)  At  /  =  /,  >0.  the  bolt-carrier  starts  to  move  leftward 
with  the  striker-plug  and  retainer  (body  4).  In  the  meantime, 
the  bolt  and  striker  (body  3)  remain  stationary.  I  his  is 
because  the  bolt  is  locked  on  the  gun  barrel  and  the  striker 
sustains  a  spring  force.  The  buffer  spring  is  gradually  released 
to  its  precompressed  length.  The  sear  moves  at  the  same  speed 
as  the  bolt-carrier  (Tig.  13). 

(c)  At  /;,  the  bolt-carrier  captures  mass  3  (striker  and  boll), 
and  together  they  move  toward  the  left.  The  sear  is 
mechanically  stopped  by  the  gun  body  (body  I )  (Tig.  8). 

(a/)  At  /,,  the  striker-plug  and  retainer  (body  4)  contact  the 
gun  body  (body  1)  and  they  are  stopped  mechanically  (i.e.. 
impact).  The  buffer  spring  is  then  compressed.  The  bolt- 
carrier,  the  striker,  and  the  boll  continue  to  move  leftward. 
The  sear  remains  motionless  against  the  body  1  (Tig.  7). 

(e)  At  i, ,  the  movement  of  the  bolt -carrier  is  ceased  and 
goes  into  the  dwell  stage.  In  the  meantime,  all  other  members 
are  also  in  dwell  (Fig.  6). 

l  ogical  spring-damper  pairs  are  defined  between  body  2 
and  bodies  3,  4,  and  6  to  simulate  the  mass  caput. e  and 
release.  In  addition,  they  are  defined  between  body  1  and 
bodies  2,  4,  and  6  to  simulate  the  mechanical  stop. 

The  dynamic  response  during  the  dwell  stages  is  not  of 
concern  in  the  modeling.  This  is  because  at  the  dwell  stages  of 
the  motion  cycle,  the  driving  chain  stud  moved  along  the 
vertical  direction  ol  the  chain  path  and  there  is  no  telative 
motion  in  the  A -direction. 

A  summary  of  the  time  events  and  the  corresponding 
graphical  representation  of  the  system  are  given  in  I  able  1 . 


4  Chain  Displacement  Function 

The  expression  of  the  chain  displacement  function  depends 
upon  the  particular  chain  derived  system  in  question.  In  the 
following,  chain  displacement  functions  aie  formulated  for 
l lie  Model  X  gun  I  he  approach,  however,  ts  general  and  can 
be  extended  to  any  chain-driven  system. 

Denoting  the  motion  parameter  of  the  bolt -carrier  along  the 
longitudinal  coordinate  as  A',  and  the  displacement  of  the 
chain  stud  along  the  chain  path  as  8  reference  positions  of 
both  -V  and  8  arc  shown  in  Tig.  5.  The  chain  path  can  be 
specified  by  eight  positions  or  stations  which  divide  the  path 
into  four  straight  and  four  curved  portions.  I  here  are  two 
ways  to  formulate  the  chain  displacement  function:  I)  the 
kinematic  displacement  function;  and  2)  the  dynamic 
displacement  fund  ion . 

(1)  The  Kinematic  Displacement  Function.  In  (his  ap¬ 
proach  the  chain  sprocket  is  physically  replaced  by  the 
sprocket  pitch  circle.  Polygonal  effect,  or  the  dynamic  effect 
of  the  chain  is,  therefore,  neglected.  The  displacement  func¬ 
tion  formulated  with  this  assumption  is  called  the  kinematic 
displacement  function  and  it  is  described  for  each  operational 
stage  as  follows; 

I .  Stage  I,  0  <  S  <  3  81  cm 
S  r  (ll  «,).«,  (Meg 
X  -  6.35  t  2/-, 


Table  1  Time  events  representing  the  motion  cycle  ot  the  gun 
mechanism 


r  i  t ure 

PorvarT 
Event  s*c. 

’tirua;-1. 

•  Event  No . • 

'  seu&ert » 
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Sxtrtun* 

?c«-t:cr. 

3 
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c 

r.-.t*rs»  -  .a'-* 

+ 

d 

1  or. 

txz rijh: 

petition 


•R*  fe  r  to  Section  3. 


X 


Fig.  14  Displacement  curve  of  (he  boh  carrier  in  kinematic  motion 


2.  Stage  2  (buffering),  3.81  <  8  <  6. 35  cm 

8=  3.81  *  r, (()-(!.),  <>.  135  deg 

.V  =  r\  +6. 35  +  /\cos</'  -ft- ) 

3.  Stage  3.6.35  <  ,S  *  12.7  cm 
8=12.7  t-MO-fh).  tlx  -  225  deg 
•V  =  r,  +6.35  -  r.lfi  -  tl ,) 

4.  Stage 4,  12.7  <  S  <  15.24cm 
8=  12. 7  >-rje-H,).U4  =450 deg 
A  =  r,  -r,sin(fi  -  to) 

5.  Stage  5,  15.24  <  8  v  If  (is  0„ 

8=15.24  t  r\(fl-(MA  540 deg 
X=0 

6.  Stageb,  19. fv  <  S  <  21.59cm 
8=  19.05  +  r,  <0  — (UA  =675  deg 
X  =  r,  -  r,cos(tf  -  «,. ) 

7.  Stage  7.  21 .59  <  5  <  27.94  cm 
8=21.59  +  rXH-0-  )A  =  ’65  deg 
X=r.  +/\  (/)-»-) 

8.  Stage  8,  27.94  s  8  30.48  cm 

8=27.94  +  r.  {II -IK). II.  =99(ldcg 
X  =  8.0  4  r, siti( (I  -  6. ) 

By  specifying  the  crank  angle,  ».  the  location  ot  the  vhain 
stud  and  the  displacement  of  the  bolt-carrier  along  the  A  avis 
can  be  calculated.  The  calculated  displacement  cutve,  shown 
in  Tig.  14,  matches  very  closely  to  the  experimental  result  ( I  Of 
which  is  shown  in  Tig.  15,  using  an  Optron  d\ namic  tracker  1 

(2)  The  Dynamic  Displacement  Function.  The  dynamic 
displacement  function  implies  that  the  dvnanm  of  the  ih.un. 

'  Mr*  fto  <»|*ll»  ,tl  xlc . pl.it  cnti-Hf  Mind*  |  ‘-f,  j  *  »|>tt.«it  l  ••!)*  \\  •.  »d 
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TIME  (MSI 

Fig.  IS  The  experimental  displacement  curve  ot  the  bolt-carrier 

namely,  the  polygonal  effect  is  incorporated  in  the 
displacement  function  using  a  model  of  four-bar  linkage. 
Considering  that  the  motion  of  chain  stud  in  the  forward  and 
rearward  direction  can  be  disided  into  three  operational 
stages  as  shown  in  f  ig.  5,  the  displacement  conations  can  be 
defined  as  follows; 

forward  Motion  and  Kearward  Motion4 

.1 .  Stage  2,  (><l)<  Hit  deg  ( Stage  6  lor  Rearward  Mot  ion  ) . 
In  litis  stage,  the  chain  stud  moves  along  the  pitch  circle  of  the 
sprocket  \S  e  have 

A  —  1  cos  H)  (16) 

where  r,  denotes  the  pitch  circle  of  spiocket  and  it  denotes  the 
crank  angle. 

It  Stage  3,  hi i  deg  <tl<M5  deg  [Stage  '  tot  Rearward 
Motion).  The  chain  stud  in  this  stage  is  in  the  span  of  chain 
links.  As  mentioned  ill  Section  2,  the  chain  dine  can  he 
modeled  as  an  equivalent  lout  hat  linkage,  lhe  chain  stud, 
therefore,  can  be  considered  as  a  point  that  moves  along  the 


coupler  link 

and  is  subject  to  periodic  motion.  This 

can  be 

seen  by  referring  to  l  ig.  3.  We  have,  therefote 

(i)  If 

tf,  >22.5  deg 

v  =  r,  +  (,V+  1)/. -r.sintT*  deg 

(P) 

(ii)  If 

it,  <22.5  deg. 

v  =  r.  +  .V  L  +  r,sin  0, 

(IS) 

where  Ht 

It  90  deg  45  deg  •  A',  ,V  is  the 

integer 

patt  of 

(II  ‘HI  deg  ■ 

15  deg)  and  /  is  the  link  length. 

1  oi  eve 

iv  22. 6 

degiccs  ot  It 

, ,  the  chain  stud  advances  lot  a 

distance 

of  one 

chain  link. 

C.  Stage  4,  M 5  deg  <t><  405  deg  (Stage  $  for  Rearward 
Motion).  In  this  stage,  the  path  ot  the  chain  stud  coincides 
with  (he  pitch  circle  again,  and  is  moving  along  the  pilch 
circle.  We  have,  therefore 

x  =  r.  +  6.35  +  r, sin  (0-315  deg)  (Id) 

The  displacement  function  of  the  chain  is  used  as  a 
kinematic  constraint  which  serves  as  the  input  to  the  DAI3S 
code  and  it  is  coupled  with  the  system  equations  ot  motion  to 
represent  the  dynamic  state  of  the  chain-driven  system. 

5  The  Dynamic  Response  of  the  Model  X  (inn 

Hie  output  of  (he  DADS  code  provides  information 
concerning  essentially  two  phases  of  study:  the  initial  phase 
and  the  dynamic  analysis  phase. 

1  III  lilt-  ll-.IIW.lll!  Illiltintl,  llic  Ulsgl.K  0111*11!  1-1(11.1111111'  -III-  llll-Hllt.il  .1.  itlflM-  111 
fht‘  Ifl'itfofl,  CAvE'pf  the  (lllflfMMI  ol  rllOlKDI 
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Fig  16  The  displacement  curves  of  the  bolt-carrier 


5.1  The  Initial  Phase.  The  initial  stage  refers  to  the 
assembly  phase  of  the  system  equations  of  motion  in  the 
DADS  program.  This  involves  the  formulation  o!  Jacobian 
mattiN  of  partial  derivatives  for  both  the  forward  and 
rearward  motions.  The  following  information  must  he 
supplied:  the  position  (v.r.ol  and  velocity  (v.i.o)  ot  bodies  2 
through  7;  the  constraint  forces  which  arc  generated  to 
maintain  these  initial  conditions  for  bodies  at  specified 
positions;  information  on  the  spring-damper  pairs  (such  as 
velocity,  spring  length,  torce  in  v  direction  and  lorcc  in  )• 
direction)  and  the  spring  and  damping  cocllicicnts  for  the 
logical  v|»  ingdanipei  pan  -  w  tilt  t  espe,  t  to  time:  ami  the  usei  - 
supplied  constiaim  value. 

5.2  lhe  Dynamic  Analysis  Phase.  I  he  dynamic  analysis 
phase  is  initiated  when  the  static  initial  condition  is  removed. 
Physically,  this  means  the  chain  stud  begins  to  drive  (lie  boll 
carrier  to  move  along  the  chain  path.  Certain  conclusions  can 
be  made  as  follows. 

(1)  Hie  displacement  cuivcvol  the  bolt  c.ittiet  (body  2)uie 
plotted  in  l  ig.  Ih  lot  totw.ud  and  tearward- motion.  I  hesc 
curves  represent  the  constrained  motion  ot  the  bolt-cat  tier 
along  die  A-diiection.  I  lies  match  with  the  displacement 
curve  reported  in  (l(l|.  litis  is  expected  because  the 
displacement  function  serves  as  i fie  external  constraint 
equation  of  the  mot  on  of  the  bolt-carrier. 

(2)  The  velocity  cutvcs  of  the  boll-carrici  are  shown  in 
lig..  17  and  IS.  I hese  curves  represent  the  solution  of  system 
equations  of  motion  at  each  time  increment  using  different 
chain  displacement  functions. 

(3)  Ttie  velocity  ol  other  components  can  be  readily  ob¬ 
tained  liom  the  computet  output,  l  ot  instance,  the  fotw.tid 
motion  of  the  sttiker  is  shown  in  I  ig.  If.  The  curve  represents 
the  characteristics  of  intermittent-motion.  The  velocity  of 
body  3  initially  is  the  same  as  body  2  because  of  the  mass 
capture  effect.  After  a  scry  short  time  interval,  until  t  0.025 
s,  body  3  is  stopped.  This  is  done  through  the  logical  spring- 
damper  pair  to  simulate  the  impact.  The  sudden  change  of 
velocity  gives  rise  to  a  jump  condition  and  demonstrates  the 
behav  ior  of  intermittent-motion, 

The  effect  of  chain  dynamics  on  the  contact  forces  between 
components  is  significant.  This  can  be  seen  by  a  comp.uisou 
of  the  logical  springer-damper  forces  between  the  hoi. -earner 
(body  2)  and  the  striker  plug  (body  4),  the  scar  (body  6)  and 
the  bold  (body  7)  using  different  displacement  (unctions. 
Numerical  results  and  some  typical  comparisons  arc  shown 
graphically  in  bigs.  20  and  21.  f  or  instance,  the  contact  lorcc 
between  bodies  2  and  4  (refer  to  1  ig.  5,  during  the  chain 
operational  stage  3  and  a  portion  of  stage  4)  deviates  con¬ 
siderably  for  both  approaches.  This  is  because  the  fluctuation 
of  velocity  is  totally  neglected  in  the  analysis  using  lhe 
kinematic  displacement  I  unction. 

Transactions  of  tho  ASME 
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Fig.  1 7  The  velocity  ot  the  bolt-carrier  in  (onward  motion 


TIME  ,  sec 
ooi  002 


TIME  ,  sec 

Fig.  t8  The  velocity  of  the  bolt-carrier  in  rearward  motion 


X  ,  cm 

Fig.  20  Contact  force  between  the  bolt-carrier  (body  2)  and  the  striker 
plug  end-retainer  (body  4)  in  forward  motion 


(4)  The  present  chain  drive  of  the  Model  X  gun  is  a  good 
design.  It  apparently  has  already  utilized  some  of  the  design 
criteria  discussed  and  found  in  this  investigation.  The  driving 
and  driven  sprockets  arc  of  the  same  size.  In  addition,  the 
chain  link  length  to  pitch  ratio  is  even.  I  herefoic.  a  major 
part  of  the  polygonal  effect  is  essentially  eliminated. 

(5)  t  he  dynamic  load  of  the  chain  is  estimated  in  the  next 
section.  It  can  he  reduced  by  increasing  the  number  of  teeth  in 
the  sprocket.  Such  design  should  minimize  noise  and 
vibration,  and  prolong  the  life  of  the  chain.  Moreover,  the 
vibration  transmitted  through  the  bolt-carrier  assembly  can 
also  be  reduced. 

5.3  Dynamic  load.  One  of  the  important  pieces  of  in¬ 
formation  is  the  dynamic  load  generated  by  the  chain  links 
and  this  can  be  easily  estimated  based  upon  the  theory 
developed  before. 

In  the  Model  X  gun,  the  chain  has  24  links  and  it  is  30.58 
cm  or  12  in.  in  length  and  has  a  mass  of  214.32  gm  (or  weighs 
0.4725  lb).  The  dynamic  load  at  an  operational  speed  of  the 
gun,  say  600  shots  per  minute,  may  be  predicted  with  the  help 
of  equation  (15).  We  have  the  following: 

I  cl  /V  =  1875  rpm 

I  ~  1 .27  cm  (0.5  in.) 
n  =  8 

Id  !•',  /K )  720.4  gm 

I  hen,  ci|uuli<»t  (IM  gnex  a  maximum  dynamic  load. 


8j  - WITH  CHAiN  DYNAMIC  EFFECT 


- WITHOUT  CHAIN  DYNAMIC  EFFECT 


f 


Fig.  21  Contact  force  between  the  bolt-carrier  (body  2)  end  the  striker 
(body  3)  In  reerwerd  motion 


h.,  =  352. 74  N. 

The  magnitude  of  this  dynamic  load  is  large.  Therefore,  its 
effect  is  significant.  It  is  this  load  which  generates  part  of  the 
noise,  v  ibration  and  stresses. 

6  ( 'oneltisions 

A  general  approach  utilizing  the  DADS  code  to  analyze  the 
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dynamics  of  chain-driven  s)sicms  has  been  presented.  A 
model  of  the  chain  is  formulated.  It  uses  the  equivalent  four- 
bar  linkage  to  derive  a  displacement  function  which  represents 
the  polygonal  effect  of  chain.  This  function,  in  turn,  serves  as 
the  kinematic  constraint  which  is  used  as  an  input  into  the 
DADS  to  be  coupled  with  the  system  equations.  Using  this 
approach,  the  dynamic  analysis  of  a  complicated  chain-driven 
mechanical  system  is  made  possible.  As  an  illustration,  an 
automatic  weapon  system  (i.e.,  the  Model  X  gun)  is  modeled 
and  analyzed. 

The  dynamic  responses  of  the  Model  X  gun  are  given. 
Forces  between  components  as  well  as  motion  characteristics 
of  this  weapon  system  are  fully  explored,  especially  those  for 
the  impact  forces  between  various  components  and  the  chain 
dynamic  loads,  which  are  difficult  to  estimate.  Information 
such  as  this  can  be  utilized  to  provide  a  durability  study  and 
design  optimization  of  the  gun. 

The  result  of  this  investigation,  it  is  believed,  has  extended 
our  capability  on  the  automated  dynamic  analysis  to  a  large 
class  of  mechanical  systems  involving  chain  drives. 
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On  the  Workspace  of  Mechanical 
Manipulators 

This  paper  presents  an  analytical  investigation  of  manipulator  workspace.  The  first 
part  presents  the  derivation  of  a  set  of  recursive  equations  in  terms  of  motion  and 
design  parameters  representing  the  workspace.  These  / ormolus  are  basic  for  the 
determination  of  the  characteristics  as  well  as  the  shape  of  workspace  The 
remuiniler  of  the  paper  concerns  the  investigation  ot  the  existence  of  hole  and  void 
in  workspace.  .1  set  of  criteria  is  formulated.  Algorithms  Tor  implementing  these 
criteria  arc  developed.  Manipulators  with  both  unlimited  and  limned  rerolute  joints 
are  studied  in  this  investigation. 


1  Introduction 

Otic  of  the  basic  questions  encountered  in  manipulator 
design  is  the  determination  of  ihe  shape  of  workspace.  Given 
the  structure  of  a  manipulator,  can  one  specify  the  shape  of  its 
workspace  and  study  its  characteristics  analytically?  The 
problem  is  a  difficult  one,  primarily  because  of  the  large 
number  of  degrees  of  freedom  involved  in  the  manipulator 
system  and  the  inherent  complexity  on  spatial  geometry. 
Therefore,  very  few  investigations  of  the  subject  are  on  the 
record  |l  -6|  Roth  [1|  presented  a  pioneer  study  in  1175  on  the 
relationship  between  the  kinematic  geometry  and  manipulator 
performance  including  workspace.  Kumar  and  Waldron  (2,  3| 
presented  a  thoughtful  investigation  on  workspace  and 
developed  the  theory  and  algorithm  for  tracing  the  bounding 
surfaces  of  manipulator  workspace.  Sugimoto  and  Duffy  |4| 
proved  the  theories  for  determining  the  extreme  distances  of 
manipulator's  hand.  Tsai  and  Soni  15)  studied  the  accessivc 
region  of  manipulator  with  two  and  three  resolute  joints  in 
planar  motion.  Recently,  Gupta  and  Roth  |6j  in  their 
illuminating  paper,  further  extended  the  work  of  Roth  and 
investigated  certain  design  considerations  on  workspace, 
including  the  effect  on  hand  sire.  Basic  concepts  and  some 
investigations  of  primary  (or  dextrous)  and  secondary- 
workspace  have  been  presented  |2,  3,  6|. 

This  paper  and  a  companion  paper  (7|  represent  an  attempt 
to  treat  the  foregoing  question.  The  objectives  of  the  present 
paper  arc:  I )  to  derive  a  set  of  recursive  equations  to  represent 
the  workspace  analytically;  2)  to  formulate  a  set  of  criteria  to 
investigate  certain  characteristics  of  this  workspace,  namely, 
the  existence  of  hole  and  void;  3)  to  develop  an  effective 
algorithm  for  implementing  these  criteria;  and  4)  to  extend  the 
aforementioned  studies  to  the  practical  case  ot  limited 
revolute  joint. 

2  Manipulator  Description 

A  manipulator  with  n  unlimited  revolute  joints  in  scries  can 

Conlnhuied  hv  the  Mechanisms  Commiuee  and  presented  .it  the  Design  and 
Production  I  ngineermg  technical  Conference.  Washinglon,  tl  t  .  September 
12  IS.  I*W2  ol  tin  Ami  an  an  Sorrily  oi  Miihaniiai  I  ci.ini  i  as 
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be  represented  schematically .  as  show  n  m  I  ig.  1 .  1  here  ate  n 
coordinate  frames  to  specify  the  configuration,  for  any 
coordinate  frame,  sav  ft.  the  /.  -axis  is  always  the  joint  axis 
and  the  A,  -axis  is  in  the  direction  of  common  normal  between 
axes  /(  and  .  , .  I  he  /,  -axis  is  attached  to  a  fixed  trame. 
The  last  link,  link  n.  is  associated  with  the  hand  or  the  end 
effector  of  the  manipulator.  Any  link,  say  link  k.  is  always 
considered  to  be  connected  to  joint  k.  Three  parameters  are 


r  ' 


Fig.  1  A  manipulator  with  n  revolute  joints  In  series 


Fig.  2  Geometrical  relationships  between  Joints  n  and  n  -  1(e„  ,  * 
common  normal,  b„  ,  «  aslal  distance.  n„_  ,  =  twist  angle  and  hand 
H) 
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needed  to  speedy  the  geometrical  relation  of  two  consecutive 
joints:  Referring  to  Fig.  2,  they  are  ak,  the  common  normal 
between  two  consecutive  axes  Z*  and  Zk . , ;  ut( ,  the  twist  angle 
of  these  two  axes;  and  bk ,  the  axial  offset  of  joint  A  +  1  on  the 
axis  Zk. , .  For  convenience,  the  hand  twist  angle.  is 
sometimes  considered  to  be  zero.  In  this  investigation, 
however,  it  is  included  as  one  of  the  design  parameters.  In  the 
following,  theory  and  criteria  on  manipulator  workspace  with 
unlimited  revolute  joint  are  first  developed;  they  are  then 
extended  to  the  limited  case  which  is  of  practical  concern. 

3  Analytical  Representation  of  W  orkspace 

The  workspace  of  a  manipulator  is  defined  as  the  region 
which  can  be  reached  by  the  center  of  the  manipulator's  hand, 
II.  Following  the  notion  used  in  |6|,  let  W,  (//)  denote  the 
workspace  generated  by  the  point  //  by  turning  all  the 
revolute  joints  k,  k  +  I,  .  .  .  ,  n  with  complete  rotations  while 
holding  the  axis  k  as  fixed.  Referring  to  Fig.  2,  if  the  revolute 
joint  n  makes  a  complete  rotation,  the  locus  of  the  center  of 
the  hand  H(x„,  y„,  z„)  is  a  circle  with  respect  to  the  coor¬ 
dinate  frame  n.  The  circle  can  be  expressed  in  a  parametric 
form  by 

x„  =  r„  cost9„ 

v„  =  r„  sind,  (1) 

I,  =  b,  cos<«„ 

where 

r"  =  +  (b„  sin<«„ ) :  (2) 

and  d„  denotes  the  rotation  angle  of  the  joint  n.  In  the  case 
that  the  hand  axis  is  parallel  to  the  /-axis  of  joint  n,  i.e.  <«„  = 
0,  equation  ( I )  then  reduces  to 

x„  =  a,  oosd„ 

v„  =  a„  sind„  (3) 

I,  =  b„ 

The  geometrical  relationship  between  any  two  consecutive 
coordinate  frames,  say  n  and  n  -  I ,  can  be  represented  by  a 

_ Nomenclature  _ 


Circular  projection:  It  is  the  cross-section  made  by  the 
physical  profile  of  a  rotating  body  on  a  screen  containing  the 
rotating  axis;  it  depends  on  the  shape  of  the  body  and  the 
geometrical  configuration  of  the  rotating  link,  such  as  the 
twist  angle,  etc. 


Coordinate  frame  A:  Cartesian  coordinate  with  origin  at  joint 
k,  the  /-axis  aligning  along  the  joint  axis,  and  (he  A'-axis 
being  the  common  normal  between  the  axes  k  and  k  +  I 


Space  A:  the  space  in  which  all  points  are  specified  by  the 
coordinate  frame  k 


Ml* 


bk 

PL,  PH 


a  homogeneous  transformation  which 
represents  the  geometrical  relation¬ 
ship  between  the  coordinate  frames  k 
and  k  +  I 

the  common  normal  between  the  axes 

/,  and  Zk , ,  (Fig.  2) 

the  axial  offset  distance  of  the  joint 

*  +  1  on  link  k  (Fig.  2) 

the  distances  between  the  axis  Zk  and 

points  on  W'J,,(/f)  and  !♦'!',,(//) 

(Fig.  7) 

workspace  with  respect  to  (w.r.t.)  the 
joint  A;  or  the  set  of  reachable  points 
of  manipulator's  hand,  //.  w.r.t.  the 
joint  k  (Fig.  3) 


homogeneous  transformation  Ml*-  by  first  turning  the 
coordinate  frame  n  -  1  an  angle  -  <»„  ,  with  respect  to  the 
axis.A’„  i,  then  followed  by  two  pure  translations  along  the 
axes  X „  i  and  Z„  with  magnitudes  <i„  ,  and  b„  , .  respec¬ 
tively,  we  have 

Ml,  = 

I  0 

0  eos«„ 

0  -  sinor,, 

0  0 


0  <i,  i 

sin<«„  |  0 

cos r»„  ,  0 

0  I 


;  i  oo  o 

1  1  0  I  0  0 

i  ; 

too  I  b„ 

*11  0  0  I 


1 

0 

0 

Or  1 

0 

COsr»„  , 

sin<i„  i 

/>„ 

,  sino, 

0 

-sino„  i 

COS<|„  | 

b„ 

(  CO^U, 

0 

0 

0 

i 

(4| 


Let  space  k  be  the  space  that  all  its  points  aic  specified  by  the 
coordinate  frame  k.  Using  equation  (4).  the  circle  represented 
by  equation  (I)  in  space  «  can  then  be  formulated  in  space 
n  -  I  as  follows: 

x  ‘  \  ■  !  r,  cost/, 

I 

!  i  '  I 

v  >  v  r„  smtf„ 

Ml,,  -  Ml,  i  (?) 

z  ■  z  i  b  cosu, 

I  I 

I,,  i  I  „  |  I 


and 


»?(//> 


AX,  A K 
Xk ,  Yk,Zk 


(*i .  y* .  ii ) 


(xk.y„.zk)' 


«i 


A 


0k 

Vi 

»>* 


the  circular  projections  of  »,..(//) 
on  the  left-  and  right-hand  side  of  the 
plane  Xk  Zk .  They  represent  also  left  - 
and  right-hand  side  cross-sections  on 
IT,  (//)  cut  by  the  plane  A'*  /,  (Fig. 
5(uj). 

distances  between  the  axis  /,  and 
points  on  H'i  (//)  and  H'f  (//)  (Fig. 
6) 

(he  coordinate  axes  of  the  frame  A .  7.k 
is  also  the  axis  of  joint  A  (Fig.  1 ). 
the  coordinate  of  a  point  in 
workspace  H*  (//)  w.r.t.  the 
coordinate  frame  A  (Fig.  3) 
the  coordinate  of  a  point  in 
workspace  W\{H)  w.r.t.  the 
coordinate  frame  A  (Fig.  3) 
the  twist  angle  between  the  axes  /, 
and  /, .  i  (Fig.  2) 

(he  location  angle  of  the  link  A  +  1 
w.r.t.  the  link  A;  the  angle  between 
the  axes  A',  and  A  * .  t  in  case  of 
limited  joint  A  +  I  (Fig.  14) 
rotating  angle  of  the  joint  k  (Fig.  12) 
x/2  «», ;  an  angular  quantity  (Fig.  6) 
iotaiioit.il  range  ol  l lie  loinl  A 
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(6) 


f*  i  "  r„  cos*.  -r-a,  i 

y"„  ,  =  r„sint5„cosa„  , 

+  b„cosa„sina„  ,+b„  ,sina„  ,  (6) 

z’  I  =  -r„sini)„sina„  , 

+  bn  cos<*„  cosa„  i  +  b„  |COs«„  , 

The  workspace  W„  ,  (If)  can  be  obtained  by  turning  the  circle 
now  represented  by  equation  (6)  with  respect  to  the  Z„  i  -axis. 
The  shape  of  W„  t(H)  is  a  torus,  as  shown  in  Fig.  3,  and  can 
be  represented  by 

v„  ,  =  /■„  ,cosd„  , 

y„  \  =  r„  ,sini?„  I  (7) 


where 

r"  I  =  ^X„'  |  +  i 

Equations  (I),  (6),  and  (7)  completely  specify  the  evolution 
process  of  generating  a  workspace  from  the  space  n  to  the 
space  n  -  1 . 

A  general  recursive  formula  for  workspace  lik(fh  can 
therefore  be  formulated  as  follows: 

x4  =  r4  cos  if  4 

>’»  =  r4simf4  (8) 

v*  —  Zk 


Jx„-  +  v4- 


X 

• 

* 

y 

=  W«.i 

y 

* 

A. 

z 

1 

k 

1 

rW^’H):  CIRCLE.' 

*H 


V^CH):  TORUS 


•  •  I 

Fig.  3  The  evolution  ot  workspace 

circular  projections  of  the  workspace  H\  .  ,(//)  on  the  plane 
X k7.k.  The  term  “circular  projection"  means  the  projection 
obtained  on  the  plane  XkZk  as  the  workspace  W'4,:<//) 
rotates  with  the  link  k  and  with  respect  to  the  axis  7.k . 

In  the  following  sections,  criteria  are  established  to  examine 
cross-sections  with  regard  to  hole  and  void. 

5  Hole  in  \\  orkspace 

A  workspace  that  has  a  hole  means  there  exists  at  least  one 
straight  line  which  is  surrounded  by  the  workspace  yet 
without  making  contact  to  it.  For  a  manipulator  with  resolute 
joints,  its  workspace  W\(H)  has  a  hole  if  and  only  if  the 
preceding  workspace  IF, . ,  (/f )  does  not  contact  the  axis  /, 
This  implies  that  lift//)  does  not  contact  the  axis/,. 

The  mathematical  criteria  for  the  hole  in  workspace  can  be 
established  based  on  the  minimum  distance  between  ff£f ' " ' 
and  the  axis  /, ,  denoted  by  XH . . .  as  follows: 


1.  Without  hole 


2.  With  hole 


XM .  0 


Equations  (8)  and  (9)  form  the  analytical  representation  of  the 
manipulator  workspace  and  are  the  basic  formulas  for  the 
investigation  of  the  existence  of  hole  and  void  in  workspace, 
and  for  the  determination  of  the  shape  of  workspace  17j. 

4  Cross-Section  of  Workspace  Cuf  by  XZ  Plane 

Having  the  analytical  representation  of  a  workspace,  it  is 
now  possible  to  investigate  its  characteristics.  In  this  paper  we 
are  concerned  with  characteristics  such  as  hole  and  void.  In  a 
companion  paper  (7|,  the  investigations  on  the  shape  of  the 
workspace  including  boundary  and  volume  are  reported.  A 
basic  approach  to  investigate  these  characteristics  involves  a 
close  examination  of  cross-sections  of  a  workspace.  An 
important  and  convenient  cross-section  for  this  purpose  is  the 
cross-section  of  Wk(H)  cut  by  the  plane  XkZk.  This  cross- 
section  can  be  obtained  by  setting  yt  equal  to  zero,  or 
equivalently,  d4  =  0  or  180  deg,  from  equation  (8). 
Mathematically,  the  cross-section  may  be  expressed  as 


z»  =  z;  00) 

When  xk  =  +  rk,  it  means  that  the  cross-section  is  on  the  right- 
hand  side  of  the  axis  Zt,  and  is  denoted  as  Wf  ( // ) . 
Similarity,  the  cross-section  formed  by  jx4  =  -  r4  and  ;4  =  zk 
is  denoted  as  tf^  (//).  It  is  worth  noting  that  and 

W\  ( II)  are  mirror  images  of  each  other  and  represent  the 


.YK„„„>0  (11.2) 

Figures  4(a)  and  4</>)  illustrate  these  cases. 

6  Void  in  W  orkspace 

A  void  is  defined  as  a  region  buried  within  a  reachable 
workspace,  which  however,  is  not  itself  reachable  by  the 
manipulator's  hand.  Ft  out  a  practical  point  of  vieu,  it  is 
always  desirable  to  shun  a  void  in  the  workspace,  except  in  the 
case  where  a  void  contains  a  fixed  pan  of  the  manipulator. 
The  necessary  condition  for  workspace,  IF,  ( //) .  containing  a 
void  is  that  the  preceding  workspace  II  j  ,,(//)  has  a  hole. 
However,  it  is  by  no  means  a  sufficient  condition,  as  shown  in 
Fig.  5.  It  is  found  in  this  study  that  the  shape  of  IF,  ,  ,(lf)  and 
the  design  of  link  k  also  play  important  roles  in  determining 
the  existence  of  void  in  workspace  IF4  ( H) . 

Depending  upon  the  relative  position  ol  the  axis  /4  with 
respect  to  H three  distinct  and  exclusive  cases  arc 
considered.  First,  the  workspace  is  apart  from  the  axis  /4 ; 
this  means  that  the  axis  Z4  neither  contacts  nor  passes 
through  the  hole  of  the  workspace  W'4.,(/f).  Second,  the 
workspace  H'4 .  ,(W)  contacts  the  axis  Zk .  Third,  the  axis  Z, 
passes  through  the  hole  of  workspace  IF, , , .  The  criteria  are 
different  for  these  three  situations  and  they  ate  discussed 
separately  as  follows: 

6.1  Void  in  Workspace:  When  Axis  /,  Is  Apart  from 
Workspace  W  4,  ,(fl).  I  he  condition  of  having  void  in 
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Fig.  4(a)  A  case  ol  Wk{H)  with  hole 


tf,  (//)  in  this  case,  as  shown  in  Fig.  8,  can  be  analyzed  from 
two  aspects:  firstly,  for  given  shape  of  (he  workspace, 
tf, ,  i  <//),  to  determine  the  effect  of  the  twist  angle  a,  and 
secondly,  to  sec  the  effect  of  the  depth  of  workspace  W, , ,(//) 
along  its  axis.  Criteria  for  each  case  are  established.  If  either 


Fig.  5(a)  The  axistanca  ol  void  in  workspace:  W,  ,  t  (H)  has  hols  and 
Wt(H)  has  void 


Fig.  5(6)  Tho  exlstenca  ot  void  in  workspace:  W,  +  ,(M)  has  hole  and 
Wk(H)  has  no  void 

one  of  these  criteria  is  satisfied,  a  hole  in  workspace  If,  , ,  (//) 
will  not  cause  a  void  in  workspace  W  ,  (//) . 

6.1.1  Criterion  Based  on  Twist  Angle.  One  way  to  determine 
the  existence  of  void  in  H (//)  is  by  an  examination  of  the 
cross-section  of  W't„, {H).  By  viewing  from  a  plane  per¬ 
pendicular  to  the  axis  of  the  joint  k,  one  can  examine  the 
relative  clearance  between  the  circular  projections,  namely 
H'J',  |  (//)  and  Wj [.,(//).  This  clearance  depends  upon  the 
twist  angle  a,  and  usually  increases  with  increasing  or, .  If 
(here  exists  no  clearance,  then  there  is  no  void  in  If,  (//). 

For  the  sake  of  illustration,  let  the  cutting  plane  A, ,  i  Z, , , 
for  obtaining  the  circular  projections  Iff. ,(//)  and  W '{ . ,  (//) 
be  perpendicular  to  the  plane  A',  T, .  Referring  to  Fig.  6,  by 
rotating  the  axis  Z, . ,  with  an  angle  y,  =  x.  2  -  a, ,  the  new 
position  Z,' , ,  should  then  be  parallel  to  the  plane  A',  K, .  The 
extreme  distances  between  /, , ,  and  Iff, ,(//)  and  , ,(//) 
are  shown  as  ,Y Kmm  and  A respectively,  as  shown  in  Fig. 
6.  It  is  clearly  seen  that  if  the  difference  between  these  two 
quantities  vanishes,  there  would  be  no  void  in  W,(/f).  We, 
therefore,  have  the  following  criterion: 

Criterion  I:  The  workspace  Wk(H)  has  no  void  if  the 
orthogonal  projection  o/W[ .  ,(H)  and  Wf, ,  (H),  with  respect 
to  a  plane  perpendicular  to  the  7.i , ,  axis,  overlap.  More 
specifically,  the  difference  between  the  extreme  distances. 
XKmin  and\lm„  is  less  than  or  equal  to  zero,  i.e., 

A/?m,n-Af.m>><0  (12) 

Mathematically  it  is  not  necessary  that  the  cutting  plane  A, .  | 
Z, ,  |  be  perpendicular  to  the  plane  A,  V, .  This  is  because  the 
workspace  IT, . ,(//)  is  always  symmetric  with  respect  to  axis 
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Fig.  6  Illustration  ot  the  effect  of  the  twist  angle  on  void 


characteristics  ol  the  shape  ol  li ',  .  ,(/r).  li  affects  the  relative 
distance  of  the  W'f, \(H)  and  with  respect  to  the 

joint  k ,  i.e.  Z, -axis,  which  in  turn  affects  the  existence  of  void 
on  Wt{H).  These  relative  distances  are  discussed  as  follows. 

Referring  to  Fig.  7,  let  the  axis  .V,  , ,  he  perpendicular  to  the 
plane  A,  Z,  and  PR  Ip)  and  /’/.(</)  be  the  distances  ol  the 
axis  Z,  with  respect  to  points  p  and  q  on  W, '.,(//)  and 
H'l ,,(//),  respectively.  Denoting  PRmm  -  PR  [a)  and 
=  PLlb)  as  the  minimum  and  maximum  values  of  PR  ip) 
and  PL(q),  respectively,  we,  therefore,  have  the  following 
criterion: 

Criterion  2:  The  workspace  W,(H 1  has  no  void  if  the 
circular  projections  of  Wf. ,  and  WJ.,(H),  on  any  plane 
containing  the  rotating  axis  Z, ,  overlap.  More  specifically, 
the  difference  between  the  extreme  distances  PR,,  .,  and  PL,,,., 
is  less  than  or  equal  to  zero.  i.  c. , 

PR . s (is. 

The  distance  PRix.  y)  and  PL  lx,  y)  can  be  formulated  as 
follows: 

PRlx.y)  =  \  v:  +  (use  v,  .  ,  -  r, .  , ) 

Pl.(x.y)  -  v  v;-  +  (use  v, . ,  -  r.  ,  , )  (161 


7 


Fig.  7  Illustration  of  the  effect  ol  the  workspace  depth,  d»  +  ,  on  void 


Z, ,  | .  The  criterion  is  true  whenever  the  cutting  plane  contains 
the  axis  Z, ,  | . 

Analytical  expression  for  XRmm  and  AX„,JV  can  be  derived 
as  follows: 

XR„„„  =  Min|r, . , cosy>,  —  c* .  ,  sin*,  ] 

XL,„t ,  =  Max{  —  /•*,! cos*,  —  2* , » sin*,  |  (13) 

where 


where 


I  .<■* 


1  1 


6.2  Void  in  Workspace:  When  Axis  Z,  Contacts  In 
Workspace  II  ,,  ,(//>.  l'lieie  are  iwo  types  ol  contact  between 
the  axis  Z»  and  workspace  M', ,,(//):  single-portion  contact 
(Tig.  9)  and  double-portion  contact  dig.  10).  The  ter¬ 
minology  "double-portion  contact"  means  that  the  contact 
between  the  axis  Z,  and  the  workspace  II', .  ,(/f)  has  a  range 
of  discontinuity. 

In  the  case  of  single-portion  contact,  the  criteria  of  no  void 
in  workspace  II',  (//)  are  the  same  as  the  previous  case  of  no 
contact.  Therefore  when  either  equation  (12)  or  equation  (15) 
is  true,  there  is  no  void  in  II  ,  (//).  In  the  case  of  double¬ 
portion  contact,  it  is  obvious  that  a  void  always  exists  in  the 
workspace  >♦',(//)  regardless  of  whether  the  criteria  of 
equations  (12)  and  (I?)  are  satisfied  or  not. 

6.3  Void  in  Workspace:  When  Axis  Z,  Passes  Through 
The  Hole  of  Workspace  H  , . ,(//).  An  example  of  this  case  is 
shown  in  Fig.  11.  Apparently  there  is  always  a  hole  in 
workspace  H',  (//);  however,  there  is  no  possibility  of  having 
a  void. 


±r.  .1 

0 


6.2.2  Criterion  Based  on  the  Depth  of  Workspace  W, ,  ,(H) 
Along  its  Axis.  The  depth  of  the  workspace  IT, ,  ,(/f)  along 
the  Z, ,  ,-axis  is  defined  as  dk, , ,  as  shown  in  Fig.  7.  This 
quantity  represents,  however,  only  one  aspect  of  the 


7  An  Algorithm  for  Delecting  Void  in  Workspace 

There  are  three  distinctive  geometrical  relationships  be¬ 
tween  the  axis  Z,  and  the  workspace  IF, . ,(//)  as  described 
previously.  Consequently,  corresponding  to  each  category, 
the  cross-section  of  »’,.,(//)  cut  by  the  plane  A',  Z,  are 
different,  as  shown  in  Figs.  8,  9(a)  9(b),  10.  and  1 1 . 

Given  the  structure  of  a  manipulator,  the  question, 
therefore,  is  how  to  classify  it  according  to  these  categories.  A 
mapping  technique  is  presented  in  this  investigation  for  this 
purpose  and  an  effective  algorithm  is  developed  to  implement 
the  technique. 

7.1  A  Mapping  Technique.  Realizing  that  the  circular 
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Fig.  10  AxisZ,  contacts  to  W,  »  ,(M):  double  portion  contact 
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Zk 

t  SECTION  B-B 
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Fig.  9(a) 

Section  A-A;  cut. through  solid 


Fig.  9(b) 

Section  B-B;  cut-through  hols 


Flg.1t  Asians  passed  through  the  hole  of  workspace  Wk  *  t(H) 


projection,  either  or  H'J. ,(//).  can  be  readily  Xt , ,  7.k . , -plane  to  the  plane  XtZk .  Referring  to  Fig.  12,  all 

obtained  (equation  (10)),  the  problem  now  is  to  find  the  points  on  H'f,  ,</f)  are  mapped  onto  the  plane  Xk  ?.k  through 
appropriate  procedure  to  map  this  cross-section  from  the  a  rotation  of  link  A  t  I  at  various  angles  denoted  by  d, , , . 
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This  angle  can  be  derived  using  a  homogeneous  trans¬ 
formation  formula,  equation  (9),  and  then  setting  v*  equal  to 
zero,  we  have, 

„  .  ,/  ,sin«A +b*sin<«M 

0k .  i  =sin  ( - )  (18) 

\  rk,  |  cos<*4  > 

and 

*»•  =  r, .  icosd*.  i  +fl» 

y\  =  rA,|Sini>A,ieosa*-t-zA,|Sina*  +  bAsinaA  =0  (19) 

z*  =  -  rk ,  |Sint)A ,  ,sinaA  +z4 ,  ,cos<«A  +  />4cos«A 

Equation  (18)  shows  that  different  points  (r, , , ,  zA.,)  on 
M'f.i(W)  have  different  mapping  angles  d4 .  i .  Knowing 
ftk . , ,  the  corresponding  point  on  A',  Zk  -plane  in  space  k  can 
be  calculated  from  equation  (19). 

7.2  A  Numerical  Algorithm.  An  algorithm  is  developed  to 
analyze  data  obtained  from  the  mapping  function,  equations 
(1 8)  and  (19).  It  is  worth  noting  that  equation  (18)  is  not  single 
valued;  lor  each  point  of  (r4 , zA .  i),  there  correspond  two 
values  of  it, . , ,  namely,  dA , ,  and  its  complementary  angle, 
*-  i )A , , .  Consequently,  two  values  of  xA  are  obtained  from 
equation  (19).  This  results  in  two  groups  of  points:  |.rA  (*?*., ), 
zA  |  and  [jrt  ( *  -  dA ,  i ),  zA).  From  the  sign  of  the  extreme 
values  of  .rA  for  each  group  denoted  by  (max(x, ),  min(.v,)| 
and  [maxU:),  minU2 )],  classifications  of  different  categories 
based  on  general  relationship  of  Zk  and  IV A . ,  (ff)  can  be 
established  as  follows: 

I .  Axis  Zk  is  apart  from  IVA  ,,(tf): 

Sign|max(jr, ))  =  Sign[min(jr,  )| 

=  Sign[max(x2  )|  =  Sign(min(x2  )|  (20) 

2A.  Single-portion  contact: 

SignfmaxU,  )|  =  Sign|min(jr, )) 


Sign|max(.v- )]  z=Sign|min(A  •  )| 
or 

Sign[max(.r,  )J  *  Sign[min<» , )) 

Sign[max(.v2))  --  Signjmimv.  )|  (21 1 

2B.  Double-portion  contact: 

Sign|max(.f|  )1  *  Sign|min(x ,  )J 
Sign|max(.v:  )]  *  Sign(min(A':  )| 

3.  Axis  Zk  passes  through  the  hole: 

Signlmaxt  r, ))  -  Sign[min(  v, )) 

*Sign(max(.v2)|  -  Sign[inin(.v:))  (23) 

7.3  Special  Conditions.  There  arc  two  special  conditions  of 
equation  (18):  I)  When  the  numerator  is  greater  than  the 
denominator,  there  will  be  no  solution;  and  2)  when  cos,,, 
equals  zero. 

This  first  condition  implies  that  mapping  does  not  exist. 
Physically  this  means  that  the  point  in  consideration  does  not 
come  in  contact  with  the  plane  A,/,  when  link  k  t  I  makes  a 
complete  rotation,  such  as  the  point  p  in  I  ig.  9. 

The  second  condition  refers  to  the  case  when  the  twist  angle 
between  the  axes  Zk  and  Zk .  ,  is  of  a  right  angle.  Con¬ 
sequently  the  mapping  is  meaningless.  Referring  to  f  ig.  13. 
the  points  u,  b,  and  c  cannot  be  mapped  onto  the  plane  .V,  /, . 
In  this  case  equation  ( 1 5)  should  be  used  as  a  criterion. 


8  Manipulator  With  Limited  Resolute  Joints 

The  discussions  and  development  so  far  were  based  upon 
unlimited  revolule  joints,  in  which  case  three  parameters, 
namely,  oA,  bt,  and  ak,  are  sufficient  to  specify  the 
geometrical  relationship  of  two  consecutive  joints.  However, 
in  practice,  most  manipulators  have  limited  joints.  In  this 
section,  formulas  and  modifications  are  presented  which  take 
into  account  this  effect.  It  is  found  that  an  additional  angular 
parameter  must  be  introduced.  In  the  following,  the 
justification  of  an  additional  parameter  is  presented,  followed 
by  a  discussion  on  how  to  define  this  parameter;  then  a 
formula  for  the  modified  homogeneous  transformation  ts 


T ransactlona  ol  th«  ASME 


■dt*. 


*. . 


Ml  Vol  10b,  MARCH  1983 


Character isticx  ot  woikv|ia,.c  fivudmg  .Kile  aitJ  c end  tor  a 
limited  joint  become  relatively  simple,  l  or  a  manipulator  with 
a  limited  joint,  say  joint  A+l.  there  will  be  no  hole  in 
workspace  Hj. ,(//);  consequently,  there  will  be  no  void  in 
workspace  H  ,  ( H) . 

9  Conclusion 

This  study  provides  some  basic  theories  and  formulations 
concerning  the  manipulator  workspace.  Analytical 
representation  of  workspace  is  presented.  The  difficult 
problem  of  the  existence  of  hole  and  void  in  workspace  can  be 
treated  with  a  set  of  procedures  and  criteria  established  in  this 
investigation.  Practical  considerations  arc  also  included 
Specific  examples  illustrating  the  basic  theory  presented  here 
are  reported  in  a  companion  paper  (7|,  together  with  an  in¬ 
vestigation  of  the  boundary  and  volume  of  workspace.  The 
result  of  this  research,  it  is  hoped,  will  contribute  toward  the 
basic  understanding  and  design  of  mechanical  manipulators. 


derived;  and  finally,  the  characteristics  of  hole  and  void  are 
discussed. 

For  an  unlimited  joint,  the  Cartesian  coordinate  frame 
associated  with  this  joint  is  defined  by  aligning  the  /-axis 
along  the  joint  axis,  and  by  specifying  the  A-axis  as  the 
common  normal  between  two  consecutive  axes.  Since  the  joint 
can  perform  unlimited  rotations,  the  location  of  .V-axis  on  the 
plane  A, . ,  V, . ,  can  be  arbitrary.  This  is  not  true,  however, 
for  a  limited  joint.  Since  there  is  a  portion  of  turning  range, 
physically  it  cannot  be  reached  by  the  link;  the  location 
defining  this  link  therefore  becomes  important.  Additional 
constraint  on  the  selection  of  the  A',  -axis  must  be  introduced. 

Considering  a  particular  joint  which  is  limited,  say,  the 
joint  In  I,  it  is  now  relevant  to  define  the  Cartesian  coor¬ 
dinate  frame  associated  with  this  joint.  Referring  to  l  ig.  14, 
by  specifying  the  link  k  +  1  which  is  always  in  the  mid-range  of 
the  joint  rotation,  and  consequently  the  A',., -axis  is  also 
defined.  A  new  parameter,  d,  is  introduced  and  is  defined  as 
the  angle  between  the  A, -axis  and  the  A, .  ,-axis  in  the  same 
direction  as  the  right-handed  screw  with  respect  to  the  axis 
/,  .,.  Ilus  parameter  specifies  the  physical  location  of  the 
link  k  *  I  with  reference  to  the  link  k  and,  therefore,  is  called 
the  location  angle. 

The  homogeneous  transformation  between  the  coordinate 
frames  k  v  I  and  k  can  now  be  represented  by 
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On  the  Evaluation  of  Manipulator 
Workspace 

This  paper  presents  a  theorem  regarding  manipulator  workspace  anti,  hayed  on  the 
theorem,  a  manipulator  pertormanee  index  is  introduced  li  is  tound  that  tor 
given  manipulator  structure  the  ratio  ot  the  eo/nme  o'  the  work  space  to  the  cube  ot 
its  total  link  length  is  a  constant.  Algorithms  lot  outlining  the  houndarx  ptot-.'.c  ot 
workspace  anti  for  quantitative  eradiation  ol  us  i  ohmic  arc  presented  I  compute r 
package,  KAM.  is  developed,  whnh  implements  the  theories  and  algorithms 
developed  in  this  investigation  as  sxcil  ,o  m  a  ,  i •".•/•anion  paper  |  / 1  S*'i  era I  spud  a 

examples  are  given  to  illustrate  the  app.'u  a: .  p,  r  unman,  e  index  ana  //», 

capability  ol  KAM. 


1  Introduction 

There  are  a  few  basic  problems  on  manipulator  workspace 
remaining  unsolved  or  only  partially  solved  1 1  -4] .  For  in¬ 
stance,  what  criterion  on  workspace  can  one  use  on  the 
evaluation  of  (he  kinematic  performance  of  manipulators,  ot 
more  fundamentally,  how  can  the  shape  of  a  manipulator's 
workspace  be  determined? 

This  paper  deals  with  these  problems.  In  particular,  the 
objectives  are:  firstly  to  develop  an  algorithm  to  outline  the 
boundary  of  the  workspace  and  to  evaluate  its  volume 
quantitatively;  and  secondly,  to  present  a  manipulator  per 
formance  index  which  is  based  on  workspace  and  can  Ik 
evaluated  efficiently. 

2  Manipulator  With  Kevolule  Joints 

An  analytical  investigation  of  manipulator  workspace  i' 
given  by  Yang  and  I.ee  |l|.  The  following,  taken  tiom  |l), 
represent  a  description  of  manipulator  gcomcitv  and  a 
summary  of  the  results  on  the  analytical  representation  ot 
workspace. 

A  manipulator  with  n  unlimited  resolute  toiiils  m  senes  can 
be  represented  schematically,  as  shown  in  Fig.  I .  I  lictc  are  n 
coordinate  frames  to  specify  the  configuration.  I  or  any 
coordinate  frame,  say  *,  the  /,  axis  is  always  the  joint  axis 
and  the  \\  axis  is  in  the  direction  of  the  common  normal 
between  axes  /,  and  /, , , .  The  /,  axis  is  attached  to  a  fixed 
frame.  The  last  link,  link  n,  is  associated  with  the  hand  or  the 
end  effector  of  the  manipulator.  The  link  k  is  connected  to 
joint  k.  Three  parameters  are  needed  to  specify  the 
geometrical  relation  of  two  consecutive  unlimited  joints  (joint 
being  capable  of  making  a  complete  rotation).  Referring  to 
Fig.  2(a),  these  parameters  are  a,,  the  common  normal 
between  two  consecutive  axes  Z,  and  the  twist  angle 
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Fig.  2(a)  Geometrical  relationships  between  unlimited  revolute  joints 
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Ml. 


LOV.i. 

com,  vm.i, 

vino,  void, 
0 


snip, 

0 

cost* 4  COvff* 

sinnrj 

/>,S||1<|, 

sirwi^costfi 

COS  U  4 

/tjl'OSo 

0 

0 

1 

where  d, ,  referred  to  av  the  location  angle,  denotes  the  angle 
between  the  axes  ,V,  and  A", , ,  while  the  joint  k  i  1  is  located 
at  the  center  of  its  rotational  range. 


3  Boundary  and  Volume  of  Workspace 


Having  an  analytical  representation  of  a  workspace,  it  is 
now  possible  to  investigate  the  boundary  and  volume  of 
manipulator  workspace  quantitatively.  The  basic  approach 
involves  the  construction  of  a  proper  cross  section  of  the 
workspace,  and  from  which  numerical  methods  are  developed 
to  determine  its  boundary  and  subsequently  the  volume  ot  the 
workspace. 

Following  (I |,  an  important  and  convenient  cross-section 
for  this  purpose  is  the  cross-section  on  IF,  (//)  cut  by  the 


plane  V, /,.  It  van  be  ob'.a  n„.;  !  v  ..  I  ..  ..  ,.  .qua'.  :o  <ti.  .  .•! 
equivalently,  if,  -0  oi  t Si i  dec,  Horn  j-quaiior.  Mi 
Mai  hemal  ically,  this  cross  million  may  be  expressed  as. 

V .  ir. 


When  v,  =  +  r, .  n  implies  dul  the  cross. iivt.o"  i«  on  tlie 
right-hand  side  o  -tic  axis  /,  .  and  i-  dci'.oicd  .o  H.1//1 
Similarly,  the  cross  section  toi  mid  b>  1,  \  and.-.  „ 

denoted  as  H  (  (//).  It  is  worth  noting  1  ri.it  H't/fi  and 
U[(H)  are  mirror  images  ot  each  other  and  icprc'cm  r t,e 
circular  projections  of  the  workspace  it  .  .  ill  1  on  the  plane 
X,  /, .  The  term  "circular  protection"  |i|  mean'  the 
projection  obtained  on  the  plane  A./.  .0  tie  woihsp.ui 
H , .  ,  (//)  rotates  with  the  link  k  and  with  ic'pic:  to  the  axis 
/,. 

Because  the  workspace  I  *  .  1  /  /  >  n  wnmietn,  wuhic'pi  t  to 
the  axis  /, .  It  '  (I/I  contains  all  the  gconicm, al  mlortiiat.on 
regarding  the  shape  of  the  wotkspace.  I  lie  mathematical 
expression  of  It  ft//),  in  equation  (5i.  can  be  tewmten  :u 

parametric  form  as. 

•v,  =r,  =/(!»,.,. H,.: . th  I 

:,  =;T  =m  if. . . .if. . . fS ,  >  (fo 

By  scanning  sequentially  the  ioim  angles  ,<  .  o.  ,  and 
d, ,  1  according  to  their  rotation  limits,  at  a  lej'Oiubly  small 
interval,  most  of  the  points  ot  It  .*  ( //)  can  thus  he  obtained 

3.1  Boundary.  An  nlgoiiilim  lor  diiiMing  ilie  bonn.la’v 
of  the  workspace  is  developed  It  is  given  as  lollows 

1.  Dividing  the  plane  A./,  into  /  x  ./  small  rectangles  ot 

width  Air  and  height  A.-,  where  /  and  ./  .tie  the  numbei  ot 
divisions  along  •  lie  V  axo  and  ilie  /  avis,  ic-pc.inclv  I  a, l; 
rectangle  is  denoled  as  /'  (  ,  1  winch  is  used  to  puw.de  a 

binary  imagine  of  II (//1. 1  e  . 

/*  I .  if  /’  is  a  portion  of  H  t  H) 

/'„  -  0.  otherwise. 

2.  Initiali/alioti:  set  all  /*  )  0. 

3.  Scanning  process; 

tu)  Scanning  all  ifx  beginning  with  if... 

t/>l  For  any  d.  say  if. .  scanning  up  to  O..  with  step  ,i/(  a,* 

4.  Construction  of  the  binary  image  of  it  ;  t  // ) 

P,  t.v  I ,  if  any  point  (\,;i  in  equation  (b)  falls  into  the 
following  range  |.v  -A\  ’.x  *Ai  2]  and  (r  Ae  2.* 
+  Ac  21  where,  Av  =  v . .  1  andA:  =  C.  . 

Record  the  extreme  v allies  ol  within  each  uv (angle.  P 
I  his  gives  ;  ,„„(v  )  and  ( x  J 

(do  io  step  3  until  the  scanning  pioccss  is  ,omple(c  A  data 
base  for  H'f  [//)  is  now  stored. 

5.  Construction  of  the  boundary  of  lift  IP 

The  extreme  points  of  those  rectangle'  which  have 
/’,  (v,,c,)=l  and  are  adjacent  to  the  lectangles  with 

( v,,:,)  =  0.  define  ihe  boundary  of  H‘2(//l  (referring  to 
F'ig.  3).  This  construction  is  done  by  a  scanning  process  which 
searches  all  /  x  J  rectangles  systematically 

This  algorithm  is  capable  of  determining  the  boundary  ol 
workspace  with  irregular  shape  having  eons  exit  v  and  void  If 
»!■(«)  is  concave  or  it  has  void  with  respect  to  the  /  axis  at 
certain  values  along  the  A  axis,  say  \  .  then  theie  is  a 
discontinuity  of  the  value  of  P  ,  i.c..  P  changes  trom  1  to  0 
and  returns  to  1 .  This  characteristic  can  also  be  seen  front  Fig 
3. 

3.2  Volume.  Knowing  the  boundary  of  ll'ft//!,  the 
volume  ot  w 01  k space  It  .  (II)  ,  .111  then  be  .  ah  ul.il ed  In 
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Fig  3  Construction  of  the  boundary  of  wf(H) 


Fig.  4  An  illustration  of  the  integration  of  workspace  volume 


'"‘-T  |  2r<r|(:,„„|jr)  ;„„„(v))Av] 

v*  =  I,  •!  <7> 

2irx[(;,„.u(.v)  Illim,<v)  ♦  c . (OIAx| 

referring  to  f  ig.  4,  where  A'„M,  and  .Vn„„  denote  he  maximum 
and  minimum  of  the  x  value  of  Wf  (//);  (  x )  and  cml„(.v) 

are  extreme  values  of  z  at  .v;  C,ll>ln( v)  and  (x)  are  local 
extremes  of  I  value  at  the  discontinuity  while  Wf(//)  is 
concave  or  has  void  at  x;  and  Ar  represents  the  width  of  the 
dividing  rectangle.  I  or  a  certain  value  of  x,  if  H  f  (//)  has  no 
discontinuity,  the  upper  equation  of  equation  (7)  is  used;  if 
there  exists  discontinuity,  the  lower  one  should  he  used.  If  the 
joint  k  is  a  limited  revolute  joint,  then  the  angle  2*  in 
equation  (7)  should  be  replaced  by  the  physical  turning  range 
O,  of  the  joint.  However,  in  this  situation,  error  may  occur  in 
equation  (7)  from  the  undetermined  end  shape  of  the 
workspace,  fortunately,  for  most  commercially  available 
robots  the  first  two  or  three  links,  which  iisunllv  determine  the 
volume  of  the  workspace,  arc  coplanar  ’herefore  the  error 


Fig.  5  Example  1  -  boundary  outlined  by  K  AM 


rxxr^Lt  i  2 


Fig.  6  Example  1  -  workspace  with  void 


due  io  the  undetei mined  end  shape  will  merely  he  caused  by 
the  Iasi  less  short  links,  and  shall  not  he  seis  sigmticanl  in 
practice. 

4  Computer  Program  and  Kxamples  for  the  Boundary 
and  Volume  of  Workspace 

A  computer  program  called  KAM  (kinematic  Analysis  of 
Manipulators)  and  written  in  l-OR  IRAN  is  developed  I  lie 
code  implements  the  theories  and  critcna  presented  in  this 
paper  as  well  as  in  the  companion  paper  [1]  l  or  a  given 
manipulator  structure,  the  progiam  first  defines  analvtuully 
the  workspace,  then  explores  the  existence  ol  hole  and  void, 
outlines  the  boundary,  and  finally,  calculates  the  total  volume 
of  the  workspace.  The  inputs  of  the  program  arc  the 
geometrical  parameters  specifying  the  manipulator  structure 
These  are,  the  number  of  revolute  joints,  the  common  nor¬ 
mals,  twist  angles,  and  the  axial  offset  distances.  Tor  a  limited 
joint,  two  additional  parameters,  the  location  angle  J.  and  the 
physical  turning  range  (.),  of  the  joint  are  needed.  At  present, 
kAM  can  handle  manipulators  up  to  six  revolute  joints; 
however,  the  extension  of  a  higher  number  of  joints  does  not 
present  a  problem.  It  can  be  easily  done  if  it  is  needed. 

Three  problems  which  ate  believed  to  be  representative  arc 
chosen  to  demonstrate  the  effectiveness  as  well  as  the 
capability  of  kAM. 

Ixample  I.  I  firs  example  gises  some  indications  of  how 
well  the  prediction  of  kAM  compares  wiili  the  theoretical 
tcsiilt  on  a  special  ease.  The  problem,  taken  from  Ciupta  and 
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Fig  7  Example  2  -  boundary  outlined  by  KAM 

Roth  1 5 J ,  rulers  to  a  very  simple  manipulaint .  Ii  has  three 
unlimited  resolute  joints  with  special  proportions:  </-  -  10. 
u:  -  4.  u,  I ;  /j,  -  0,  It;  =  0,  b,  -  0;  n,  =  0  deg.  ...  90  deg, 

n,  0  deg.  Because  of  the  simplicity  ol  the  structure, 
analytical  formulas  for  calculating  the  volume  and  the 
condition  of  existence  of  void  in  terms  of  the  angle  u,  can  he 
derived  and  have  been  given  |5|.  A  comparative  studs  is  given 
as  follows: 

volume  |in. 'J  angle  a, 
analytical  |5|  1202.7  u,  >  14.447.5  deg 

K  AM  1215.7  a,  >14.478  deg 

I  he  comparison  looks  very  good.  The  predicted  value  of 
volume  is  within  I  percent  of  the  analytical  value,  and  the 
difference  on  the  angle  nr,  which  is  the  angle  that  the 
workspace  starts  to  have  void,  is  less  than  0.004  percent, 
figure  5  chows  the  boundary  of  the  workspace  of  this 
example. 

figure  6  shows  a  workspace  which  has  a  void  m  it.  All 
design  parameters  are  the  same  for  this  ease,  except  for  the 
twist  angle  nr,  which  equals  20  degrees.  The  volume  ol  this 
workspace  equals  2209.9  in’. 

kxumpU’  2.  In  this  example,  a  comparative  study  is  given 
with  the  method  of  Kumar  and  Waldron  |6|.  A  manipulator 
has  six  unlimited  revolute  joints  with  the  following 
geometrical  configuration:  tr,  72,  a.  16,  r i,  -  8.  a,  4, 
n.  2.  i>„  l;/>,  0,  b.  5,  />,  0.  b,  0,  />,  I.  />„  0; 
ri,  26.7deg,  ri.  127.5  deg,  n,  24.5  deg,  <i4  86.0 

deg.  n,  147  deg,  n.,  -  0  deg.  I  he  boundary  of  workspace 
outlined  by  KAM  is  shown  in  l  ig.  7,  and  matches  with  the 
result  of  (6|.  furthermore.  KAM  provides  additional  in- 
t  or  mat  ion  such  as  no  void  in  workspace  and  a  quantitative 
value  of  the  volume  ol  workspace  which  equals  476757  in' . 

kxumple  J.  In  this  example,  a  commercially  available  robot 
PUMA  600  is  used  as  an  illustration.  Of  particular  interest  is 
Ibal  this  is  a  practical  example,  like  any  industrial  robot, 
usually  having  limited  revolute  joints.  Design  parameters  of 
this  manipulator  are:  cr ,  0,  u,  18,  a ,  0.  a4  0,  <r,  0, 

</„  -  0;  />,  -  10.  />,  =4.  6,  ^  1 7,  bi  -  0,  6,  2.  b,,  0; 

■  i,  -  90 deg,  <»,  -  I80deg,  u,  =  90  deg,  •»,  -  90  deg,  u,  90 
deg.  He  0 deg;  & ,  =  -90  deg,  0.  -  -  90  deg,  d ,  ^  20  deg, 
d4  IHOdeg,  0,  =0  deg;  O,  =  220  deg,  ();  - 250  deg, 
O,  -  270  deg,  04  =  280  deg,  O,  -  200  deg,  0„  -  572  deg.  Since 
the  hand  of  a  manipulator  is  interchangeable,  depending  on 
applications,  it  is  best  to  exclude  it  from  our  evaluation  of  the 
workspace.  The  manipulator  is,  therefore,  equivalent  to  a 
live  joint  manipulator  Using  KAM,  the  win  k space  of 1*1  IMA 


8  !8  .'t  38  *  V 

Fig  810)  PUMA  600  -  real  case 
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F fg.  8ft>)  PUMA  600  -  ideal  case 


600  is  found  to  be  1 87885.6  in. ',  and  its  boundary  is  shown  in 
I'ig.  8(a). 

To  make  a  comparison  with  an  ideal  ease,  that  is.  if  all  the 
joints  on  the  I’lJMA  600  were  unlimited,  the  volume  of  the 
workspace  then  equals  to  216410  in.',  and  its  boundary  is 
shown  in  fig.  8(/>).  It  is  worth  noting  that  in  this  example  the 
workspace  of  the  real  manipulator  is  about  17.6  percent  less 
than  the  ideal  case. 

5  A  Performance  Index  Rased  on  Workspace  Volume 

Given  a  manipulator  we  are  now  capable  of  evaluating  its 
workspace  volume  quantitatively.  The  question  remaining  in 
this  investigation  is  the  development  of  a  performance  index 
based  on  workspace  volume.  In  this  section,  such  index  is 
presented.  A  tlieoietical  proof  is  given  which  shows  this  index 
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is  a  ainst.nu  and  remains  to  be  the  same  lot  sttmlai 
manipulator  structures. 

The  length  of  links  of  manipulator  with  n  resolute  joints 

are  a, ,  a: . .  a„  and  />, .  h:, . I>„  which  correspond 

to  their  common  normals  and  axial  offsets.  I  or  convenience 

in  our  discussion,  they  are  renamed  as  . /„, 

...  Based  on  the  concept  of  dimensional  analysts,  the 
solumc  ol  this  manipulator's  workspace  can  he  represented 
by  a  general  formula. 

V=  <X) 

where  i.j.  and  k  can  be  any  integer  from  I  to  2//,  and  is 
the  corresponding  coefficient  of  the  term  /, /, /* .  The  coef¬ 
ficients  tk,,L  can  be  any  real  value  depending  on  the  structure 
of  a  manipulator. 

It  is  found  in  this  investigation  that  there  is  a  definite 
relationship  between  the  volume  of  workspace  and  the  total 
length  of  a  manipulator.  This  is  given  in  the  following 
theorem: 


Theorem,  l  or  a  given  manipulator,  (i.e..  the  twist  angles, 
rotation  limits,  location  angles,  and  the  proportions  of  link 
length  are  given),  the  ratio  between  the  volume  of  the 
workspace  and  the  cube  of  thetotal  length  is  a  constant,  i.e.. 


V 

(/,  +7;  +  ....  f 


V 

-j  =  constant 


(4) 


Proof: 

Since  the  proportions  of  link  length  are  fixed,  the  length  of 
any  link  can  be  represented  by  the  length  of  another  link  with 
the  given  proportion.  Now  let  every  link  be  in  terms  of  the 
axial  offset  distance  of  the  last  link.  /.„ ,  we  have. 


l\  =  P  i 

I :  =  /’:  /> 

.  .  .  (10) 


,  =  />>  i 

where /> 's  are  known  proportions. 

Combining  equations  (8)  and  (10), 

V  y,n„il,l:h  =  £ </>,/>,/’!>- in/'.,  (II) 


where  1/’|  -•  a/’ /’./’*  =  constant. 

t./.A 

Let  /  represent  thetotal  link  length,  then. 

/.  =  /,♦/..*  ....  *  /.„ 

=  (/'i  *■/>:*  ...  .  *  />:„  ,  <  It  (Id) 

- 


where  |yi  (/’i  t-  /';  ♦  .  .  .  t  p:„  ,  t  1)  -  constant. 
hrom  equations  (II)  and  (12),  we  have 


v  =  =  i  n 

i'  \Q\'il  iyr 


(13) 


Since  |P)  and  |Q|  arc  constants,  VI  is  also  a  constant.  This 
proves  the  theorem. 

The  notation  VI  stands  for  the  Volume  Index  of 
manipulator  workspace,  which  gives  an  indication  of  the 
effectiveness  of  link  length  on  the  creation  of  reachable 
workspace.  The  index  can  be  normalised  by  dividing  its 
possible  maximum  value,  l.ctting  the  total  length  of  a 
manipulator  equal  L,  the  maximum  workspace  that  this 
manipulator  can  possibly  have  is  a  sphere  with  radius  /.  and 
centered  at  the  joint  I.  Therefore,  the  maximum  VI  a 
manipulator  can  possibly  have  is. 


4  tl  7.3 

i: 


4ir 

3 


»4. 1888 
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Fig.  9(a2) 


Fig.  9(a)  Structure  and  boundary  ot  PUMA  600 


Using  this  value  as  a  normalizing  factor,  a  normalized  volume 
index,  denoted  as  NVI,  can  be  obtained,  and  its  value  is 
between  0  and  1. 

This  index  is  a  kinematic  performance  index  which  can  be 
used  in  practice  for  quantitative  evaluation  of  manipulators 
based  on  workspace.  Using  the  algorithm  developed  in  this 
investigation,  this  performance  index  can  be  evaluated  ef¬ 
ficiently. 
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FiB.  9(62) 

Fig  9(6)  Structure  and  boundary  of  Cincinnati  Milacron  T  J 
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Fig  9(c2) 

Fig.  9(c)  Structure  and  boundary  ot  Pane  Bobo 


A  Case  Studies:  The  evaluation  of  Manipulators 
Rased  on  Volume  Index 

Five  commercially  available  manipulators,  as  shown  in  hi):. 
9,  are  evaluated  in  terms  of  VI  and  NVI.  These  manipulators 
are: 

Case  I:  PCMA  600  J7|.  The  design  parameters  of  this 
manipulator  are  given  in  example  .1  of  Section  4.  The  structure 
of  this  manipulator  and  the  boundary  of  the  workspace  from 
KAM  arc  given  in  Tig.  9 (u). 

Case  2:  Cincinnati  Milacron  r'(8|.  The  design  parameters 
of  this  manipulator  are:  a,  <l,u.  40,  a,  40,  </,  S.u.  O, 

u„  0;  />,  0.  h.  <),/;,  <>.  b, -0,  b,  13.  />„  0;  9t) 

deg.  <r:  -  Odeg,  <r,  =0  deg,  <»4  =90  deg,  nr,  =  90  deg,  n>  =  0 
deg;  di  60  deg.  d;  =90  deg,  d,  =0  deg,  d4  --  90  deg; 

O,  -  240dcg,  O.  =  120  deg.  O,  =  120  deg,  O,  120  deg, 
(>4  -  1 80 deg.  O.  =  270 deg,  0„  =  360 deg.  The  structure  of 
this  manipulator  ami  the  boundary  ol  the  workspace  are  given 
m  I  tg  91  h) 


Case  .7:  Pana-Kobo  |9).  The  design  parameters  of  this 
manipulator  are:  a,  =0.  ti,  =  23.622,  «,  31.496,  u,  (I. 

u,  =  0;  6,  =  0.  b:  =  0,  />,  =  0.  b4  =  4.724.  h.  0;  <> ,  --  90  deg. 
no  -  Odeg,  €»<  =0  deg,  or4  90  deg,  o,  =0  deg:  ,i ,  •-  82.5 

deg.  d;  =  97.5  deg.  di  =  90  deg:  O,  -  3(K)  deg.  O;  -  105  deg. 

O,  =75  deg,  04  =  180  deg,  O,  =  360  deg.  The  structure  ot  this 
manipulator  and  the  boundary  of  the  workspace  from  KAM 
are  given  in  Tig.  9(c). 

Case  4:  All)  800  1 101.  The  design  parameters  of  this 
manipulator  are:  o,  =  0,  a,  =23.622,  a ,  31.496.  a,  0. 

u,  0;  b,  0,  />,  0.  ft,  o'  h,  1.9,  h.  0;  <r,  90  deg. 

it,  -  (I deg,  r»,  0  deg,  or4  90  deg,  n.  0  deg;  d,  92  5 

deg,  d;  =  100  deg.  d.  =  90  deg:  (),  =  300  deg,  O.  -  95  deg. 

O,  =  70 deg,  04  =  180 deg,  O,  =  360 deg.  The  structure  of  this 
manipulator  and  the  boundary  of  the  workspace  from  KAM 
are  given  in  Fig.  9(r/). 

Case  5:  Cvbolech  3.70  |ll|  lire  design  parameters  ot  this 
niailipulaloi  are:  <i,  O.u.  7I.496.U,  79.37,  a,  0,  </.  0, 
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Fig.  9W1) 


I  IS  *  «  m 

Fig.  Btttt) 

Fig  9(d)  Struct ure  and  boundary  of  Automatic  AID  900 


Fig.  9(a2) 

Fig.  9(a)  Structura  and  boundary  ol  Cybotech  V30 


b,  =o.  h:  =0.  b,  =0,  b4  =  2J.622,  b,  =0;  «,  =90  deg,  <>  -  0 
deg,  ii,  0  deg,  «4  90  deg,  ,«»  0  deg,  |f,  90  deg, 

H,  =  80 deg,  Ii,---  -90 deg;  O,  =240 deg,  0.-  80  deg,  O,  = 
90 deg,  04  240 deg,  O,  -  >60 deg  [he  structure  of  this 

manipulator  and  the  boundary  of  the  workspace  from  KAM 
are  given  in  Fig.  9(e). 

It  is  worth  noting  that  there  is  a  link  structure  which  is 
common  to  all  the  abovementioned  manipulators.  Thai  is  a 
joint,  say  joint  k  +  I,  located  on  the  /  axis  of  the  previous 
joint,  joint  k  (for  instance,  joints  1  and  2  of  cases  2,  3,  4,  and 
5;  and  joints  4  and  5  of  case  1  belong  to  this  kind  of  struc¬ 
ture).  Consequently,  the  distance  between  these  two  joints 
with  respect  to  the  axis  Z4  does  not  affect  the  shape  and 
characteristics  of  workspace  ff't(W).  Referring  to  Fig.  10,  it 
is  obvious  that  the  distance  d.  only  affects  the  relative  height 
of  the  workspace  tf’4  (//)  with  respect  to  the  coordinate  frame 
k.  Therefore  d.  is  considered  to  be  zero  in  this  comparative 
study. 


Fig.  10  An  illustration  ol  tha  link  structural  the  joint  h  *  1  la  on  the  aria 
of  joint* 
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The  evaluation  of  these  manipulators  based  on  VI  and  NVI 
is  given  in  the  following  table: 

I  able  I  I’crfurmance  evaluation  of  industrial  robots  based 
on  volume  index 


I>.K 

1.  Length 

Volume 

VI 

NVI 

H:\Iamu 

6 

5  to 

18.39X6 

1 .11 

l>.  33) 

t  At.  I  1 

6 

101.0 

2295233 

2.2.1 

0.532 

I'anu  Kobo 

5 

51*4 

355937 

I.M. 

0.  VT 

All)  SOU 

5 

59.02 

267350 

1.30 

0.310 

IT  V30 

5 

94.49 

1299251 

1  54 

0.368 

7  Conclusion 

The  result  of  this  investigation  has  demonstrated  that  it  is 
possible  to  analytically  evaluate  the  manipulator  workspace 
and  most  importantly,  it  presents  a  kinematic  performance 
index  for  manipulators.  It  is  the  ratio  of  the  volume  of  a 
manipulator's  workspace  to  the  cube  of  its  total  link  length. 
Based  on  this  index,  five  commercially  available  robots  can  be 
compared  and  rated  for  the  first  time  according  to  their 
workspace.  This  study,  it  is  believed,  provides  one  of  the  very 
basic  tools  necessary  for  the  evaluation  of  manipulator 
performance. 
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ABSTRACT 

This  papr  presents  an  analytical  and  computer-aided  pro¬ 
cedures  on  the  design  of  manipulators  Tor  optimum  workspace 
The  subject  is  treated  in  two  parts  In  part  1.  a  performance  index 
ts  derived  for  evaluating  the  workspace  and  determining  the 
dimensions  of  manipulators  The  effects  ot  geometry  on  the  K.ne- 
maic  characteristics  of  manipulators  and  the  design  parameters 
are  investigated  Based  on  this  analysis  an  algonU.ni  for  the 
design  of  manipulator  workspace  is  described  in  the  second  part 
The  effects  of  the  design  parameters  on  manipulator  workspace 
are  then  investigated  and  design  considerations  arc  presented 
lising  this  procedure,  the  "computer-designed  model"  is  demon¬ 
strated,  in  examples,  to  be  capable  of  obtaining  better  workspace 
than  the  conventional  models  available  commercially  at  the 
present  lime  In  addition,  an  optimized  model  of  the  UKRd  robot  is 
found  and  compared  with  the  proportions  of  a  human  arm 

I.  Introduction 

The  subject  of  optimum  design  of  manipulators  is  an  area  of  practical 
interest  in  which  little  has  been  done  and  reported  This  is  primarily  due  to  the 
difficulty  in  establishing  performance  criteria  for  robots  and  manipulators  One 
way  to  evaluate  performance  Is  from  the  kmr malic  viewpoint  Tor  instance, 
given  the  kinematic  structure  of  a  manipulator,  what  is  its  workspace?  Mow  run 
we  design  the  proportions  of  the  structure  in  order  to  optimize  its  workspace’ 
There  have  been  a  number  of  investigations  on  manipulator  workspace  - 
5]  Recently  Yang  and  l.cc  :  1.2]  developed  an  analytical  representation  of  mani¬ 
pulator  workspace  and  presented  the  first  performance  index  on  a  mempulalor 
based  on  workspace  A  computer  program  called  KAV  iKlnomatic  Analysis  of 
Manipulators)  was  developed  which  implemented  the  theories  and  cntrr.a 
presented  in  ;  1 1  For  a  given  manipulator  structure.  KAV.  first  defines  analyti¬ 
cally  the  workspace  then  explores  the  existence  of  hole  and  void  outlines  the 
boundary,  and  finally  calculates  the  total  volume  of  the  workspace 

In  the  present  study,  a  new  approach  is  introduced  \n  effective  algorithm 
Is  developed  called  KAPOV  which  combines  the  Heuristic  Optimization  Technique 
(HOT)  with  the  KAV.  to  provide  a  systematic  procedure  to  optimize  the  manipu¬ 
lator  structure  based  on  workspace  Two  problems  are  investigated  The  first 
problem  deals  with  the  optimization  of  the  RRRa  manipulator  and  the  s-rrnil 
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deals  with  the  optimization  of  some  commercially  available  robots  A  compara¬ 
tive  stuay  between  the  human  arm  proportions  and  the  optimized  RRR3  struc¬ 
ture  is  also  investigated 


2.  Manipulator  with  Kcvolute  Joints 

An  analytical  investigation  of  manipulator  workspace  is  given  by  Yang  and 
U'e  The  following,  token  from  [ represents  a  description  of  manipulator 
geometry  and  a  summary  of  the  results  on  the  analytical  representation  of 
workspace 

A  manipulator  with  n  unlimited  revolute  joints  in  series  can  be  represented 
schematically,  as  shown  in  Ttg.l.  There  are  n  coordinate  frames  to  specify  the 
configuration  Tor  any  coordinate  frame,  say  k.  the  Z,  axis  is  always  the  joint 
axis  and  the  A,  axis  is  in  the  direction  of  the  common  normal  between  axes  2 , 
and  Z*M  The  Z\  axis  is  attached  to  a  fixed  frame.  The  last  link,  link  n.  is  associ¬ 
ated  with  the  hand  or  the  end  effector  of  the  manipulator  The  link  k  is  con¬ 
nected  to  joint  k  Three  parameters  are  needed  to  specify  the  geometrical  rela¬ 
tion  of  two  consecutive  unlimited  joints  (joint  being  a  capable  making  a  com¬ 
plete  rotation)  Referring  to  Tig  2a.  these  parameters  are  a«,  the  common  nor¬ 
mal  between  two  consecutive  axes  Z,  and  Z* , t;  a*,  the  twist  angle  of  these  two 
axes,  and  b».  the  axial  offset  of  Joint  k-»l  on  the  axis  Z lt|  For  limited  revolute 
joint  (joint  with  physical  rotational  limit),  an  additional  parameter  fit .  the  loca¬ 
tion  angle  of  link  k»l  with  respect  to  link  k.  is  needed  as  shown  in  Tig  2b  A 
position  angle  fit  is  defined  as  the  angle  between  the  ,Y*  axis  and  the  A*,,  axis 
while  both  the  link  k  and  the  link  k»l  being  in  the  mid-range  of  their  joint  rota¬ 
tion 

The  workspace  of  a  manipulator  is  defined  as  the  region  which  can  be 
reached  by  the  center  of  the  manipulator's  hand.  H.  and  &'*(//)  denotes  the 
workspace  generated  by  the  point  H,  holding  the  axis  k  as  fixed  while  all  the 
rcvolule  joints  k.  k»  l.  .  n  make  rotations 

A  general  recursive  formula  for  workspace  >»,(//)  can  therefore  be  formu¬ 
lated  as  follows 

x,  =  r,eostJ* 

y*=r,sim>»  0) 

**  =  *«' 


vJiKre  r,  =  Vx*^  +  y,'!  and 


f* 

• 

V 

i 

=  ^]..l 

y 

X 

* 

• 

,  I  , 

(2) 


where  Aj*,i  denotes  the  homogeneous  transformation  wtuch  provides  the 
georn  tribal  relationship  between  two  consecutive  coordinate  frames  k»l  and  k 
Depending  upon  whether  the  revolute  joint  is  unlimited  or  limited,  there  are 
different  expressions  for  A ]*,!  For  unlimited  revolute  joint,  it  is  given  as. 


A  ]* .  i 


10  0a* 

0  cosa*  sina*  6»sino, 
0  -sana*  cos  a*  6*c  osa* 
0  0  0 


(3) 


Tor  limited  revolute  joint  '  1  ].  we  have 

cosfl t  -sin  0t  0  a* 

cosa,sin4»  cosa*cos0,  *ina»  fc»«ma* 
-sm  at  unfit  ~m  a,  cosj9,  cos  a*  6*  cos  a* 
0  0  0  1 


(4) 


Having  an  analytical  representation,  the  boundary  and  the  volume  of  the 
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workspace  can  subsequently  be  determined  by  the  construction  of  the  cross- 
secton  on  W,{ll)  cut  by  the  plane  X,'/,  [2)  Mathematically  this  cri  ss-section 
ni y  be  expressed  us 


plane 

*k  r  *rM  -  /  '$») 

z„  =  i,*  =  .d.) 


V) 


Because  the  H»(//)  is  symmetric  with  respect  to  the  axis  Zt.  therefore,  Eq  (') 
contains  all  geomcrical  information  regarding  the  shape  of  the  workspace  By 
scanning  sequentially  the  joint  angles  ,i5*(1  in  Eq  (5)  according  to 

their  rotation  limits,  at  a  reasonably  small  interval,  most  of  points  of  »',{//)  can 
thus  be  obtained 


3.  A  Performance  Index  Based  on  Workspace  Volume 


Given  a  manipulator  we  are  now  capable  of  evaluating  its  workspace  volume 
quantitatively  Furthermore,  the  development  of  a  performance  index  based  on 
workspace  volume  is  presented  '2]  A  definite  relationship  between  the  volume 
of  workspace  and  the  total  length  of  a  manipulator  is  found  This  is  given  as 
for  a  given  manipulator,  (i  e  the  twist  angles,  rotation  limits,  location  angles 
and  the  proportions  of  link  length  are  given),  the  ratio  between  the  volume  of 
the  workspace  and  the  cube  of  the  total  length  is  a  constant. "  i  e  . 


V 


JT~  constant  -  17 


;fi) 


The  notation  VI  stands  for  the  Volume  Index  of  manipulator  workspace  winch 
gives  an  indication  of  the  effectiveness  of  link  length  on  the  creation  of  reach¬ 
able  workspace  The  index  can  be  normalized  by  dividing  its  possible  maximum 
value  [21  Let  the  total  length  of  a  manipulator  equal  L.  the  maximum 
workspace  that  this  manipulator  can  possibly  have  is  a  sphere  with  radius  1.  and 
centered  at  the  joint  Therefore  the  maximum  VI  a  manipulator  c..n  jir.---iblv 
have  Is 

if) 


Using  this  value  as  a  normalizing  factor  a  normalized  volume  ir.drx.  di  nolei  as 
NV1.  can  be  obtained,  and  its  value  is  between  0  to  : 


Five  commcricallv  available  manipulators,  namely  I’t  .VA  CCO  ,  t.nmtal.on 
Inc).  ra  (Cincinnati  V.ilacron).  Pana-robo  (Matsushita  i.TD).  AID  BCD  .Automatix 
Inc.),  and  V30  (Cybotech  Corp  )  are  evaluated  based  on  VI  and  W1  The  results 
are  given  in  the  follow uyg  table 


- 

<UJ 

T  l.sngth 

Ibtume 

VI 

VI 7 

pm  a  f>  oo 

6 

i>\  0 

-.61906 

:  39 

o  3i: 

CU  I* 

6 

:f>:  o 

2290231 

2  21 

0  532 

Pana  -hobo 

5 

59  04 

305937 

:  66 

0  397 

Mil  000 

5 

59  02 

267150 

:  io 

03;C 

ct  no 

f> 

9z.  z9 

: 299251 

:  54 

0  160 

Table  1  Pi*rf  ynnancs  rtfUunfion  n_f  lndtesfrini  robots 
based  on  Ibfume  fnrfej 


4.  Optimization  Algorithm 

In  the  previous  actions  a  performance  index  based  on  manipulator 
workspace.  NYI.  is  pi c"  .ted  This  index  provides  an  objective  function  m  the 
optimization  of  manipulator  design  In  this  section  a  computer  program  named 
Kinematic  Analysis  and  Design  of  Manipulators  (KAPOV)  is  developed  which 
searches  for  the  optimal  design  of  manipulators  based  on  the  volume  of 
workspace  ft  combines  the  Heuristic  optimization  Technique  (HOT)  of  fzie  and 
Frcudenstcin  [8j  and  the  manipulator  analysts  program  IKAV)  of  lam  and  Vang 
[2|  The  design  para/m-li-rs  for  the  optimization  are  those  basic  kinematica! 
properties  of  a  manipulator  which  includes  the  common  normals  twist  angles, 
offset  distances,  and  location  angles  among  joints  The  objective  function  of  the 
optimization  is  the  normalized  volume  index  (VM)  The  remainder  of  this  sccton 
gives  a  brief  introduction  to  the  Heuristic  Optimization  Technique 
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The  heuristic  algorithm.  HOT  program,  developed  especially  for  the  pur¬ 
pose  of  mechanism  optimization  differ  from  most  of  the  conventional  optimiza¬ 
tion  methods  in  many  ways  First,  it  is  a  non-numerical.  effective  discrete  optim¬ 
ization  technique,  so  designed  to  integrated  well  with  the  mathematical  struc¬ 
ture  of  a  general  mechanism  design  problem,  which  is  usually  characterized  hi¬ 
king  iterative  procedures  and  nonlmcarities  The  algorithm  is  built  wilh  the  spe¬ 
cial  capability  of  by-passing  the  usual  obstacles  of  nonlmoariti  <■;. countered  in 
mechanism  design  problems  The  method  is  therefore  not  an  exhaustive  scare n 
method,  rather  it  is  a  systematic,  intelligent  guessing  technique  and  is  proba¬ 
bilistic  in  nature  Secondly,  the  algorithm  is  very  efficient  ar.d  simple  to  use 
The  attainability  of  a  solution  does  not  depend  on  continuity  and 
differenliub.iity.  and  it  is  essentially  independent  of  the  starting  point  of  the 
search 

The  method  is  a  discrete  one  and  therefore  requires  the  discretization  of  a 
problem  In  a  mechanical  system,  this  would  naturally  refer  to  the  conversion  cf 
a  continuous  parameter  variation  into  a  discrete  variation,  that  is.  me  conver¬ 
sion  of  the  solution  space  into  a  finite  number  of  states  of  combinations  called 
the  selection  space  The  parameters  are  allowed  only  on  certain  prescribed 
discrete  values  or  states,  chosen  in  accordance  with  what  seem  to  be  reasonable 
ranges  or  proportions  This  is  advantageous,  since  the  design  engineer  generally 
knows  reasonable  upper  and  lower  bounds  on  key  parameters  The  selection 
space  can  e  represented  mathematically  by  a  matrix  ;NxL)  m  which  there  are  \ 
rows  and  1.  columns  corresponding  to  S'  discrete  states  of  each  of  the  1,  design 
parameters 

To  apply  the  H  O  T  algorithm,  the  designer  needs  to  supply  only  I  »i>  subrou¬ 
tines  one  defines  the  feasibility  criteria,  the  other  specifies  the  objective  func¬ 
tion  Once  the  designer  selects  S  and  provides  the  lower  and  upper  bounds  of 
each  pamaters.  the  solution  space  is  dpflncd  by  the  input  lection-space 
matrix.  (A xL)  The  algorithm  then  generates  as  many  locally  optimum  solutions 
as  computation  time  permits  The  best  of  the  locally  optimum  solutions  may  be 
chosen  as  a  solution  to  the  problem  The  underlying  principles  of  the  combina¬ 
torial  heuristics  are  given  in  '6.7) 

In  the  following  sections  KADOV.  is  employed  to  optimize  two  special  mani¬ 
pulator  structures  and  some  cnmmercially  available  robots  to  deua  nslrate  its 
capability  and  potential 


b  Optimum  Design  of  Manipulators  with  KlUCi  Structure 

Two  manipulator  structures  which  are  commonly  adopted  in  industry  are 
optimized  through  the  KAPOV  program  The  result  is  then  compared  with  an 
analytical  solution  It  is  clear  that  usually  for  a  limited  revolule  joint  the 
greater  the  rotation  range  the  larger  the  corresponding  workspace  It  is  there¬ 
fore  reasonable  to  consider  all  revolute  joints  be  unlimited  in  a  general  study 

The  first  structure,  as  shown  in  Fig  3.  is  a  manipulator  with  six  revolute 
joints  and  with  design  parameters  a,  =  0.  a,  =  0  a,  =  0.  b ,  =  0  6- =  0.  6,  =  0 
6,  =  0.  6,  -  0  6,  =■  0.  n,  =  90° .  ne  =  0* ,  a,  =  0* .  a4  =  90* .  o,  =  90*  "d,  =  0* .  and 
a(.  Oj  and  a3  ore  dimensions  not  equal  to  zero  If  the- axes  of  the  last  three 
joints  intersect  at  one  point,  sav  point  H,  and  two  pairs  of  them  arc  orthogonal, 
then  these  three  revolute  joints  are  equvalent  to  a  spherical  joint  at  point  H 
Therefore,  this  manipulator  w,!l  have  a  RRRS  structure,  as  shown  m  Fig  4  The 
advantage  of  this  structure  is  that  its  entire  workspace  is  primary  workspace 
3.‘jj  The  primary  workspace  is  that  every  point  in  this  workspace  is  rrachable 
by  the  manipulator  s  hand  from  all  directions  Contrarily.  a  workspace  contains 
points  which  can  not  be  reached  by  the  manipulator  hand  from  all  directions  is 
named  as  secondary  workspace 

Referring  to  Fig  4.  let  design  parameters  0|.  and  a3  be  variables,  and  all 
other  part  of  the  structure  be  fixed,  the  task  is  then  to  find  the  optimal  propor¬ 
tions  between  as.  at  and  aa  in  order  to  obtain  maximum  workspace  volume  The 
result  from  KADOtl  shows  that  the  optimal  proportion  between  o,.  a?  and  oj  is 
Oil,  and  the  manipulator  corresponding  to  this  proportion  is  given  in  Fig  b 
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The  second  structure  studied  in  this  section  is  shown  in  he  6  The  design 
parameters  of  this  manipulator  are  a t  =  0.  Oj  =  0.  a,  =  0  a5  =  0  o*  -  0.  bj  -  3 
b,  ~  0.  fcs  =  0.  b0  =  0.  a,  =  90“ .  oj  =  O’  u,  =  90“ .  a,  =  90“  n,  =  90“  a,  =  O' .  and 
aj.  b,.  and  b-j  arc  dimensions  not  equal  to  zero  This  is  also  a  KH1I.S  strut  lure 
The  parameters  for  optimization  are  now  b,.  a2  and  b3  The  result  frem  KVK/V 
shows  the  optimal  proportion  among  b,.  a2  and  63  is  also  0  I  I  The  res.ut  : 
exactly  the  same  as  the  previous  example 

The  question  then  raised  is  whether  the  structure  obtained  from  KADuV  is 
really  a  global  optimal  or  just  a  local  extreme  In  order  lo  answer  it  the  follow¬ 
ing  concept  is  introduced  Referring  to  Fig  7  The  workspace  created  by  tins 
manipulator  equals  a  sphere  with  radius  1.  and  centrred  at  joint  Hecau«e  the 
hand  of  a  manipulator  can  reach  any  point  beyond  its  total  link  length,  there¬ 
fore.  theoretically,  the  possible  maximum  workspace  will  bo  the  sphere  with 
radius  the  total  link  icnglh  and  centered  at  the  first  joint  joint  This  proves 
that  a  RKR3  manipulator  with  a,  =  0.  a2  =  L/2.  a?  =  L/2.  b,  3  C.  t?  =  0.  bi  =  C, 
and  at  =  90“.  a2  =  0“,  aj  =  0*.  is  the  global  optimum  based  on  the  volume  of 
workspace  It  is  worth  noting  that  the  structure  will  remain  the  same  when  the 
design  parameters  a3  =  0.  b,  =  L/2  and  a3  =  90“ 

The  optimal  design  of  the  RRRS  manipulator  discussed  ir.  this  secticn  is 
kinematically  equivalent  to  a  THS  structure  ;here  T  denotes  a  hccke  joint) 
Since  Q(  =  b,  =  0  and  aj  =  90“ ,  the  first  two  revolute  joints  can  therefore  be 
replaced  by  a  hookc  joint 


6.  Comparative  Study  of  the  Proportions  of  Human  Arm*  with  the  Optimum 
KKKS  Structure 

Perhaps  it  is  interesting  and  noteworthy  to  consider  h,r  a  moment  ti  c 
structure  of  human  arms  Is  the  human  arm  an  "optimal  desugr.  1  hs  there  any 
similarity  between  a  human  arm  and  the  optimal  KHKS  structure  d.-cus-i  r! 
previously’’  There  are  apparent  differences  between  a  human  arm  and  !hc 
optimal  THS  structure  for  instance,  usually  anatomical  joints  ,.f  human  arm 
have  six  degrees  of  freedom  [Qj.  they  are  '"designed  1  for  two  arms  cooperation, 
and  ail  anatomical  joints  are  limited,  etr  However  kinematically  the  dm.-rv 
sions  are  surprisingly  close 

The  dimensional  comparison  brtwren  an  average  human  area  <ij  and  trie 
optimal  THS  manipulator  is  shownin  Fig  0  The  ratio  between  the  upper  «ji — 1  and 
the  sorcorm  is  ’.  0  02  However,  if  the  size  9of  the  hand  is  included  m  Ihc  r-n- 
sideration  and  the  grasp  position  of  the  hand  is  estimated  as  fiCZ  “  cosin’ 
(referring  to  Fig  0  where  0  27  is  the  average  length  proportion  cf  the  back  of  a 
hand  and  30*  refers  to  the  average  grasp  angle)  then  the  proportion  hctwi  or. 
the  upper  arm  and  the  distance  between  elbow  and  grasp  point  becomes  On 
TTiis  ratio  is  very  close  to  the  dimensions  of  the  optimal  TKS  structure  wh.ch 
equals  l  l.  This  finding  is  interesting  It  suggests  a  surther  analytical  investiga¬ 
tion  of  the  subject  of  the  workspace  problem  of  human  arms 


7.  A  Dimensional  Study  of  Some  Commercially  Available  Robots 

In  this  section  numerical  rxpenmrnls  of  the  link  lengths  of  fim  commer¬ 
cially  available  robots  discused  in  Chapter  V  are  performed  using  the  K-VDf'V 
program  The  purpose  is  to  search  for  the  best  proportions  of  the  link  lengths 
for  the  existion  robots  To  avoid  major  redesign,  all  non-zero  link  lengths  are 
allowed  too  vary  only  within  t  20*  of  their  original  length  and  ether  design 
parameters,  such  as  twist  angles,  location  angles  and  zero  length  links  remain 
unchanged 

The  results  of  the  optima!  proportions  found  by  KADOV  arc  shewn  m  Tub  e 
2  From  a  comparison  of  the  NTT  values  between  the  original  and  r»  »  designs  •> 
is  concluded  that  for  a  small  variation  of  link  length,  the  workspace  volume  for 
the  robots  CV  T 3.  I'ana-Robo  and  AID  000,  do  not  shew  significant  increase 
However,  for  PL’VA  600  and  CT  Y30  robots,  relatively  large  improvement  can  be 
achieved  Particularly  in  the  Pl'VA  600  robot  .  a  betler  ratio  o'  1  chime  to  link 
length  is  ohtamed  by  reducing  both  design  parameters  b,  ami  b...  whirr  b,  amt 
bj  are  offset  distances  of  link  and  2  respectively  The  parameter  6,  ;6,  =  1C 
in),  which  is  commonly  referred  lo  as  the  shoulder  of  the  Pl’VA  arm.  actually 
does  not  contribute  much  to  the  volume  of  the  workspace  Hv  eliminating  the 
shoulder  the  total  link  length  in  this  case  would  be  about  20“  less  than  the  ori- 
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gmal  one  However,  (he  workspace  volume  would  be  only  OS'!  less  and  the  NY1 
value  wouldbe  nearly  O'.Z  better  than  the  original  one  (Table  3)  There  may  be 
reasons  to  ineludr  the  shoulder  m  the  I’L’V.A  structure  However,  from  a 
kinematic  point  of  view,  eliminating  the  shoulder  would  optimize  the  workspace 


Robot 

Volume 

Link  Length 

NAT 

(in  *'3) 

(in  ) 

l’UV.A  600 

•.  H39U6 

t>:  0 

0  331 

I’L'MA  600 

102349 

a:  o 

0  632 

(No  shoulder) 

C  V  TO 

229S2A0 

-.0’.  0 

0  632 

Table  .1  .1  comparison  of  Pi’.UA  600  Ho  hot 
with  and  without  shoulder 


6.  Conclusion 

An  approach  for  the  optimum  design  of  manipulators  based  on  workspace  is 
presented  and  a  computer  code  KADO.V  is  developed  Some  practical  examples 
are  giben  and  an  interesting  comparative  stuky  with  the  human  arm  proportions 
is  provided  The  study  deals  with  only  one  of  many  objectives  for  the  optimiza¬ 
tion  of  manipulator  design  Some  other  important  engineering  aspects  such  as 
the  accessibility  (dexterity)  of  the  workspace,  dynamics  and  controls,  etc  need 
to  Le  considered  also  It  is  hope,  however,  this  study  provides  a  basic  approach 
and  an  algorithm  for  further  research  in  this  area 
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Table 


Improvement  nf  workspace  volume  1  a-cd  on 
in  the  link  1»  i.qth  proportions 


•  original  dc*»i<jn 
*•  KADOH'n  result :« 


Poirot 

Urii 

iqn 

Paramrti-i 

Lenoth  Proport  ions 

('  f 

0.56:1:0. 22:0.  ''4:0.11 

.331 

• 

i  PUMA  COO 

*V 

•>2 

h  ,  b 

1  r% 

•  3  5 

0.37:1:0. 05:0. Ql:0. 14 

•  • 

! 

1  1 

1:1:0. 2:0. 33 

0.532 

• 

C.M.  T 

&20 

*3 

f\  . .  h. 

4  > 

0.02:1:0.20:0.4 

0.545 

*  * 

2.4* 

0.75:1:0.15 

0.  3  *7 

. 

• 

Pana-KoliO 

V 

J3 

b4 

5.5% 

0.67:1:0.19 

.41) 

•  • 

0.75:1:0.13 

0.  U  .l 

* 

AID  BOO 

a 2 • 

V 

b4 

0.58:1:0. 14 

.■>.33’ 

•  • 

>.A% 

0. 8:1:0. 6 

0. 3f>R 

* 

CT  V30 

V 

V 

b4 

0.51:1:0.67 

0.44  3 

*  * 

z  c  % 
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Feasibility  Study  of  a  Platform 
Type  of  Robotic  Manipulators  from 
a  Kinematic  Viewpoint1 * III 
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1  Introduction 

I  radii  tonally .  lobot  aims  arc  anllnopomoi pine  open  chain 
mechanisms.  I  Ins  type  of  maitipulalor  usually  has  longer 
reach,  larger  workspace,  and  more  dcMious  maneuverability 
then  closed  knicmulic  chain  mechanisms.  Ilowcscr.  there  arc 
also  disadvantages.  The  cantilever-type  <>l  sirncimc  loi  the 
open  manipulators  inherently  has  poor  stillness  and  iheielore 
has  undesirable  dynamic  characteristics,  especially  at  Inch 
speed  and  Inch  loading  operating  conditions.  In  addition, 
will)  (he  exception  of  certain  configurations  |l|,  lire  piohlcm 
ol  solving  the  inverse  tnnetion  in  manipulator  control  is 
always  difficult. 

Perhaps  there  is  a  possible  alternative  design  tor 
manipulators  using  the  closed  kinematic  chain  su.  Ii  as  ihc 
Sleward  plalloiui  |2|  and  olhcis  |.1J.  I  Ins  would  be  p.u 
licnl.ulv  Icasiblc  in  applications  where  dviiami.  loading  is 
seven'  and  vei  the  demand  oil  vvoi  kspacc  and  ui.mcuvi  i.U’thlv 

is  low.  I  or  nisi, in. e.  il  could  he  used  lor  dcvclopimii . in 

"intelligent”  multi  deerecot  ■fieedom  posuionei  w..ich 
would  be  versatile  in  many  uses  such  as  Inch  piccision 
machining.  welding.  painting.  and  mans  military  ap 
plications. 

There  have  not  been  many  investigations  o!  manipulators 
will)  closed  kinematic  chain  on  l tie  resold.  I  he  Steward 
pl.Uloim  |’|  winch  is  consirucicd  by  connecting  two  plates  to 
sis  a<l|iisi.ihlc  lees  and  is  a  six  deercc-ot  lieedom  6  SI’S 
pl.illoiin  inch.niiMi.  (S  and  I’  denote  the  spherical  and 
prismatic  loinls.  respectively),  was  originally  designcvl  as  an 
aitaall  sumilaioi .  and  was  also  suggested  lor  the  application' 
of  machine  loot,  space  vehicle  simulaioi.  iraustu  machine, 
etc  I  atet,  Holtm.in  and  McKinnon  |4|  tried  to  smuilaie  ihc 

1  H.t  .nl  Ml  f.nl  MM  ill,  .ll.s.lMlHMI  Ml  III.  Ill'l  .IMlhMI  Ml  p.Mti.il  In  I  f  tllilu'li:  .'I 
tin  M'v  |  (I  tl  yl  |  H '  1 1 1  li«t  lilt  Ihulill  l»t  III  ll»-  Sxfioul  o|  \  lil'tlli  1  l  till' . 
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ail  crall  motion  via  tins  type  plat  lot  in  t’v  appiv  me  an  s  \  P  I  \ 
finite  element  progtam.  Met  allian  and  I'uoug  |5|  used  this 
device  as  an  auiom.uic  assembly  table.  I  hev  dc'eiibed  1?|  die 
mechanical  hardvv.uc  and  the  probiem  ol  pad;  'vntlie'is  "a- 
als*'  inv C'lieaicd  1 1  ;tiis  book .  II  ml  | 'I  adopted  ’.be  V ew at d 
plaltoim  as  a  mechanism  lot  a  iv'bonc  ,um  I  olhnvm.  'h.u 
idea.  I  ichtci  and  MeDv'well  |('l  presented  a  review  and  some 
pi chmm.u v  design  cotuvpis  on  iltt'  type  ol  mampula'.o: . 
Recently.  Mum  |”|  iinJcnook  a  osienianc  'tuvlv  oi>  die 
paiallel  aelualeil  lobv'lic  arm.  in  which  manv  [vossii'le  ap 
plicablc  ut  paiallcl  'Itucinrcs  «eic  tevicwcd  and  vvtti.ii  u 
chided  the  Steward  plaitoi  tit. 

Despite  its  potential  usetulness.  i Mis  class  ot  closed 
kmemaltc  chain  niantpnlalv't '  vtoes  have  pltvst.al  luuuatioio 
which  it  is  important  lo  uvoeiu/c  so  that  kinem.mc  svitthe'is 
methods  in. iv  he  developed  vvuli  these  .onsiiaiuls  and 
chaiacteiisiies  in  iiiuul.  Ii  is.  theieloie.  ne.essaiy  to  earn 
deeper  uiulei standing  ol  ihi'  ivpe  ol  mechani'iu  and  to  obt.mt 
some  design  and  apphealion  guidelines.  (  otiseviucntly .  a 
leastbiliiy  study  ol  tho  mechaiusiu  tit  (lie  lielil  v>l  tl '  poiential 
application  as  a  ivtbvme  positionct  is  hvilh  impotram  and  of 
practical  interest 

Ihc  purpose  ol  this  investigation  is  to  develop  analytical 
methods  and  computer-aided  pi occdures  capable  ol  analv  /mg 
the  basic  kinematic  eharaetet  t'lics  v>l  the  pl.nlo-tti 
mechanism,  such  as  the  exl rente  ranges  ol  motion  as  well  as 
workspace.  In  pailieul.tr,  the  objectives  ate:  II  to  petlonn  a 
number  synthesis  on  lire  esiMeiisC  ol  a  single  (or  manv  I 
platform  sirueliuetsl  which  is  suitable  tor  luneliomng  as  a 
robotic  positionct;  ')  to  fointulate  an  aualviic.il  model  lot 
litis  class  of  mechanism  and  to  provide  a  displacement 
analysis  including  the  solvability  ot  its  mvetse  conttol 
dilutions;  1)  lo  analytically  lormulale  the  phvsic.il  .on 
Miami'  of  llns  mechanism,  such  as  the  ioi.it. Utility  ol  legs  and 
ball-and-socket  toiiiis;  and  41  ultimately  to  investigate  the 
wot  kspacc  and  inuttcm  ci  ability  ol  the  mcch.iiusni. 

Although  pievtom.  mv  e  a  teal  u»n  vat  tin  me.  h.uu-iu  It.iv, 
uichulcd  ceil.un  aspects  ol  I  he  ptoblcms  niciiiioucd  Itetc.  sit.  It 
,rs  mobihtv  [2.  ft,  7|  amt  thspla. eruertl  felaiion.lup  |2.  V  r.| .  it 
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is  believed  that  a  comprehensive  investigation  ol  these  suh- 
jecis,  especiallv  from  the  consideration  of  physical  constraints 
and  workspace  using  analytical  and  computet  anled  means,  is 
much  needed  and  the  effort  would  cei  lainlv  he  iiotewoilhv 

2  Number  Synthesis 

It  is  well  known  that  a  mechanism  which  is  used  as  a  general 
robotic  manipulator  should  have  sis  ot  more  controllable 
degrees- of-trecdoin  (d.o.f.).  lor  manipulators  with  ar¬ 
ticulated  configuration,  in  many  cases  ot  simple  applications, 
the  number  of  controllable  d  o  t.  can  be  less  than  sis  if  cost 
and  simplicity  are  a  consideration.  The  Steward  platform  |2) 
has  sis  extensible  links  and  has  sis  degrees-of-freedom.  I  here 
are  a  few  basic  questions  one  can  ask  regarding  the  platform 
mechanism  or  similar  mechanisms,  l  or  instance,  whs  is  it 
necessary  to  have  sis  connecting  legs,  or  equivalently,  what  is 
the  relationship  betw  en  the  number  of  legs  and  the  number 
of  controllable  d. o.l .  i  this  type  of  mechanism?  Is  it  possible 
to  introduce  more  legs  into  the  mechanism  and  yet  still 
maintain  sis  conltollable  d.o.f.?  What  is  the  possibility  of 
having  controllable  d.o.f.  with  more  than  or  less  than  sis? 
I  his  section  of  the  paper  deals  with  these  problems. 

2.1  n-SS  Platform  We  begin  with  a  discussion  of  the  n-NS 
platform  as  shown  in  Tig.  1.  Since  the  connecting  joints 
between  the  legs  and  the  plates  ol  (lie  plaltoim  ate  hall  joints 
and  the  d.o.f.  of  hall  joints  are  not  controllable,  we  iherefote 
first  determine  the  number  “n"  so  that  the  mechanism  would 
has  e  zero  degrecs-ol -freedom. 

A  general  form  of  the  degree-ol'-freedom  equation  tor  both 
planar  and  spatial  mechanisms  can  he  vv  i  it  ten  as  lollows  IS): 

/•'-  Ml-J-  II  ►  /.,  (11 

where 

/  the  effective  degree-ol'-freedom  ot  the  assembly  ot 
mechanism, 

X  the  degree  ot  freedom  ol  the  space  in  which  the 
mechanism  operates  (lor  spatial  motion  X  (•.  and  for 
plane  motion  and  motion  on  a  surface  X  A  t . 


Flq  I  A  fi  SS  platform  rvmrhmtUfn 


/  -  number  of  links, 
i  -  number  of  joints, 

./,  degreo-of- freedom  of  all  joint, 

/  idle  or  passive  dogtcc  ol  l  iccdom. 

It  is  known  that  genetally  there  are  special  gcomctiis.il 
coiulilisms  which  need  to  he  taken  mts>  consideration  m  the 
deter niinaiion  ol  tbs'  d.o.t.  on  spatial  mechanism..  spsctaJ 
cases  arc  olieii  associated  wuh  parallel,  mieisectmg.  ot 
per  pendicular  joint  avis.  (  ouscqucntly .  the  modi  I  v  mg  ts't  in  ! 
takes  eaie  ol  this  idle  d  o  t 

Referring  to  Tig.  I.  ws  have  x  n.  /  n  *  2  (iiishishng  n  legs 
and  2  plates).  /  2/t  and  t  >  ttoi  each  ball  |omt,  the  d  o  l  i- 

?).  It  is  also  noted  that  each  pair  ot  ball-and-socket  iomt- 
which  associate  wuh  the  same  leg  can  (inn  with  ts.pe.t  to  tits 
avis  ot  the  leg  freely:  tlieictoie.  tins  is  an  idle  d  o  I.  whuh 
dsvs's  not  com ti hute  to  the  mobility  ol  the  mechanism  u-sli 
C  onsequently ,  only  fiv  e  s>1  the  siv  d.o  I .  ot  the  ball  iomt  pan 
associate  with  the  mobility  ot  the  mcsh.mi'm  \k  e  have, 
therefore,/,  •-  «  which  corresponds  to  the  tmitiKt  otsxp.u: 
or  Ills'  number  ol  legs.  I  sing  equation  1 1 ). 

/  r>  n  i  h  i 

1  sjliation  (2l  show .  (hilt  w  hen  n  equals  so  .  the  do!  . . 

platform  osiuals  .vie.  i  e  .  the  t.-ss  plat f . o  t:r  i  a  ■. 

When  n  is  less  than  st\.  we  will  has  s-  <r>  \t  u:;.  t"  ;  .o- 

d.o.f.  atul  when  >:  is  etc. not  ill. in  sis,  the  . . 

constraitis'sl 

2.2  n-SI’s  on  n-s(  si  Mechanism.  Sm.s  1,  i-  " 
ms'chanism  t»  a  iicul  sttn.tmc.  u  is  not  ni.uun  i 

older  to  ptovislc  coiiliollcsl  nioiion  on  llns  nu.i  n.  ; 
cylindrical  (s>t  prismaiist  joint  is  added  to  cash  leg  Ksis- 
to  Tig.  2.  each  s> t  these  cvlmdiic.il  tonus  can  be  .oturo  ,  !  • 
a  tns'lor  or  linear  iistuaiot .  and  the  length  ol  c.n  h  K  e  ,.,;i  :•  ■ 
be  adjusisil  itisliv  isfti.illv  \  svliiivliical  iomt  ha-  t«  >  d> 
ol-liccdotn.  Ilows'vct.  using  a  'imilai  .iigumcn'  a  n  ",  . 

section,  onlv  the  iianslaiiona!  degree  ot  tis'Cvlom  alts, 
mobility  ol  the  mcsliati.-m  an  I  the  ioialion.il  'uvdot; 
redundant  one  I  Ins  MIiIlVM*  I  lk.lt  ilk  SC  S  i^'j  m  ;  ;i,j' 

ticulai  c:»nC  i **  kiiKMiuiikjllv  s  kjui k dlcui  u*  dii  sps  ^  V\  ,  i 
.Hu)  l>  represent  the  e\ limit u.il  .»nJ  priMiuiiK  p.nts.  ’sS” 
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Inch.  I  or  a  n  SI’S  pkiilnnn,  »o  li;nc  /  2/i  •  2  amt  /  I n, 
thcrcloic. 

/-6|<2n-2)  l«  " 


6/it6*  6 rt  /  /i  n  6  (  I) 

I  lie  result  liom  equation  (.1)  slums  dial  a  n  SI’S  pkillnim 
always  has  si\  d.o.f.  and  among  ilicin  n  d.o.f.  are  con¬ 
trollable.  It  n  is  less  than  six,  there  exists  (6  n)  un¬ 
controllable  d.o.l.  If  n  is  greater  than  six,  there  exists  (n  -  6) 
wasted  controls  which  are  needed  to  compensate  the  oxer- 
constraints  as  can  be  seen  in  equation  (2).  Only  when  n  equals 
six,  there  are  six  controllable  d.o.l".  Therefore  the  6-SPS 
plat  form  is  the  only  mechanism  of  its  type  w  hich  can  be  used 
as  a  general  robotic  manipulator. 

Another  interesting  question  could  be  asked.  Can  a  plat¬ 
form  with  n- SI’S*  (6  -  n)- SS  structure  sene  as  a  n  d.o.f. 
manipulator  for  n  less  than  6?  Theoretically  the  answer  is  yes, 
but  practically  it  is  no.  The  reason,  which  will  become  clear 
later .  is  that  m  specifying  the  leg  lengths,  the  problem  of 
wlcntitung  the  corresponding  position  (both  location  and 
orientation)  ot  the  top  plate  would  become  a  very  difficult 
task,  on  the  other  hand,  it  the  desired  position  of  the  top  plate 
is  prescribed,  it  is  rather  easy  to  obtain  the  corresponding  leg 
lengths  Howe'er,  this  is  under  the  condition  that  thcic  must 
be  'is  adiiist.ible  legs  We  may  therefore  conclude  that  it  is 
impi, isiic.il  lot  the  plait  or  m  type  ol  manipulators  with  d.o.f. 
less  itian  si\,  whereas,  in  contrast,  the  situation  often  is 
allowed  hi  i  tie  nincutaied  type  of  manipulators. 


.1  Displacement  Analysis  ami  Controllability 

h  w.r.  > on,  Inded  hi  the  previous  section  that  the  6  SI’S 
pi  ut. mu  ss Ins l>  Imp  sis  sonliollahlc  degtees  ot  liccdom  is 
poiciniall'  iiselul  as  a  gcncial  manipulating  devise.  Ills' 
purpose  ol  this  section  is  to  present  a  mathematical  analysis  of 
Mic  mpiii  output  relationship  and  subsequently  to  demon- 
sir,  ue  the  soieioll.tbilily  ol  the  mechanism. 

.VI  Displacement  Analysis,  lo  analytical!'  model  a  h  SI’S 
plat l or m  (I  ig.  2k  we  shall  at  the  beginning  define  the  con¬ 
trollable  elements  and  the  desired  outputs  ot  the  mechanism. 
I  irst.  let  ilic  sontrolling  of  the  position  of  the  center  point.  />. 
ol  die  lop  plate  be  the  desired  objective,  i.c..  the  point 
represents  the  hand  lot  lire  end  el  lector  I  of  lire  manipulator. 
Second,  lei  each  ol  the  legs  of  lire  plallorm  be  adjustable  and 
indisidualls  controllable,  i.c..  for  each  prismatic  ioim  (or 
c x linclric.il  lomll.  thcic  associates  with  it  a  sei'omotoi  (or  lire 
likd  llurd,  I oi  simplicilv  and  wiihom  losing  geneiahlv.  lei 
both  lire  lop  and  lire  honour  plates  be  citculai  in  shape  with 
the  top  one  assumed  lo  be  movable  and  the  bottom  one  fixed 
arid  has  mg  its  center  at  It  f  urthermore,  let  the  distribution  ol 
each  pair  ot  hall  joints  on  the  bottom  plate  be  symmetrical 
with  respect  lo  each  ol  the  three  radii  located  at  120  deg 
apart  from  one  another  on  the  plate.  1  he  ball  joints  on  the  top 
plate  have  similar  arrangements.  The  platform  mechanism  is 
then  an  octahedron.  Denoting  R  and  r  as  the  radii  ot  the 
bottom  and  the  lop  plates,  respectively :  I,  lliiougli  /„  as  the  leg 
lengths;  /„  as  lire  nominal  length  lot  all  legs;  and  ■'/  as  the 
distance  between  the  two  plates  when  all  leg  lengths  ate  equal 
lo  the  nominal  length 

I  o  give  an  analytical  representation  ol  l lie  kinemalics  ol  l lie 
h  SI’S  mechanism,  lire  lived  (  artesian  cooidinatc  frame  is 
selected  at  point  H  with  the  /  axis  pointing  vertically  upward 
and  the  X  asis  passing  through  point  II, .  I  he  locations  of  the 
ball  ioiiiis  on  the  lised  boilnm  plate  can  l hen  Ik-  established 

II,  IK.0.01' 
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where  O  represents  the  angle  between  ball  joints  H  and  H  . 

/(,.  and  H,.  and  H.  and  H„.  and  and  Y,  leptesem  ills  a" 
Oand  sin  O,  respectively. 

Similarly,  let  a  moving  ('artesian  coordinate  frame  1 7")  be 
associated  with  the  top  plate,  having  its  origin  at  the  point  />. 
its  r  axis  normal  to  the  plate,  and  the  v  axis  passing  through 
the  point  h,.  Therefore,  the  relative  locations  of  the  ball 
joints,  /),' s.  with  respect  to  the  mov  ing  frame  are 
/>,  -  (r.O.OJ ' 

b:  -  (/•(  „. rV... ()J ' 
b,  -“-  |rC,;„.r.S,;„.())' 
b*  '  |rr,:„....r.Si;„. 
b .  |r(  Ajo.r.S  -4...()| 

/),,  =-  \rC:il . r.S.4,  ..or  1 5 ) 

where  0  denotes  the  angle  between  ball  joint'  b.  and  b  .  b 
and  b 4 .  and  Ik  and  />,. . 

l  et  the  point  />„  and  the  Cattest.m  coordinate  I  /j.  be  the 
position  of/)  when  the  top  plate  is  at  its  original  (the  nominal) 
position  as  shown  in  fig  4  I  he  geometrical  relationship 
between  the  coordinate  frame  |  /I  and  the  fixed  coordinate 
frame  |/<|  can  be  represented  N  a  4  •  4  homogeneous 

transloi mation.  and  n  is  as  lollows; 
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w  here  >  -  t )  H  2 . 

Now  assume  the  center  ol  the  top  plate  moves  away  Horn  n- 
nomiual  position  />,  to  a  new  desired  position  /’  and  base  r he 
transformation  |71  with  respect  to  i he  fixed  coordinate 
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Where  the  (  v^.i,, .;,, )  are  the  ("artesian  locations  of  the  point  /’ 
.  and  (r/n  ,rfM  ,dw ).  Ir/,'. </■.(/,),  and  (</, ,  .i/m.i/,,)  are 
direction  cosines  of  the  axes  s ,  e,  and  with  respect  to  the 
fixed  coordinate  |/f|.  correspondingly .  I  herefoie.  the  locaium 
of  each  ball  joint  on  the  mining  top  plate  with  respect  to  the 
coordinate  frame  \U)  can  be  obtained  from  equations  (51  and 
(7).  and  tiles  are 
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J.2  Controllability.  The  controllability  of  a  manipulator 
could  be  qualitatively  assessed  by  the  degree  of  difficulty 
involved  in  the  solution  of  the  inverse  kinematic  problem 
which  is  the  determination  of  the  corresponding  cootdinales 
of  the  manipulator  by  prescribing  the  desired  position  of  the 
end-effector  in  the  lived  Cartesian  coordinate  linme.  I  or 
manipulators  with  articulated  configuration,  the 
soi responding  coordinates  are  the  rotational  angles  of  the 
resolute  joints.  |-oi  manipulators  with  platloim  con 
figuration,  they  are  the  corresponding  leg  lengths. 

l  et  the  desired  position  of  the  end-el tector  ot  a  platform 
manipulator  he  |  /]r.  The  absolute  locations  of  the  ball  toint 
b,' s  on  the  top  plate  can  be  derived  from  equation  (SI.  I  he 
required  leg  lengths  can  he  obtained  by  calculating  the 
distance  between  each  corresponding  pair  of  ball  joints,  h  \ 
and  It,  V  from  equations  (4)  and  (X),  we  have 

/;  (r/n r  r  t„  *)•’♦  (r/:lr +»•,,»-’*  (</Mr  •  (Vo) 

I]  r  (r/|,  rC„  *  r/, , rS„  4  v,.  -  tt( 

•  „  t  t  vr  K. S,,) 

•  (r/i,  r(  t  it , . rS„  t  '  ('>b) 

l\  <(/,  r(  ,  i  ft | . rS |  t  v,.  K(  ,,,)' 

•  I </  ,i<  | i  it  rS,  „  i  »■„  It S,  ,,l 

i  (r/,,r(  ,  i  r/,.r.S|.„  r  .-  ,|  ('It  ) 


=  (r/urC;,. . „  i  *  \r  «(  ;j. 

+  (r/:,/-C-4„  ♦  rf.«r.V;4ll . •  r,.  tt<  .  ,.) 

t  (f/j.rC.jH.  -  4  r/i.-r\;4,, ,  „  *  I'l/I 

I  quatioir  (V)  represents  a  set  ot  evplicit  input  output 
equations  of  the  h  SI’S  platform  mechanism.  I  or  a  pic'ctibed 
end-effeeior  posiiion.  the  coricsponding  leg  lengths  can  be 
readily  obtained;  however,  the  direct  kincmatis  control 
problem  is  a  rather  dilfienlt  task  since  it  involves  the  solution 
s>l  a  set  of  highly  nonlinear  simultaneous  equations.  I  he 
situation  is  found  to  be  exactly  contrary  to  the  case  of  a  r  obot 
with  open  kinematic  chain. 

4  Holahililv  of  Ball-and-Sockcl  Joints 

A  ball  joint,  theoretically,  can  move  lively  svrili  respect  to 
all  three  Cartesian  axes.  In  practice,  however,  this  is  not  tine. 
The  motion  of  the  ball  joint  is  always  restricted  because  ot  us 
physical  dimensions.  A  ball  joint  usually  includes  three  parts; 
the  ball  head,  socket,  and  the  connecting  leg.  Referring  to  lig 
T.  let  the  radii  of  the  hall  head  and  the  connecting  leg  be  It 
and  it/ 2,  respectively.  To  physically  hold  the  ball  head  in  the 
socket,  the  holding  width  c.  as  shown  in  the  tigurc.  must  be 
grealci  than  zero  I  el  ■<  denote  the  rotation  angle  of  the  toint 
oil  the  V/  plane,  the  iot.ilroii.il  Imnis  ot  ihc  angle  ,.iu  be 
determined  as  follows 
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Fig.  3  The  principal  cross-section  ol  a  ball-and-socket  joint 

z 


Fig.  4  A  6-SPS  plattorm  mechanism  with  bail  joints  evenly  distributed 
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(11) 
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I  or  Ira  rotation,  »>  should  be  It'"  Ilian  or  equal  lo  ir  2  l; 
Irom  cqu.il ion-- 1 10)  ami  (II),  we  have 

it 


C  -) 


i  sin 


N,. 


"?s"'  \  (U> 
the  complete  rotational  range  ol  ■  <  on  the  V/  plane  is. 
Iherelore. 


(  ■  ll  ,  •'  \ 
t  I  III  I  Mil  I  "  • 

'  IK,  K,.  > 

I  (illation  (P>  is  nseltil.  ll  tepresenls  the  physical  eonsirainis 


sill  1  i  Mil  1  (H) 

2K.,  K, 


Fig.  5  A  planar  four  bar  mechanism  with  two  adjustable  finks 

lit  ball  joints  auj  prowJcs  ihc  pt  .tci  k\tl  Josiytt  tor 

these  joint''. 

5  W  orkspace  anti  Muni-uxcrahilitx 

I'lie  problem  ol  anals/ing  the  shape  ol  the  vtoikspace  ol  -i 
6-SI’S  plat  lot  in  manipulator  i'  a  rather  difticulf  one  lr 
addition  to  the  inherent  complexify  in  dealing  ssilh  'patia: 
geometry,  the  mechanism  i'  chaiacicn/cd  In  Ms  tnultiloop 
Mi  net  lire  and  its  displacement  behas  lor  sy  Imh  ms  ol\e->  a  'd  o' 
highly  nonlineai  simultaneous  equation'.  lor  'implkiis. 
consider .  iherelore.  a  special  case  in  almli  all  hall  loitst'  are 
ex  only  distributed  along  the  pcuphcucs  ol  both  it1e  Kmon 
and  the  top  plates  (i.e..  both  angles  i)  and  "  arc  ol  60  dec!  u- 
shown  in  I  ig.  4.  lire  I'.imc  approach  employed  i'  lo  in 
vestigaie  die  woikspace  ol  din  mechanism  ttom  n-  co" 
section  on  th:ee  particular  planes.  Kcfcning  lo  I  ig.  4.  these 
are  the  A/  plane  uhich  contain  the  ball  joint'  /</>,/'./<.  and 
the  planes  found!  by  mining  die  \/  plane  with  i C'pc, I  io  die 
/  asis  with  -in  angle  ol  60  and  l^d  deg.  re'Pcclixels.  lo 
further  'imphty  the  ptoblent  sse  coii'ider  a  .pecial  cum 
yy herein  the  motion  ot  the  top  plate  ot  the  inechani'in  i 
allosyed  only  to  rotate  sxuli  tespeet  lo  the  )  axis,  t.c  .  die 

lol.n.il'ihl y  w  i ill  I 'i >di  \  ami  /  ase'  is  i es) i  icted  I  fie  pi  o I’ lei n 
llieieloie.  i'  tedmej  lo  limhue  die  le.uh.iPle  -pu.c  ol  die 
ceniet  point,  /<.  on  the  couplet  ot  a  plaiuu  4-bar  meeli.iiiiMi. 
tas  'Mown  in  t  ig.  5)  ssuh  adjustable  link  lengths.  /  ami  /.. 
subiecled  lo  s annus  eonMiaints  Mich  as  die  extieme  linin'  o! 
all  legs  and  the  lotatabilus  ol  all  joint'.  In  the  folloxxtug  die 
displacement  equation  ot  t  Ills  planat  adjustable  4-ba; 
mechanism  is  lirsi  dctiseJ.  I  lie  inlet lerence  condition'  ate 
i hen  considered. 

5.1  Displacement  Analysis  of  A  IManar  4-ll-ar  With  Ad¬ 
justable  I  .inks,  figure  5  shows  a  planar  4  bar  mechanism  with 
adjustable  link  lengths,  /,  and  /4.  A  coordinate  frame.  |/<|.  i' 
located  at  the  point  It  with  /  axis  pointing  sonic-alls  upssaul 
and  the  A  axis  pointing  horizontally  led  II,.  b.,  h,  and  H, 
represent  font  ball  joint'.  I  he  link  lengths  ol  the  coupler  ami 
the  fixed  link  are  die  diameters  ol  the  top  and  -tie  bottom 
plates,  i.e.,  2/  and  2K.  respectixels.  I  he  coordinates  ol  tin- 
fixed  hall  joints,  /),  and  /<,.  on  the  bottom  plate  are  I  A’,0)  ami 
(  K,0).  res  peed  \  els ,  and  coordinate  ol  die  moving  ball  lom: 

t) ,  on  the  top  plate  is  [K  /,cos<i| ,  /,  sin<>  ).  Speeds  iag  s  aim-' 
ot  the  input  parameters  u,.  ,  ami  /4.  then  the  output'  which 

are  the  corresponding  coordinates  of  points  h4  an. I  />.  and  the 
cortcxpoiiding  angles  of  o  and  >■,  can  be  louml  It  i'  out 
propose  lo  dense  die  closed  toiin  mpm-outpul  displ.ueimm 
relationship. 

fquating  the  distance  betsseen  the  two  ball  joints,  b  and 
/>4,  to  the  length  ot  the  coupler.  2r.  sse  base 

'll  Is  (M’llll  IIKlllli'  lull-  ih.u  thru*  f'|xl‘  .1  -py.  '.ll  sinftil.inl  V  >  Olid '  "  *' '  '  >' 
this  OCItlv  ll|\U  |b||K'd  y.lM’.  I'M  wll'cll  I  hr  1 1  iVI  flN  i  •  t  'lu.  It'll)  It. .11 
pfCM’llU'd  in  SOillOII  '  |X  tK»l  .ipi'ltc.iMi*  (i  I'-  tht.  pt'-H'.'U  th.|f  ,|lt  k'f  i  - --  .r 
riju.il  >u ill  bi»l It  the  hollom  .m  l  die  top  pl.itcv  .nr  p.i'dlil  I  his  ..ondtiton  - 
Kfiff.tr  fit  tiff  Hi  If  l  now  ff  *  I  »  *  ll  I  hi  l  1  •  if  p.ti  till-!*'-  *  •!('  "fill  fir  i.I.f" 

hill  nl  i  i|it.lf  y  f  .ink  l«  fi»*l If  '  •n,n  .  i  »tir  • 1 . . t*  1 1  pi  •  i»*  - ''  'hi  ■  »Kpf<  1  Hi!  t  'i 

It. tint  In  ihn  c.isr  ot  sjmimI  ptnini •  urn  o*  ImV-  Mu  rit.i.tl  *U  u.  . 
tun  loin  i*i|ii.in<*ii  n- n. »lk  iUs  tm'  ipp1\ 
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Fig.  7  A  method  to  determine  the  distance  B2D2 


Fig.  6  Geometrical  parameters  ot  the  platform  mechanism 


646,  '■=  (*4  -/?  +  /,  C„,  )2  +  (Z4  -  /,  S„, ):  =  4r:  (14) 

Equating  the  distance  between  the  ball  joints  64  and  B 4  to 
link  length /4,  we  obtain 

(at4  +  /?)*  +  (z4)“  —  /4  (15) 

Subtracting  equation  (14)  from  equation  (15),  we  have 
/tWi-^-2R,C"1_2  «~/,C 
2/(5„,  /,.V„, 

Let 

..  /?  +  /4  ~4r2  -  2/?/i  C„, 


2/? -/,<:,„ 


then  we  have 

Zt-P-  Qxt  (19) 

the  coordinates  of  point  64(je4,;4)can  be  solved  explicitly  and 

-(/?  -  PQ)  ±  vT*  -  /><?)-  -  ( 1  +  Q- ){P2  +  R :  -  /; )  ,,A4 

'■ - TiTQT, - 1201 

Subsequently,  the  coordinate  of  the  point  p(.xn,z,,)  is  given  by 
x,  +  K~t,C,n  } 


Extensibility  of  l.egs.  For  a  6-SPS  mechanism,  each  leg  is 
limited  by  its  extreme  extensible  length,  i.e..  llu, . 

Assuming  all  legs  have  the  same  limits  of  length,  the  motion 
of  the  point  p  (equation  (21))  is  then  constrained  by  the  ex¬ 
treme  lengths  of  the  legs,  6,  B,  b:B:  ,hxB\  and  64W4,  i.e., 

— /|  >(;./>. (4  (24) 

Since  the  top  plate  is  symmetric  with  respect  to  the  X/.  plane, 
the  legs  6SB<  and  h^Bk  are  mirror  images  of  the  legs  6,fl,  and 
h;B:,  respectively.  Therefore  they  are  not  included  in  this 
discussion. 

For  a  given  leg  length  /,  and  /4,  it  is  necessary  to  know 
whether  the  corresponding  leg  lengths  /.  and  />  are  within  the 
extensible  limits.  Referring  to  Fig.  6  we  have, 

/;  =  B,a, '  r  b,a, '  (25) 

where  the  a  points  represent  the  projections  of  the  b  points  on 
the  XY  plane.  The  distances  //,</,-'  and  B, are  given  as 
follows: 

=  (K-2  rC0-.v4): 

=  (0.51?  -  l.5rC0  -  v4):  +  {RSH, -rSH,): 

BxOx 2  =  (  -  0.5  -  0.5/-C,  -  .v4 );  +  (RSM  -  rS„,  F 

B,o4:  =  (  —  /?  —  ,v4 ) '  (26) 

and 

6,0,  =  C4  +  2rS,, 

b:a:  =  c4  +  l.5r.S'., 

6,o,  =  z4  +  0.5r.S'o 


--  2 

The  corresponding  angle  n4  can  be  obtained  from 

«4=tan  1  (22) 

.v4  +  B 

The  corresponding  tilt  angle  of  the  top  plate,  <j>,  is  given  by 

*  =  ,an'^fe^  (2?) 

Equations  (l9)-(23)  represent  a  set  of  explicit  equations 
defining  the  displacement  of  a  planar  4-bar  with  adjustable 
link  lengths/,  and/4. 


where  the  tilt  angle,  0,  of  the  top  plate  is  gisen  by  equation 
(23).  The  desired  corresponding  leg  lengths,  /.  and  /,,  can  then 
be  calculated  from  equations  (25)-(27). 

5,3  Rotatability  of  the  Ball  Joint  on  the  Bottom  Plate.  I  he 

workspace  of  the  center  point,  p,  of  the  top  plate  (equation 
(21))  is  the  space  bounded  by  all  the  constraints  attributed  to 
the  ball  joints,  as  given  by  equation  (13),  in  addition  to  the 
constraints  due  to  the  interferences  of  links.  The 
corresponding  angles  of  the  ball  joints  on  the  bottom  plate 
can  easily  be  derived  in  terms  of  the  known  parameters.  ««,, 
/1 ,  and  /4 .  Referring  to  Fig.  6. 

«,  =tan  1  — -  (28) 

B,a, 

where  the  angle  u4  can  be  calculated  from  equation  (22). 


5.2  Motion  of  the  Top  Plate  With  Constraints  on  the  5.4  Rotatability  of  the  Ball  Joint  on  the  Top  Plate. 
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Fig.  8  Th#  boundaries  ol  workspace  on  th»  XZ  plane 


Referring  (o  Tig.  6,  (he  position  angles  ol  ihc  ball  joinis  on  ihc 
iop  plate,  i.e.,  the  /3's,  arc  the  angles  formed  by  each  of  the 
legs  with  the  top  plate.  We  have 

li  i  =  «|  +  <t> 

Ht  =ABSUu  -<*>)  (29) 

where  <t>  is  the  tilt  angle  of  the  top  plate.  The  position  angles 
and  lit,  however,  do  not  have  a  simple  relationship.  The 
derivation  of  them  first  requires  determining  the  distance 
between  the  ball  joint  B,  (or  B, )  and  the  plane  containing  the 
top  plate,  i.e.,  B,D2  (orfl,£M  (Fig.  1).  One  way  to  find  the 
distance  B, />.  is  to  analyze  the  geometrical  relationship  of  the 
project  of  the  triangle  B:b2D:  on  the  plane,  denoted  as  plane 
A,  which  is  parallel  to  the  XZ  plane  and  contains  the  line 
B;D The  distance  B, D:  can  therefore  be  obtained  from 

B:Di  =  B:b:'  sink  (30) 


Similarly,  the  equation  of  the  position  angle  of  the  ball  loim 
b,  can  be  obtained  as  follows: 


where 

B,D,-V(:4  *  <  (0.3rf'„  »  v4  4  0.5fl)'sinX 

Having  analytically  developed  the  necessary  kinematic 
relationship  of  this  special  6-SI’S  mechanism,  we  are  now  m  a 
position  to  develop  an  algorithm  to  numerically  outline  the 
reachable  boundary  of  the  mechanism  on  the  .V/  plane 


where 


B:b:‘  -  *  I  5rS):  +  (l.5rT+  xt  -0.5 R)'- 

and 

X=»-/tflsftan  1  |  ♦  rt 

l  hj£2  J 

The  position  angle  of  the  ball  joint  b:  can  then  be  calculated 
front 


6  A  Numerical  Algorithm 

An  algorithm  for  detecting  the  reachable  boundai  y  ot  the  6- 
SPS  mechanism  on  the  XZ.  plane  is  presented  as  follows: 

I.  ( a I  Scan  the  angle  <»,  within  the  physical  rotational 

limits  of  the  given  ball  joint, 

(6)  -Scan  the  lengths  /,  and  I,  within  the  specified 
extensible  limits  of  the  legs. 
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2.  For  a  particular  set  of  a, ,  ,  and  A,: 

(a)  -Check  the  limits  of  the  corresponding  leg  lengths 
l:  and  / 1 , 

( b )  -  Check  the  limits  of  the  corresponding  ball  joint 
angles,  i.e.,  a’s  and  fJ’s. 

If  out  of  the  range  then  go  to  Step  1,  otherwise  continue. 

3.  Calculate  the  desired  position  of  the  center  of  the  top 
plate,  i.e.,  ptxl.,zr)  and  record  the  extremes  on  a  data 
base.  Go  to  Step  I  until  all  scanning  processes  are 
complete. 

4.  Outline  the  workspace  boundary  from  the  recorded 
extremes. 

A  computer  program  which  is  based  on  the  analytical 
formulations  and  criteria  derived  in  this  investigation,  and 
which  is  written  in  FORTRAN,  is  developed.  Given  a  6-SPS 
platform  with  evenly  distributed  ball  joints,  the  program 
outlines  the  boundary  of  the  workspace  on  the  XZ  plane  of 
the  center  point  of  the  top  plate.  Inputs  of  this  program  are: 
the  width  of  the  hold  of  the  ball  joints  (e),  the  radius  of  the 
ball  head  (/<*),  the  diameter  of  the  connecting  legs  (<•/),  the 
extensible  limits  of  the  connecting  legs  (/lnal  and  /„„„),  the 
radii  of  the  top  and  the  bottom  plates  (r  and  K),  and  the  leg 
lengths/, ,  /4.  The  scanning  interval  of  the  joint  angle  t«,  needs 
to  be  specified. 

Four  practical  examples  are  given  to  demonstrate  the  range 
of  the  workspace  and  the  maneuverability  of  this  type  of 
mechanisms. 

Case  I :  A  platform  mechanism  has  the  following  design 
parameters:  /mjx  =  15”,  /mm  =  10”,  r= 2,  K  =  5”,  and  the  ball 
joint  dimensions  are  e  =  0.i ",  Rh  =0.5",  and  d  =  0.2" .  Using 
the  algorithm,  the  cross  section  of  the  workspace  on  the  XZ 
plane  is  shown  in  Fig.  8(<r). 

Case  2:  All  design  parameters  remain  the  same  as  Case  I 
except  the  width  of  the  hold  and  the  diameter  of  the  con¬ 
necting  leg  of  the  ball  joint  are  reduced  to  half  of  the  original 
values,  i.e.,  e  =  0.05"  and  r/  =  0.I '.  The  cross  section  of  the 
workspace  on  the  XZ  plane  is  given  in  Fig.  8(h). 

Case  3\  All  design  parameters  are  the  same  as  Case  2 
except  the  radius  of  the  top  plate  is  reduced  to  half  of  its 
original  value,  i.e.,  r=  I  ",  The  cross  section  of  the  workspace 
on  the  XZ  plane  is  given  in  Fig.  8(c). 

Case  4\  Repeat  the  same  design  parameters  as  Case  I 
except  the  radius  of  the  top  plate  is  now  reduced  to  half  of  its 
original  si/e,  i.e,.  r-  1 ".  The  cross-section  of  the  workspace 
on  the  XZ  plane  is  given  in  Fig.  8(r /). 

7  Conclusion 

This  paper  provides  a  basic  kinematic  investigation  on 
platform  type  manipulators.  It  is  found  that  the  6-Sl’S,  or  the 


Steward  platform  appears  to  be  the  only  mechanism  of  its 
type  that  can  possibly  be  adopted  as  a  general  maneuverable 
device.  The  analytical  model  and  inverse  control  equations  are 
derived.  Although  the  complexity  of  the  model,  and  the 
solvability  of  the  inverse  equations  are  more  favorable  as 
compared  to  a  six-joint  articulated  robot,  the  workspace  and 
the  maneuverability  seem  to  be  relatively  restricted.  This  is 
expected  and  due  mainly  to  the  multiloop  structure  and  the 
interference  between  the  ball  joints,  legs,  and  plates.  Unlike 
some  previous  studies  [5,  6)  in  which  ball-and-socket  joints  of 
the  platform  mechanism  were  replaced  by  two  revolute  joints 
in  their  analytical  model,  this  investigation  considers  the 
rotatability  of  the  ball  joints  and  its  formulation  is  believed  to 
be  of  new  and  practical  significance.  The  result  of  the  in¬ 
vestigation  on  workspace,  performed  in  a  special  cave, 
represents  an  upper  bound  for  our  design  considerations. 

It  should  be  noted,  however,  that  this  investigation  only 
concerns  a  kinematic  viewpoint.  An  investigation  of  the 
statics  and  dynamics  of  this  mechanism,  which  are  currently 
underway,  is  essential  for  the  practical  application  of  this 
platform.  The  result  of  this  research,  it  is  hoped,  will  con¬ 
tribute  toward  a  basic  understanding  of  the  limitations  as  well 
as  the  potential  usefulness  of  this  type  of  manipulator. 
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